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Abstract 

In the past ten years, the ideas of supersymmetry have been prof- 
itably applied to many nonrelativistic quantum mechanical problems. 
In particular, there is now a much deeper understanding of why cer- 
tain potentials are analytically solvable and an array of powerful new 
approximation methods for handling potentials which are not exactly 
solvable. In this report, we review the theoretical formulation of su- 
persymmetric quantum mechanics and discuss many applications. Ex- 
actly solvable potentials can be understood in terms of a few basic 
ideas which include supersymmetric partner potentials, shape invari- 
ance and operator transformations. Familiar solvable potentials all 
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have the property of shape invariance. We describe new exactly solv- 
able shape invariant potentials which include the recently discovered 
self-similar potentials as a special case. The connection between in- 
verse scattering, isospectral potentials and supersymmetric quantum 
mechanics is discussed and multi-soliton solutions of the KdV equation 
are constructed. Approximation methods are also discussed within the 
framework of supersymmetric quantum mechanics and in particular 
it is shown that a supersymmetry inspired WKB approximation is 
exact for a class of shape invariant potentials. Supersymmetry ideas 
give particularly nice results for the tunneling rate in a double well 
potential and for improving large N expansions. We also discuss the 
problem of a charged Dirac particle in an external magnetic field and 
other potentials in terms of supersymmetric quantum mechanics. Fi- 
nally, we discuss structures more general than supersymmetric quan- 
tum mechanics such as parasupersymmetric quantum mechanics in 
which there is a symmetry between a boson and a para-fermion of 
order p. 
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1 Introduction 



Physicists have long strived to obtain a unified description of all basic inter- 
actions of nature, i.e. strong, electrowea k, and gravitational interactions. 
Several ambitious attempts have been made in the last two decades, and it is 
now widely felt that supersymmetry (SUSY) is a necessary ingredient in any 
unifying approach. SUSY relates bosonic and fermionic degrees of freedom 
and has the virtue of taming ultraviolet divergences. It was discovered in 
1971 by Gel'fand and Likhtman jl], Ramond and Neveu and Schwartz 
U and later was rediscovered by several groups f|, ||, |J. The algebra in- 
volved in SUSY is a graded Lie algebra which closes under a combination of 
commutation and anti-commutation relations. It was first introduced in the 
context of the string models to unify the bosonic and the fermionic sectors. 
It was later shown by Wess and Zumino how to construct a 3 + 1 dimen- 
sional field theory which was invariant under this symmetry and which had 
very interesting properties such as a softening of ultraviolet divergences as 
well as having paired fermionic and bosonic degrees of freedom. For particle 
theorists, SUSY offered a possible way of unifying space-time and internal 
symmetries of the S-matrix which avoided the no-go theorem of Coleman and 
Mandula []7J which was based on the assumption of a Lie algebraic realization 
of symmetries (graded Lie algebras being unfamiliar to particle theorists at 
the time of the proof of the no-go theorem). Gravity was generalized by 
incorporating SUSY to a theory called supergravity || ||. In such theories, 
Einstein's general theory of relativity turns out to be a necessary consequence 
of a local gauged SUSY. Thus, local SUSY theories provide a natural frame- 
work for the unification of gravity with the other fundamental interactions 
of nature. 

Despite the beauty of all these unified theories, there has so far been 
no experimental evidence of SUSY being realized in nature. One of the 
important predictions of unbroken SUSY theories is the existence of SUSY 
partners of quarks, leptons and gauge bosons which have the same masses 
as their SUSY counterparts. The fact that no such particles have been seen 
implies that SUSY must be spontaneously broken. One hopes that the scale 
of this breaking is of the order of the electroweak scale of 100 GeV in order 
that it can explain the hierarchy problem of mass differences. This leads 
to a conceptual problem since the natural scale of symmetry breaking is 
the gravitational or Planck scale which is of the order of 10 19 GeV. Various 
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schemes have been invented to try to resolve the hierarchy problem, including 
the idea of non-perturbative breaking of SUSY. It was in the context of this 
question that SUSY was first studied in the simplest case of SUSY quantum 
mechanics (SUSY QM) by Witten and Cooper and Freedman |nj]. In 



a subsequent paper, a topological index was introduced (the Witten index) 
by Witten |L2| for studying SUSY breaking and several people studied the 
possibility that instantons provide the non-perturbative mechanism for SUSY 
breaking. In the work of Bender et al. [13] a new critical index was introduced 
to study non-perturbatively the breakdown of SUSY in a lattice regulated 
theory. Thus, in the early days, SUSY was studied in quantum mechanics as 
a testing ground for the non-perturbative methods of seeing SUSY breaking 
in field theory. 

Once people started studying various aspects of SUSY QM, it was soon 
clear that this field was interesting in its own right, not just as a model 
for testing field theory methods. It was realized that SUSY gives insight 
into the factorization method of Infeld and Hull JE4J which was the first 
method to categorize the analytically solvable potential problems. Gradually 
a whole technology was evolved based on SUSY to understand the solvable 
potential problems. One purpose of this article is to review some of the major 
developments in this 

Before we present a brief historical development of supersymmetric quan- 
tum mechanics, let us note another remarkable aspect. Over the last 10 
years, the ideas of SUSY have stimulated new approaches to other branches 



of physics [15]. For example, evidence has been found for a dynamical SUSY 
relating even-even and even-odd nuclei. The Langevin equation and the 
method of stochastic quantization has a path integral formulation which em- 
bodies SUSY. There have also been applications of SUSY in atomic, con- 
densed matter and statistical physics |T3| . 

In SUSY QM one is considering a simple realization of a SUSY algebra 
involving the fermionic and the bosonic operators. Because of the existence 
of the fermionic operators which commute with the Hamiltonian, one obtains 
specific relationships between the energy eigenvalues, the eigenfunctions and 
the S-matrices of the component parts of the full SUSY Hamiltonian. These 
relationships will be exploited in this article to categorize analytically solv- 
able potential problems. Once the algebraic structure is understood, the 
results follow and one never needs to return to the origin of the Fermi-Bose 
symmetry. In any case, the interpretation of SUSY QM as a degenerate Wess- 
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Zumino field theory in one dimension has not led to any further insights into 
the workings of SUSY QM. 

The introduction by Witten |T2] of a topological invariant to study dy- 
namical SUSY breaking led to a flurry of interest in the topological aspects of 
SUSY QM. Various properties of the Witten index were studied in SUSY QM 
and it was shown that in theories with discrete as well as continuous spectra, 
the index could display anomalous behavior 



16, 17, 18, 19, EOL 21, 22 



Using the Wigner-Kirkwood h expansion it was shown that for systems in 
one and two dimensions, the first term in the h expansion gives the exact 
Witten index [^3J . Further, using the methods of SUSY QM, a derivation of 
the Atiyah-Singer index theorem was also given 



25 



|, £7], P8| . In another 
development, the relationship between SUSY and the stochastic differential 
equations such as the Langevin equation was elucidated and exploited by 
Parisi and Sourlas [215 and Cooper and Freedman [IT|. This connection, 
which implicitly existed between the Fokker-Planck equation, the path inte- 
grals and the Langevin equation was then used to prove algorithms about 
the stochastic quantization as well as to solve non-perturbatively for the 
correlation functions of SUSY QM potentials using the Langevin equation. 

A path integral formulation of SUSY QM was first given by Salomonson 
and van Holten |30j . Soon afterwards it was shown by using SUSY methods, 
that the tunneling rate through double well barriers could be accurately 
determined 



m 
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At the same time, several workers extended 
ideas of SUSY QM to higher dimensionsal systems as well as to systems with 
large numbers of particles with a motivation to understand the potential 
problems of widespread interest in nuclear, atomic, statistical and condensed 
matter physics |35|, |36|, |37|, ^ 



91 TO, 41, H 



In 1983, the concept of a shape invariant potential (SIP) within the struc- 
ture of SUSY QM was introduced by Gendenshtein jPJ . This Russian paper 
remained largely unnoticed for several years. A potential is said to be shape 
invariant if its SUSY partner potential has the same spatial dependence as the 
original potential with possibly altered parameters. It is readily shown that 
for any SIP, the energy eigenvalue spectra could be obtained algebraically 
Much later, a list of SIPs was given and it was shown that the energy 



eigenfunctions as well as the scattering matrix could also be obtained alge- 
braically for these potentials [[4^, [46], [47], HH[. It was soon realized that the 



formalism of SUSY QM plus shape invariance (connected with translations of 
parameters) was intimately connected to the factorization method of Infeld 
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and Hull |TJ . 

It is perhaps appropriate at this point to digress a bit and talk about the 
history of the factorization method. The factorization method was first intro- 
duced by Schrodinger j|9 to solve the hydrogen atom problem algebraically. 
Subsequently, Infeld and Hull |14| generalized this method and obtained a 
wide class of solvable potentials by considering six different forms of factor- 
ization. It turns out that the factorization method as well as the methods 
of SUSY QM including the concept of shape invariance (with translation 
of parameters), are both reformulations f50| of Riccati's idea of using the 
equivalence between the solutions of the Riccati equation and a related sec- 
ond order linear differential equation. This method was supposedly used for 
the first time by Bernoulli and the history is discussed in detail by Stahlhofen 
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The general problem of the classification of SIPs has not yet been solved. 
A partial classification of the SIPs involving a translation of parameters was 
done by Cooper et al. |52], |53|. It turns out that in this case one only gets 
all the standard analytically solvable potentials contained in the list given by 
Dutt et al. |>4j] except for one which was later pointed out by Levai j5l|. The 



connection between SUSY, shape invariance and solvable potentials j56, 5? 
is also discussed in the paper of Cooper et al. [[5^] where these authors 



show that shape invariance even though sufficient, is not necessary for exact 
solvability. Recently, a new class of SIPs has been discovered which involves 
a scaling of parameters [58 ] . These new potentials as well as multi-step SIPs 
[dE] have been studied, and their connection with self-similar potentials as 
well as with q-deformations has been explored |60], |6l], |63|| . 



In yet another development, several people showed that SUSY QM offers 
a simple way of obtaining isospectral potentials by using either the Dar- 
boux or Abraham-Moses |65| or Pursey [^6| techniques, thereby offering 
glimpses of the deep connection between the methods of the inverse quantum 
scattering [§7], and SUSY QM |7D|, |7TJ, |72j|. The intimate connection 

between the soliton solutions of the KdV hierarchy and SUSY QM was also 
brought out at this time J7j| [TJ, [75] [TjJ. 

Approximate methods based on SUSY QM have also been developed. 
Three of the notable ones are the 1/N expansion within SUSY QM [77], 5 ex- 
pansion for the superpotential [[TJ| and a SUSY inspired WKB approximation 
(SWKB) in quantum mechanics for the case of unbroken SUSY []79|, ^(J jSljl . 



It turns out that the SWKB approximation preserves the exact SUSY rela- 
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tions between the energy eigenvalues as well as the scattering amplitudes of 
the partner potentials f82fl . Further, it is not only exact for large n (as any 
WKB approximation is) but by construction it is also exact for the ground 
state of V\{x). Besides it has been proved that the lowest order SWKB 
approximation necessarily gives the exact spectra for all SIP (with transla- 
tion). Subsequently a systematic higher order SWKB expansion has been 
developed and it has been explicitly shown that to 0{U & ) all the higher order 
corrections are zero for these SIP f83| . This has subsequently been gener- 
alised to all orders in h |]84|, 85]. Energy eigenvalue spectrum has also been 
obtained for several non-SIP |86|, ^7] [88], ^ and it turns out that in many 
of the cases the SWKB does better than the usual WKB approximation. 
Based on a study of these and other examples, it has been suggested that 
shape invariance is not only sufficient but perhaps necessary for the lowest 

9T|. Some attempts have 



order SWKB to give the exact bound state spectra 
also been made to obtain the bound state eigenfunctions within the SWKB 
formalism plj 92] , |93f| . 

Recently, Inomata and Junker pil have derived the lowest order SWKB 



quantization condition (BSWKB) in case SUSY is broken. It has recently 
been shown that for the cases of shape invariant three dimensional oscillator 
as well as for Poschl- Teller I and II potentials with broken SUSY, this lowest 
order BSWKB calculation gives the exact spectrum [^, 0. Recently, Dutt 
et al. have also developed a systematic higher order BSWKB expansion 



and using it have shown that in all the three (shape invariant) cases, the 
higher order corrections to 0(h e ) are zero. Further, the energy eigenvalue 
spectrum has also been obtained in the case of several non-SIP and it turns 
out that in many cases BSWKB does as well as (if not better than) the usual 
WKB approximation 



96 



Recursion relations between the propagators pertaining to the SUSY part- 
ner potentials have been obtained and explicit expressions for propagators of 
several SIP have been obtained . 

Several aspects of the Dirac equation have also been studied within SUSY 
QM formalism [IPO , |101| , |102| |. In particular, it has been shown using the re- 
sults of SUSY QM and shape invariance that whenever there is an analytically 
solvable Schrodinger problem in 1-dimensional QM then there always exists 
a corresponding Dirac problem with scalar interaction in 1+1 dimensions 
which is also analytically solvable. Further, it has been shown that there is 
always SUSY for massless Dirac equation in two as well as in four Euclidean 
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dimensions. The celebrated problem of the Dirac particle in a Coulomb field 
has also been solved algebraically by using the concepts of SUSY and shape 
invariance [103|. The SUSY of the Dirac electron in the field of a magnetic 
monopole has also been studied [104, |105||. Also, the classic calculation of 



Schwinger on pair production from strong fields can be dramatically simpli- 
fied by exploiting SUSY. 

The formalism of SUSY QM has also been recently extended and models 
for parasupersymmetric QM ||106| , |107| , |108|| as well as orthosupersymmetric 
QM |109| of arbitrary order have been written down. The question of singular 
superpotentials has also been discussed in some detail within SUSY QM 
formalism ||110| , |1 1 1| , |1 12| , |113| , 114] . Very recently it has been shown that 
SUSY QM offers a systematic method | |115| | for constructing bound states in 
the continuum [IT?|, |TT8| . 



As is clear from this (subjective) review of the field, several aspects of 
SUSY QM have been explored in great detail in the last ten years and it is 
almost impossible to cover all these topics and do proper justice to them. We 
have therefore, decided not to pretend to be objective but cover only those 
topics which we believe to be important and which we believe have not so 
far been discussed in great detail in other review articles. We have, however, 
included in Sec. 14 a list of the important topics missed in this review and 
given some references so that the intersested reader can trace back and study 
these topics further. We have been fortunate in the sense that review articles 
already exist in this field where several of these missing topics have been 

3, 0, [123], |T2|, [125], |10g]. We must also apologize 



discussed [119. 120, 121 



to several authors whose work may not have been adequately quoted in this 
review article in spite of our best attempts. 

The plan of the article is the following: In Sec. 2, we discuss the Hamil- 
tonian formalism of SUSY QM. We have deliberately kept this section at a 
pedagogical level so that a graduate student should be able to understand 
and work out all the essential details. The SUSY algebra is given and the 
connection between the energy eigenvalues, the eigenfunctions and the S- 
matrics of the two SUSY partner Hamiltonians are derived. The question of 
unbroken vs. broken SUSY is also introduced at a pedagogical level using 
polynomial potentials of different parity and the essential ideas of partner po- 
tentials are illustrated using the example of a one dimensional infinite square 
well. The ideas of SUSY are made more explicit through the example of one 
dimensional SUSY harmonic oscillator. 
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In Sec. 3, we discuss the connection between SUSY QM and factorization 
and show how one can always construct a hierarchy of p > 1 Hamiltonians 
with known energy eigenvalues, eigenf unctions and S-matrices by starting 
from any given Hamiltonian with p bound states whose energy eigenvalues, 
eigenfunctions and S-matrices (reflection coefficient R and transmission co- 
efficient T) are known. 

Sec. 4 is in a sense the heart of the article. We first show that if the 
SUSY partner potentials satisfy an integrability condition called shape in- 
variance then the energy eigenvalues, the eigenfunctions and the S-matrices 
for these potentials can be obtained algebraically. We then discuss satisfying 
the condition of shape invariance with translations and show that in this 
case the classification of SIP can be done and the resultant list of solvable 
potentials include essentially all the popular ones that are included in the 
standard QM textbooks. Further, we discuss the newly discovered one and 
multi-step shape invariant potentials when the partner potentials are related 
by a change of parameter of a scaling rather than translation type. It turns 
out that in most of these cases the resultant potentials are reflectionless and 
contain an infinte number of bound states. Explicit expressions for the energy 
eigenvalues, the eigenfunctions and the transmission coefficients are obtained 
in various cases. It is further shown that the recently discovered self-similar 
potentials which also statisfy g-SUSY, constitute a special case of the SIP. 
Finally, we show that a wide class of noncentral but separable potential prob- 
lems are also algebraically solvable by using the results obtained for the SIP 
[ |131|| . As a by product, exact solutions of a number of three body problems 
in one dimension are obtained analytically ||132|| . 

In Sec. 5, we discuss the solvable but non-shape invariant Natanzon and 
Ginocchio potentials and show that using the ideas of SUSY QM, shape 
invariance and operator transformations, their spectrum can be obtained 
algebraically. We also show that the Natanzon potentials are not the most 
general solvable potentials in nonrelativistic QM. 

Sec. 6 is devoted to a discussion of the SUSY inspired WKB approxima- 
tion (SWKB) in quantum mechanics both when SUSY is unbroken and when 
it is spontaneously broken. In the unbroken case, we first develop a system- 
atic higher order h -expansion for the energy eigenvalues and then show that 
for the SIP with translation, the lowest order term in the ^-expansion gives 
the exact bound state spectrum. We also show here that even for many of 
the non-SIP, the SWKB does as well as if not better than the WKB approx- 
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imation. We then discuss the broken SUSY case and in that case too we 
develop a systematic ^-expansion (BSWKB) for the energy eigenvalues. We 
show that even in the broken case the lowest order BSWKB gives the exact 
bound state spectrum for the SIP with translation. 

Sec. 7 contains a description of how SUSY QM can be used to con- 
struct multiparameter families of isospectral and strictly isospectral poten- 
tials. As an illustration we give plots of the one continuous parameter family 
of isospectral potentials corresponding to the one- dimensional harmonic os- 
cillator. From here we are immediately able to construct the two as well as 
multisoliton solutions of the KdV equation. 

In Sec. 8, we discuss more formal aspects of SUSY QM. In particular, we 
discuss the path integral formulation of SUSY QM as well as various sub- 
tleties associated with the Witten index [|12| A = (— 1) F . We also discuss 
in some detail the connection of SUSY QM with classical stochastic pro- 
cesses and discuss how one can develop a systematic strong coupling and S 
-expansion for the Langevin equation. 

Sec. 9 contains a description of several approximation schemes like the 
variational method, the 5-expansion, large- N expansion, energy splitting in 
double well potentials within the SUSY QM framework. 

In Sec. 10, we discuss the question of SUSY QM in higher dimensions. 
In particular, we discuss the important problem of a charged particle in a 
magnetic field (Pauli equation) in two dimensions and show that there is 
always a SUSY in the problem so long as the gyromagnetic ratio is 2. 

In Sec. 11, we show that there is always a SUSY in the case of massless 
Dirac equation in two or four Euclidean dimensions in the background of 
external electromagnetic fields. Using the results of SUSY QM we then list 
a number of problems with nonuniform magnetic field which can be solved 
analytically. We also show here that whenever a Schrodinger problem in 1- 
dimensional QM is analytically solvable, then one can always obtain an exact 
solution of a corresponding Dirac problem with the scalar coupling. We also 
show how the calculation of the fermion propagator in an external field can 
be simplified by exploiting SUSY. 

Sec. 12 contains a comprehensive discussion of the general problem of 
singular superpotentials, explicit breaking of SUSY, negative energy states 
and unpaired positive energy eigenstates. We also show here how to construct 
bound states in the continuum within the formalism of SUSY QM. 

Quantum mechanical models relating bosons and parafermions of order 
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p are described in Sec. 13. It is shown that such models encompass p su- 
persymmetries. Various consequences of such models are discussed including 
the connection with the hierarchy of Hamiltonians as well as with strictly 
isospectral potentials. We also discuss a quantum mechanical model where 
instead there is a symmetry between a boson and an orthofermion of order 
p. 

Finally, in Sec. 14, we give a list of topics related to SUSY QM which we 
have not discussed and provide some reference s for each of these topics. 

2 Hamiltonian Formulation of Supersymmet- 
ric Quantum Mechanics 

One of the key ingredients in solving exactly for the spectrum of one dimen- 
sional potential problems is the connection between the bound state wave 
functions and the potential. It is not usually appreciated that once one 
knows the ground state wave function (or any other bound state wave func- 
tion) then one knows exactly the potential (up to a constant). Let us choose 
the ground state energy for the moment to be zero. Then one has from the 
Schrodinger equation that the ground state wave function ipo( x ) obeys 

HlMx) = ~t^d^ + = ' W 

so that 

K ' 2m^ (i) 

This allows a global reconstruction of the potential V\(x) from a knowledge 
of its ground state wave function which has no nodes (we will discuss the 
case of using the excited wave functions later in Sec. 12). Once we realize 
this, it is now very simple to factorize the Hamiltonian using the following 
ansatz: 

ifi = AU (3) 



where 



A = JLA + W(x) , A* = - A A + W(x) . (4) 
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This allows us to identify 

V 1 (x) = W 2 (x) - -?L=W'(x) . (5) 

This equation is the well-known Riccati equation. The quantity W(x) is 
generally referred to as the "superpotential" in SUSY QM literature. The 
solution for W(x) in terms of the ground state wave function is 

W(x) = . (6) 

This solution is obtained by recognizing that once we satisfy Ai[) = 0, we 
automatically have a solution to Hitp = A^Aip = 0. 

The next step in constructing the SUSY theory related to the original 
Hamiltonian Hi is to define the operator H2 = AA^ obtained by reversing 
the order of A and A* . A little simplification shows that the operator H2 is 
in fact a Hamiltonian corresponding to a new potential V 2 (x). 

The potentials Vi(x) and V 2 (x) are known as supersymmetric partner poten- 
tials. 

As we shall see, the energy eigenvalues, the wave functions and the S- 
matrices of Hi and H 2 are related. To that end notice that the energy 
eigenvalues of both Hi and H 2 are positive semi-definite {E^^ > 0) . For 
n > 0, the Schrodinger equation for Hi 

Hi^ = A t Aij W = £(ity(D (8) 

implies 

H 2 (A^) = AA^A^ = E^(A^) . (9) 
Similarly, the Schrodinger equation for H 2 

H 2 ^ = AAW = ESH® (10) 

implies 

HiiAtyV) = AtAAtyW = E^(A^) . (11) 
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Prom eqs. (|J)-(0) and the fact that Eq = 0, it is clear that the eigenvalues 
and eigenfunctions of the two Hamiltonians H% and H 2 are related by (n = 
0,1,2,...) 



E 



(i) 

n+l) 



2# } 



0. 



Notice that if V'i+i 



(12) 

[E^^l, (13) 

[ E W)- 1/2 AW- (14) 
?/4 2 ^) of Hi (H 2 ) is normalized then the wave function 

tfj^ (ipn+i) m ec L s - dH) an d (0) ^ a l so normalized. Further, the operato r A 
(A^) not only converts an eigenfunction of Hi (H2) into an eigenfunction of 
H2(Hi) with the same energy, but it also destroys (creates) an extra node in 
the eigenfunction. Since the ground state wave function of Hi is annihilated 
by the operator A, this state has no SUSY partner. Thus the picture we get is 
that knowing all the eigenfunctions of Hi we can determine the eigenfunctions 
of H 2 using the operator A, and vice versa using we can reconstruct all 
the eigenfunctions of Hi from those of H 2 except for the ground state. This 
is illustrated in Fig. 2.1. 

The underlying reason for the degeneracy of the spectra of Hi and H 2 can 
be understood most easily from the properties of the SUSY algebra. That is 
we can consider a matrix SUSY Hamiltonian of the form 



H 



Hi 
H 2 



(15) 



which contains both Hi and H 2 . This matrix Hamiltonian is part of a closed 
algebra which contains both bosonic and fermionic operators with commu- 
tation and anti-commutation relations. We consider the operators 



Q 



"0 


0" 




A 


0_ 







AU 











(16) 



(17) 



in conjunction with H. The following commutation and anticommutation 
relation s then describe the closed superalgebra s/(l/l): 
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[H,Q] = [H,Q*}=0, 
{g,gt} = H, {Q,Q} = {Q\Q^} = . (18) 

The fact that the supercharges Q and commute with H is responsibl 
e for the degeneracy. The operators Q and can be interpreted as oper- 
ators which change bosonic degrees of freedom into fermionic ones and vice 
versa. This will be elaborated further below using the example of the SUSY 
harmonic oscillator. 

Let us look at a well known potential, namely the infinite square well and 
determine its SUSY partner potential. Consider a particle of mass m in an 
infinite square well potential of width L. 

V{x) = 0, < x < L , 

= oo, — oo < x < 0, x > L . (19) 

The ground state wave function is known to be 

4 1} = (2/L) 1/2 sin(7rs/L), < x < L , (20) 

ft 2 2 

and the ground state energy is E = 2m ^ a . 

Subtracting off the ground state energy so that we can factorize the 
Hamiltonian we have for Hi = H — E that the energy eigenvalues are 

«> = =££W (2!) 

and the eigenf unctions are 

^ = (2/L) 1 / 2 sin {n + 1)nX , < x < L . (22) 

The superpotential for this problem is readily obtained using eq. (^) 

W(x) = ^-cot(7rx/L) (23) 

V2m L 

and hence the supersymmetric partner potential V<i is 

h 2 n 2 

V 2 (x) = ^^[2cosec 2 (7rx/L) - 1] . (24) 
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The wave functions for H 2 are obtained by applying the operator A to the 
wave functions of H±. In particular we find that 

ipQ si 11 ( 7TX /L), ip\ oc sin(7ra;/L) sin(27ra;/L) . (25) 

Thus we have shown using SUSY that two rather different potentials 
corresponding to Hi and H 2 have exactly the same spectra except for the fact 
that H 2 has one fewer bound state. In Fig. 2.2 we show the supersymmetric 
partner potentials V\ and V 2 and the first few eigenfunctions. For convenience 
we have chosen L = n and % = 2m = 1. 

Supersymmetry also allows one to relate the reflection and transmission 
coefficients in situations where the two partner potentials have continuum 
spectra. In order for scattering to take place in both of the partner potentials, 
it is necessary that the potentials Vi,2 are finite as x — > —00 or as x — > +00 
or both. Define: 

W(x -> ±00) = W± . (26) 

Then 

Vi,2 -+Wl asx^ ±00. (27) 

Let us consider an incident plane wave e tkx of energy E coming from the 
direction x — > —00. As a result of scattering from the potentials Vi t2 (x) one 
would obtain transmitted waves T h2 (k)e ik ' x and reflected waves R 1<2 (k)e %kx . 
Thus we have 

V> (1,2) — > -00) -> e fc + R l)2 e~ ikx , 

^' 2 \k', x -> +00) -> T li2 e ifc,:c . (28) 

SUSY connects continuum wave functions of i?i and if 2 having the same 
energy analogously to what happens in the discrete spectrum. Thus we have 
the relationships: 

e ikx + R ie - ikx = N[{-ik + W_)e ikx + {ik + W.)e~ lkx R 2 ], 

Tie ik'x = N [(-ik' + W + )e ik ' x T 2 ], (29) 

where iV is an overall normalization constant. On equating terms with the 
same exponent and eliminating N, we find: 
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where k and k! are given by 

k = (E-Wl) 1 / 2 , k' = {E-WlY' 2 . (31) 

A few remarks are now in order at this stage. 

(1) Clearly |i?i| 2 = |_R 2 | 2 and | T\ | 2 = |T 2 | 2 , that is the partner potentials 
have identical reflection and transmission probabilities. 

(2) Ri(Ti) and R 2 (T 2 ) have the same poles in the complex plane except that 
i?i(7i) has an extra pole at k — —iW^. This pole is on the positive imaginary 
axis only if W- < in which case it corresponds to a zero energy bound state. 

(3) In the special case that W + = W_, we have that Ti(k) = T 2 (k). 

(4) When W. = then R±(k) = -R 2 (k). 

It is clear from these remarks that if one of the partner potentials is a con- 
stant potential (i.e. a free particle), then the other partner will be of necessity 
reflect ionless. In this way we can understand the reflect ionless potentials of 
the form V(x) = Asech 2 ax which play a critical role in understanding the 
soliton solutions of the KdV hierarchy. Let us consider the superpotential 

W(x) = A tanh ax . (32) 

The two partner potentials are 

V 1 = A 2 - A (A + a^=)sech 2 ax , 

V2m 

V 2 = A 2 - 

We see that for A = cc;^=, V 2 (x) corresponds to a constant potential so that 
the corresponding V 1 is a reflect ionless potential. It is worth noting that V 1 is 
^-dependent. One can in fact rigorously show,though it is not mentioned in 
most text books, that the reflectionless potentials are necessarily ^-dependent. 

So far we have discussed SUSY QM on the full line (— oo < x < oo). Many 
of these results have analogs for the n-dimensional potentials with spherical 
symmetry. For example, in three dimensions one can make a partial wave 
expansion in terms of the wave functions: 

*Pnim(r, 9, 0) = -R nl {r)Y lm (9, 0) . (34) 



A(A — a )sech .ax 



2m' 



(33) 
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Then it is easily shown [126|| that the reduced radial wave function R satisfies 



the one-dimensional Schrodinger equation (0 < r < oo) 

n 2 d 2 tfj(r) n , * zrz + m 2 .,,, 

- stsf 2 + [v(r) + h^-^ r) = Ei ' (r) (35) 

We notice that there is an effective one dimensional potential which contains 
the original potential plus an angular momentum barrier. The asymptotic 
form of the radial wave function for the I'th partial wave is 

^(rj) ^ ^[S\k')e ik ' r - (-l) l e- ik ' r \ , (36) 

where S l is the scattering function for the I'th. partial wave. i.e. S l (k) = e lSl ^ k ' 
and S is the phase shift. 

For this case we find the relations: 

= *W ■ < 37 > 

Here W + = W(r — > oo). 

We thus have seen that when H% contained a known ground state wave 
function then we could factorize the Hamiltonian and find a SUSY partner 
Hamiltonian H 2 - Now let us consider the converse problem. Suppose we are 
given a superpotential W(x). In this case there are two possibilities. The 
candidate ground state wave function is the ground state for Hi (or H 2 ) and 
can be obtained from: 

A^\x) = ^^ 1 \x)=Nexp^-y^J X W(y)dy^ , 

A^i 2 \x) = =^4 2) (^) =Nexp [+^^ f W (y) d v\ ■ ( 38 ) 

By convention, we shall always choose W in such a way that amongst 
Hi, if 2 only Hi (if at all) will have a normalizable zero energy ground state 
eigenfunction. This is ensured by choosing W such that W(x) is posi- 
tive(negative) for large positive (negative) x. This defines Hi to have fermion 
number zero in our later formal treatment of SUSY. 

If there are no normalizable solutions of this form, then Hi does not 
have a zero eigenvalue and SUSY is broken. Let us now be more precise. A 
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symmetry of the Hamiltonian (or Lagrangian) can be spontaneously broken 
if the lowest energy solution does not respect that symmetry, as for example 
in a ferromagnet, where rotational invariance of the Hamiltonian is broken 
by the ground state. We can define the ground state in our system by a two 
dimensional column vector: 



|0 >= if> (x) 



^\x) 



(2), 



X 



For SUSY to be unbroken requires 



Q\Q >= QT|Q >= 0|0 > 



(39) 



(40) 



Thus we have immediately from eq. (jl8| ) that the ground state energy 
must be zero in this case. For all the cases we discussed previously, the ground 
state energy was indeed zero and hence the ground state wave function for 
the matrix Hamiltonian can be written: 



4 1} (x) 



(41) 



where i()q~\x) is given by eq. (BS 



If we consider superpotentials of the form 

W(x) = gx n , 



(42) 



then for n odd and g positive one always has a normalizable ground state 
wave function. However for the case n even and g arbitrary, then there is 
no normalizable ground state wave function. In general when one has a 
superpotential W(x) so that neither Q nor annihilates the ground state 
as given by eq. (|3"§| ) then SUSY is broken and the potentials V\ and V 2 have 
degenerate positive ground state energies. Stated another way, if the ground 
state energy of the matrix Hamiltonian is non zero then SUSY is broken. 
For the case of broken SUSY the operators A and A* no longer change the 
number of nodes and there is a 1-1 pairing of all the eigenstates of Hi and 
i?2- The precise relations that one now obtains are: 



E^=E^>0, n = 0,1,2, 
^ = [^ ) ]- 1/2 M 1) , 



(43) 
(44) 
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€ ] = [E^\- 1/2 AH^. (45) 

while the relationship between the scattering amplitudes is still given by 
eqs. ( P0| ) or fl5T|). The breaking of SUSY can be described by a topological 



quantum number called the Witten index [T2| which we will discuss later. 
Let us however remember that in general if the sign of W{x) is opposite as 
we approach infinity from the positive and the negative sides, then SUSY is 
unbroken, whereas in the other cases it is always broken. 

2.1 State Space Structure of the SUSY Harmonic Os- 
cillator 

For the usual quantum mechanical harmonic oscillator one can introduce a 
Fock space of boson occupation numbers where we label the states by the 
occupation number n. To that effect one introduces instead of P and q the 
creation and annihilation operators a and a) . The usual harmonic oscillator 
Hamiltonian is 

Let us rescale the Hamiltonian in terms of dimensionless coordinates and 
momenta x and p so that we measure energy in units of hu). We put 

H = Hhu, q = (-^-) 1/2 x, P = (2mhu) 1/2 p. (47) 
2mw 



Then 

Now introduce 
Then 



H=(p 2 + ^), [x,p)=i. (48) 
J = (| + ip), a f = (|-zp). (49) 



[a, aJ] = 1, [N,a] = —a, [N, a)} = a\ 

N = a f a, H = N+^ . (50) 

The usual operator formalism for solving the harmonic oscillator potential 
is to define the ground state by requiring 

a|0>=0, (51) 
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which leads to a first order differential equation for the ground state wave 
function. The n particle state (which is the n'th excited wave function in the 
coordinate representation) is then given by: 

\n b >= -7—10 > ,. (52) 

where we have used the subscript b to refer to the boson sector as distinct from 
the fermions we will introduce below. For the case of the SUSY harmonic 
oscillator one can rewrite the operators Q (Q^) as a product of the bosonic 
operator a and the fermionic operator ip. Namely we write Q = aip^ and 
= afip where the matrix fermionic creation and annihilation operators 
are defined via: 

(53) 
(54) 

Thus ip and ift obey the usual algebra of the fermionic creation and 
annihilation operators, namely, 

{^} = 1, {^,^} = {^} = 0, (55) 

as well as obeying the commutation relation: 

(56) 

The SUSY Hamiltonian can be rewritten in the form 

H = QQt + gtg = + _ 1 ^ ^t]. (57) 

The effect of the last term is to remove the zero point energy. 

The state vector can be thought of as a matrix in the Schrodinger picture 
or as the state \n b ,rif > in this Fock space picture. Since the fermionic 
creation and annihilation operators obey anti-commutation relations hence 
the fermion number is either zero or one. As stated before, we will choose 



1 




— <7_ = 





1 



[*I>,V] = ^3 = 
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the ground state of H\ to have zero fermion number. Then we can introduce 
the fermion number operator 

nr = lzp.= l -V>+<\ . (58) 

Because of the anticommutation relation, rif can only take on the values 
and 1. The action of the operators a, aJ, ip, ^ in this Fock space are then: 

a\rib,rif > = \n b — l,n/>, t/>|ti&, rif >— |n&, n/ — 1 >, 
adrift, rif > = |n& + l,n/>, ^\n b , rif >= \n b , rif + 1 > . (59) 

We now see that the operator = —iaip^ has the property of changing a 
boson into a fermion without changing the energy of the state. This is the 
boson-fermion degeneracy characteristic of all SUSY theories. 

For the general case of SUSY QM, the operator a gets replaced by A 
in the definition of Q, Q\ i.e. one writes Q = Atp* and = Aty. The 
effect of Q and are now to relate the wave functions of Hi and H 2 which 
have fermion number zero and one respectively but now there is no simple 
Fock space description in the bosonic sector because the interactions are non- 
linear. Thus in the general case, we can rewrite the SUSY Hamiltonian in 
the form 

H = (-^ + W 2 )I - [V, ifi\W'. (60) 

This form will be useful later when we discuss the Lagrangian formulation 
of SUSY QM in Sec. 8. 

2.2 Broken Supersymmetry 

As discussed earlier, for SUSY to be a good symmetry, the operators Q 
and must annihilate the vacuum. Thus the ground state energy of the 
super-Hamiltonian must be zero since 

H = {Q\Q}. 



Witten |12j proposed an index to determine whether SUSY is broken in 



super symmetric field theories. The index is defined by 

A = Tr(-1) F , (61) 
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where the trace is over all the bound states and continuum states of the 
super-Hamiltonian. For SUSY QM, the fermion number rip = F is defined 
by \\\ — C3] and we can represent (— 1) F by the matrix a 3 . If we write the 
eigenstates of H as the vector: 

ipn(x) = 

then the ± corresponds to the eigenvalues of (— 1) F being ±1. For our con- 
ventions the eigenvalue +1 corresponds to Hi and the bosonic sector and the 
eigenvalue —1 corresponds to H 2 and the fermionic sector. Since the bound 
states of Hi and H 2 are paired, except for the case of unbroken SUSY where 
there is an extra state in the bosonic sector with E = we expect for the 
quantum mechanics situation that A = for broken SUSY and A = 1 for 
unbroken SUSY. In the general field theory case, Witten gives arguments 
that in general the index measures N + (E = 0) — N_(E = 0) which is the 
difference A between the number of Bose states and Fermi states of zero 
energy. In field theories the Witten index needs to be regulated to be well 
defined so that one considers instead 

A(/3) = Tr(-l) F e-^ , (63) 
which for SUSY quantum mechanics becomes 

A(/3)=Tr[e- /3Hl -e-^ Ha ]. (64) 

In field theory it is quite hard to determine if SUSY is broken non- 
perturbatively, and thus SUSY quantum mechanics became a testing ground 
for finding different methods to understand non-perturbative SUSY break- 
ing. In the quantum mechanics case, the breakdown of SUSY is related to 
the question of whether there is a normalizable wave function solution to the 
equation Q\0 >= 1 > which implies 

if, (x) = Ne-f w(x)dx . (65) 

As we said before, if this candidate ground state wave function does not 
fall off fast enough at ±00 then Q does not annihilate the vacuum and SUSY 
is spontaneously broken. Let us show using a trivial calculation that for 



(62) 
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two simple polynomial potentials the Witten index does indeed provide the 
correct answer to the question of SUSY breaking. Let us consider 



Expanding the term proportional to 03 in the exponent and taking the trace 
we obtain 



We are interested in the regulated index as (5 tends to 0, so that practically 
we need to evaluate 



If we directly evaluate this integral for any potential of the form W(x) = 
gx 2n+1 , which leads to a normalizable ground state wave function, then all 
the integrals are gamma functions and we explictly obtain A = 1. If in- 
stead W(x) = gx 2n so that the candidate ground state wave function is not 
normalizable then the integrand becomes an odd function of x and therefore 
vanishes. Thus we see for these simple cases that the Witten index immedi- 
ately coincides with the direct method available in the quantum mechanics 
case. 

Next let us discuss a favorite type of regularization scheme for field 
theory- namely the heat kernel method. (Later we will discuss a path integral 
formulation for the regulated Witten Index). 

Following Akhoury and Comtet one defines the heat kernels K±(x, y; (3) 
which satisfy 




(66) 




(67) 




(68) 



These have the following eigenfunction representation: 




(69) 



K±(x,y;(3) =< y\e 



-f3H± 



X > 




n 



(70) 



In terms of the heat kernels one has 
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, or 



poo 

A(/3) = N + (E = 0) - AL(£ = 0) + / dEe^ E ( P+ (E) - p_(E)), (71) 

•I En 



where p± corresponds to the density of states. What Akhoury and Comtet 
were able to show, was that in cases when W(x) went to different constants at 
plus and minus infinity, then the density of states factors for the continuum 
did not cancel and that A(/3) could depend on (3 and be fractional at (3 — 0. 
We refer the interested reader to the original paper for further details. 

Another non-perturbative method for studying SUSY breaking in field 
theory is to explicitly break SUSY by placing the theory on a lattice and 
either evaluating the path integral numerically or via some lattice non weak- 
perturbative method such as the strong coupling (or high temperature) ex- 
pansion. This method was studied in detail in |T3| , |2l[| and we will summarize 
the results here. The basic idea is to introduce a new parameter, namely the 
lattice spacing a. This parameter explicitly breaks SUSY so that the ground 
state energy of the system will no longer be zero, even in the unbroken case. 
One hopes that as the lattice spacing is taken to zero then the ground state 
energy will go to zero as a power of the lattice spacing if SUSY is unbroken, 
that is we expect 



where 7 is a critical index which if greater than zero should be easy to 
measure in a Monte Carlo calculation. Measuring the ground state energy 
at two different lattice spacings one studies: 



in the limit a 1 and a — > 0. The case of broken SUSY is more difficult because 
then we expect 7 = which is a hard measurement to make numerically. 
In this latter case it is easier to directly measure the ground state energy 
and show that it remains non-zero as one takes the lattice spacing to zero. 
To see how this works in quantum mechanics one can do a lattice strong 
coupling expansion of the Langevin equation which allows one to determine 
the ground state wave function of H\ as we shall show later. 

For the superpotential W(x) = gx 3 , we expect to find a positive critical 
index since here the candidate ground state wave function is proportional to 




(72) 




(73) 
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e ^ 4 /4 anc [ j s normalizable. The ground state expectation values of x n for 
the Hamiltonian Hi can be determined by first solving the Langevin equation 

^ + W(x) = V (t) (74) 

and then averaging x(r)(t)) over Gaussian noise whose width is related to h. 
Since by the virial theorem 

E = 3g 2 < x e > -3g < x 2 > , 

knowing the corelation functions will allow us to calculate the ground state 
energy. We first put the Langevin equation on a time lattice (t n = na): 

e{x n - x n -i) + gx 3 = r} n , (75) 

where e = 1/a, which allows a solution by strong coupling expansion for large 
g. The result is 

*n = (^) 1/3 + (^iC 273 " Vn 1/3 ) + 0(, 2 ). (76) 

As we will demonstrate in Sec. 8, the quantum mechanical expectation values 
of < x n (t) > are the same as the noise averaged expection values of x n (rj) 

< x "(t) >= J [Dr]]P[r]}x n (r](t)). (77) 

On the lattice the path integral becomes a product of ordinary integrals 
which can be performed with: 

P[rj] = U ie - a ^/\[^. (78) 

V Z7T 

The ground state energy on the lattice regulated theory then has the form 

m 

E = ^zJ2CnZ n , (79) 

n=0 

where 

3 2 
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is a dimensionless length. The critical index can be determined from the 
logarithmic derivative of Eq with respect to z. Using Pade approximants to 
extrapolate the lattice series to small lattice spacing we found that 

E = ca 116 (80) 

verifying that SUSY is unbroken in the continuum limit. Using the same 
methods for the case W(x) = gx 2 /2 we were able to verify that the ground 
state energy was not zero as we took the continuum limit. After verifying the 
applicability of this method in SUSY QM, Bender et al. then successfully 
used this method to study non-perturbative SUSY breaking in Wess-Zumino 
models of field theory [13[ . 



3 Factorization and the Hierarchy of Hamil- 
tonians 

In the previous section we found that once we know the ground state wave 
function corresponding to a Hamiltonian Hi we can find the superpotential 
Wi(x) from eq. (|6]). The resulting operators A\ and A\ obtained from eq. 
(fD can be used to factorize Hamiltonian Hi. We also know that the ground 
state wave function of the partner Hamiltonian Hi i s determined from the 
first excited state of H\ via the application of the operator A\. This allows a 
refactorization of the second Hamiltonian in terms of W%. The partner of this 
refactorization is now another Hamiltonian H3. Each of the new Hamiltoni- 
ans has one fewer bound state, so that this process can be continued until 
the number of bound states is exhausted. Thus if one has an exactly solv- 
able potential problem for Hi, one can solve for th e energy eigenvalues and 
wave functions for the entire hierarchy of Hamiltonians created by repeated 
refactorizations. Conversely if we know the ground state wave functions for 
all the Hamiltonians in this hierarchy, we can reconstruct th e solutions of 
the original problem. Let us now be more specific. 

From the last section we have seen that if the ground state energy of a 
Hamiltonian Hi is zero then it can always be written in a factorizable form 
as a product of a pair of linear differential operators. It is then clear that if 
the ground state energy of a Hamiltonian Hi is Eq with eigenfunction ipQ 
then in view of eq. (0), it can always be written in the form (unless stated 
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otherwise, from now on we set H = 2m = 1 for simplicity): 



H, = A{Ai + 4 1] = + V 1 (x), (81) 



where 



Al = 1. + Wl(x) , 4 = _A. + Wi(x) , Wi(x) = (82) 

The SUSY partner Hamiltonian is then given by 

H 2 = A 1 A\ + = + V 2 (x), (83) 

where 

d 2 

V 2 (x) = Wl + W[ + 4 1} = + 2^ = 14 (x) - 2— ln4 X) - (84) 

We will introduce the notation that in E^\ n denotes the energy level 
and (to) refers to the m'th Hamiltonian H m . From Sec. 2, the energy eigen- 
values and eigenfunctions of the two Hamiltonians Hi and H 2 are related 

EH = E^\ ^ = (E^-E^r 1/2 A 1 ^ 1 . (85) 

Now starting from H 2 whose ground state energy is E^f 1 = e[^ one can 
similarly generate a third Hamiltonian H 3 as a SUSY partner of H 2 since we 
can write H 2 in the form: 

if 2 = A 1 A\ + = A\A 2 + E[ 1] , (86) 

where 



M = i + m{x) ' 4 = + m{x) ' ^ = -^sr ■ (87) 

Continuing in this manner we obtain 



ff 3 = A 2 4 + = --^ + V 3 (x), (88) 
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where 



V 3 (x) = WZ + Wl + E™ = V 2 (x)-2^-\n4 2) 

dx z 

d<2 ■ ' <.(i)„/.( 2 )> 



Furthermore 



Vl {x)-2— HWW)- (89) 



p(3) _ p(2) _ p(l) 

^ = (E { s - 4 2) )- i/2 ^i 2 i 
= (4-^!V i/2 (e-4V i/2 Mfi. 00) 

In this way, it is clear that if the original Hamiltonian Hi has p(> 1) 
bound states with eigenvalues and eigenfunctions ipW with < n < (p— 
1), then we can always generate a hierarchy of (p — 1) Hamiltonians H 2 , ...H p 
such that the m'th member of the hierarchy of Hamiltonians (H m ) has the 
same eigenvalue spectrum as H\ except that the first (to — 1) eigenvalues of 
H\ are missing in H m . In particular, we can always write (to = 2, 3, ...p): 

H m = A\ n A m + E$_ x = ~ + V m (x), (91) 



where 

A 

ax ax 

One also has 



d d lm// m) 

A m = f- + W m (x) , W m (*) = (92 ) 



p(m) p(m-l) __ _ 

J -'n n+1 •'• -"-^n+rre— 1 ) 



V m (s) = \/ 1 (x)-2^1n(4 1) ...4 m - 1) ). (93) 

In this way, knowing all the eigenvalues and eigenfunctions of Hi we im- 
mediately know all the energy eigenvalues and eigenfunctions of the hierarchy 
of p — 1 Hamiltonians. Further the reflection and transmission coefficients (or 
phase shifts) for the hierarchy of Hamiltonians can be obtained in terms of 
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Ri, T\ of the first Hamiltonian H\ by a repeated use of eq. ( p0|) . In particular 
we find 



= \ :.., i-i ::.„. , !"•(*•>• 




where and fc' are given by 

fc = [E - (W^) 2 ] 172 , jfc' = [£ - (W^) 2 ] 1 / 2 . (95) 

4 Shape Invariance and Solvable Potentials 

Most text books on quantum mechanics describe how the one dimensional 
harmonic oscillator problem can be elegantly solved using the raising and 
lowering operator method. Using the ideas of SUSY QM developed in Sec. 2 
and an integrability condition called the shape invariance condition |33 



we 



now show that the operator method for the harmonic oscillator can be gen- 
eralized to the whole class of shape invariant potentials (SIP) which include 
all the popular, analytically solvable potentials. Indeed, we shall see that for 
such potentials, the generalized operator method quickly yields all the bound 
state energy eigenvalues, eigenfunctions as well as the scattering matrix. It 
turns out that this approach is essentially equivalent to Schrodinger's method 
of factorization [[29], [IJJ although the language of SUSY is more appealing. 

Let us now explain precisely what one means by shape invariance. If the 
pair of SUSY partner potentials V\ i2 (x) defined in Sec. 2 are similar in shape 
and differ only in the parameters that appear in them, then they are said to 
be shape invariant. More precisely, if the partner potentials Vi t2 (x; ai) satisfy 
the condition 

V 2 (x;a 1 ) = V 1 (x;a 2 )+R(a 1 ), (96) 

where a% is a set of parameters, a 2 is a function of a\ (say a 2 = f{a>i)) and 
the remainder R(ai) is independent of x, then Vi(x; oi) and V 2 (x; a\) are said 



to be shape invariant. The shape invariance condition ( |9q) is an integrability 



condition. Using this condition and the hierarchy of Hamiltonians discussed 
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in Sec. 3, one can easily obtain the energy eigenvalues and eigenf unctions of 
any SIP when SUSY is unbroken. 



4.1 General Formulas for Bound State Spectrum, Wave 
Functions and S-Matrix 

Let us start from the SUSY partner Hamiltonians Hi and H2 whose eigenval- 
ues and eigenf unctions are related by SUSY. Further, since SUSY is unbroken 
we know that 



E, 



(i) 



en r f x 

= 0, % ; (x;ai) = iVexp - / Wi(y\ai)dy 



(97) 



We now show that the entire spectrum of Hi can be very easily obtained 
algebraically by using the shape invariance condition (0). To that purpose, 
let us construct a series of Hamiltonians H s , s = 1,2,3... In particular, 
following the discussion of the last section it is clear that if Hi has p bound 

states then one can construct p such Hamiltonians Hi, H2 H p and the n'th 

Hamiltonian H n will have the same spectrum as Hi except that the first n — 1 
levels of Hi will be absent in H n . On repeatedly using the shape invariance 
condition (|96|), it is then clear that 

d 2 s ~ 1 
H s = -— + Vi(x-a s ) + Y / R(a k ), (98) 
ax k=l 

where a s = f^ 1 ^) i.e. the function / applied s — 1 times. Let us compare 
the spectrum of H s and H s+ i. In view of eqs. (|96D and (^) we have 



d 2 s 
H s+ i = -— + Vl(x; a s+ i) + R ( a k) 



k=l 

d 2 ^ 



dx 2 



+ V 2 (x;a s ) + Y,R(a k ). (99) 

k=l 



Thus H s and H s+ i are SUSY partner Hamiltonians and hence have identical 
bound state spectra except for the ground state of H s whose energy is 

E^ = J2R(a k ). (100) 

k=l 
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This follows from eq. (j98| ) and the fact that Eq = 0. On going back from 
H s to i? s _i etc, we would eventually reach H 2 and Hi whose ground state 
energy is zero and whose n'th level is coincident with the ground state of the 
Hamiltonian H n . Hence the complete eigenvalue spectrum of Hi is given by 

n 

KM = E Riflh)\ E ( ai ) = o. (ioi) 

k=l 

We now show that, similar to the case of the one dimensional harmonic 
oscillator, the bound state wave functions ip£'(x\ ai) for any shape invariant 
potential can also be easily obtained from its ground state wave function 
i/j^\x;ai) which in turn is known in terms of the superpotential. This is 
possible because the operators A and link up the eigenf unctions of the 
same energy for the SUSY partner Hamiltonians Hi 2- Let us start from the 
Hamiltonian H s as given by eq. ( P3| ) whose ground state eigenfunction is 
then given by iPq\x\ a s ). On going from H s to H s _i to H 2 to Hi and using 
eq. (0) we then find that the n'th state unnormalized, energy eigenfunction 
ip^\x] ai) for the original Hamiltonian Hx(x; ai) is given by 

i)n\x; a\) oc A\x\ a{)A\x\ a 2 )...A\x\ a n )i[)o\ 

), (102) 

which is clearly a generalization of the operator method of constructing the 
energy eigenfunctions for the one dimensional harmonic oscillator. 

It is often convenient to have explicit expressions for the wave functions. 
In that case, instead of using the above equation, it is far simpler to use the 
identify |46| 

oi) = A\x- ax)^lx(x; a 2 ). (103) 

Finally, it is worth noting that in view of the shape invariance condition 
(p6|), the relation ( pOj) between scattering amplitudes takes a particularly 
simple form 

*<*^ = (S££) *<*"■>■ (104) 

T ^={Wr^) Tilki " ) ' (105) 

thereby relating the reflection and transmission coefficients of the same Hamil- 
tonian Hi at ax and a 2 (= f(ai)). 
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4.2 Shape Invariance in More Than One Step 

We can expand the list of solvable potentials by extending the shape invari- 
ance idea to the more general concept of shape invariance in two and even 
multi-steps. We shall see later that in this way we will be able to go much 
beyond the factorization method and obtain a huge class of new solvable 
potentials [ ^] . 

Consider the unbroken SUSY case of two superpotentials W(x; a\) and 
W(x; ai) such that V 2 (x; ai) and V\(x\ ai) are same up to an additive constant 
i.e. 

V 2 {x-a 1 )=V 1 {x-a l )+R{a l ) (106) 

or equivalent ly 

W 2 (x; ai) + W'(x; a x ) = W 2 (x; a x ) - W'(x; a x ) + R( ai ). (107) 
Shape invariance in two steps means that 

V 2 (x;a 1 ) = V 1 (x;a 2 )+R(a 1 ), (108) 

that is 

W 2 (x; ai) + W'{x; a{) = W 2 {x; a 2 ) - W'{x; a 2 ) + R{a x ). (109) 

We now show that when this condition holds, the energy eigenvalues and 
eigenfunctions of the potential Vi(x; ai) can be obtained algebraically. First 
of all, let us notice that unbroken SUSY implies zero energy ground states 
for the potentials V%(x] ai) and V%(x] a%): 

Ej 1) (a 1 ) = 0, ^ 1) (a 1 ) = 0. (110) 

The degeneracy of the energy levels for the SUSY partner potentials yields 

= E^\ ai ) = EZ(ai). (in) 

From eq. ( |106| ) it follows that 

Ei 2 \a 1 )=E^(a 1 )+R(a 1 ), (112) 

so that for n = 0, these two equations yield 

£{ 1} (ai) =fl(oi). (113) 
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Also, the shape invariance condition ( |108| ) yields 

E^(a 1 ) = E^(a 2 ) + R(a 1 ). (114) 
From the above equations one can then show that 

^(ai) = E^M + R(ax) + ^(ai). (115) 
On solving these questions recursively we obtain (n = 0, 1, 2, ...) 

n 

E$ = y £ l [R(a k )+R(a k )), (116) 



k=l 



e£+i = J2[ R ( a k)+R(a k )}+R(a n+1 ). (117) 



i — 

k=l 

We now show that, similar to the discussion of the last subsection, the 
bound state wave functions %l)£'(x\a\) can also be easily obtained in terms 
of the ground state wave functions ip^\x; a\) and $q(x; ai) which in turn 
are known in terms of the superpotentials W and W. In particular from eq. 
QIUBD it follows that 



ipn+i( x 'i K A\x; ai)ip^\x; a{) oc A\x; a\)^\x\ ai), (H8) 
while from eq. ( |108| ) we have 

$n+x{ x \ x A\x] ai)ijj^(x; ai) oc A\x; ai)ip^\x; a 2 ). (H9) 

Hence on combining the two equations we have the identity 

Tpn+ifa a i) A\x\ ai)A\x] a\)ij)£\x\ a 2 ). (120) 

Recursive application of the above identity yields 

win { x \ a i) oc [A j (x; ai)i f (x; ai)]...[A*(x\ a n )A j (x; a n )]ip^\x; a„ +1 ), (121) 

^2n+i( x i a i) [A j (x; ai)A f (x; ai)]...[A f (x; a^i^x; a n )]A f (x; a„ +1 )^ 1} (x; a n+ i), 

(122) 
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where we have used the fact that 



ipi~\x\ ai) oc A\x; ai)ip^\x; ax). (123) 



Finally, it is easily shown that the relation fl3~0|) between the scattering 
amplitudes takes a particularly simple form 

( WJ^)+ik \ f W-(ai)+ik \ 
\W-{ax)-ik) \W-(ai)-ik) 

7i(fe;ai) = — - — Ti fc;a 2 , 125 

\W_(aiJ-z/c/ \ M/_(ai) - «A; / 

thereby relating the reflection and transmission coefficients of the same Hamil- 
tonian at a\ and a 2 . 

It is clear that this procedure can be easily generalized and one can con- 
sider multi-step shape invariant potentials and in these cases too the spec- 
trum, the eigenfunctions and the scattering matrix can be obtained alge- 
braically. 



4.3 Strategies For Categorizing Shape Invariant Po- 
tentials 

Let us now discuss the interesting question of the classification of various 
solutions to the shape invariance condition fl96|) . This is clearly an important 
problem because once such a classification is available, then one discovers 
new SIPs which are solvable by purely algebraic methods. Although the 
general problem is still unsolved, two classes of solutions have been found and 
discussed. In the first class, the parameters a\ and a 2 are related to each other 
by translation (a 2 = d\ + a) |0, Remarkably enough, all well known 
analytically solvable potentials found in most text books on nonrelativistic 
quantum mechanic s belong to this class. Last year, a second class of solutions 
was discovered in which the parameters a\ and a 2 are related by scaling 
(a a = gai) 



4.3.1 Solutions Involving Translation 

We shall now point out the key steps that go into the classification of SIPs 
in case a 2 = a\ + a ||52j| . Firstly one notices the fact that the eigenvalue 
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spectrum of the Schrodinger equation is always such that the n'th eigenvalue 
E n for large n obeys the constraint ||1 33 



l/n 2 < E n <n 2 , (126) 

where the upper bound is saturated by the square well potential and the 
lower bound is saturated by the Coulomb potential. Thus, for any SIP, the 
structure of E n for large n is expected to be of the form 

E n ~ Yl C « na ' - 2 < a < 2. (127) 

a 

Now, since for any SIP, E n is given by eq. (|101[) , it follows that if 

i?K)~E fc7 ( 128 ) 

7 

then 

-3<7<1 (129) 

How does one implement this constraint on i?(afc)? While one has no 
rigorous answer to this question, it is easily seen that a fairly general factor- 
izable form of W(x; ai) which produces the above fc-dependence in R(ak) is 
given by 

m 

W{x; ai) = J2( k i + c i)9i{x) + hi{x)/{ki + a) + fi{x) (130) 
i=i 

where 

a x = (hi, k 2 ...), a 2 = (h + a,k 2 + (3...) (131) 

with d,a,/3 being constants. Note that this ansatz excludes all potentials 
leading to E n which contain fractional powers of n. On using the above ansatz 
for W in the shape invariance condition (P^| ) one can obtain the conditions 
to be satisfied by the functions gi(x) , hi(x) , fi(x) . One important condition 
is of course that only those superpotentials W are admissible which give a 
square integrable ground state wave function. It turns out that there are no 
solutions in case m > 3 in eq. ( |130| ), while there are only two solutions in 
case m = 2 i.e. when 

W(x; a x ) = (h + ci) gi (x) + [k 2 + c 2 )g 2 (x) + fi(x), (132) 
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which are given by 

W(r; A,B) = Atanhar - Bcoth ar, A > B > 0, (133) 

and 

W(x; A, B) = A tan ax -B cot ax; A,B>0, (134) 

where < x < ir/2a and ip(x = 0) = ip(x = ir/2a) = 0. For the simplest 
possibility of m = 1, one has a number of solutions to the shape invariance 
condition ((5(J). In Table 4.1, we give expressions for the various shape invari- 
ant potentials Vi(x), superpotentials W(x), parameters a\ and a 2 and the 
corresponding energy eig envalues E$ @, |5|. 

As an illustration, let us consider the superpotential given in eq. ( |134j ). 
The corresponding partner potentials are 

Vi(x; A, B) = -(A + B) 2 + A(A - a) sec 2 ax + B(B - a)cosec 2 ax 

V 2 (x; A, B) = -(A + Bf + A(A + a) sec 2 ax + B(B + a)cosec 2 ax (135) 

V\ and V2 are often called Poschl- Teller I potentials in the literature. They 
are shape invariant partner potentials since 

V 2 (x; A,B) = V 1 (x;A + a,B + a) + {A + B + 2a) 2 - (A + B) 2 (136) 

and in this case 

{ai} = (A,B);{a 2 } = (A+a, B+a), i2(a x ) = (A+B+2a) 2 -(A+B) 2 . (137) 

In view of eq. (|101| ), the bound state energy eigenvalues of the potential 
Vi(x] A, B) are then given by 

n 

= R{a k ) = (A + B + 2na) 2 - (A + B) 2 . (138) 

k=l 

The ground state wave function of Vi(x; A, B) is calculated from the super- 
potential W as given by eq. ( |134| ). We find 

^\x-A,B) oc (cos ax) s (sin ax) (139) 

where 

s = A/a; X = B/a. (140) 
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The requirement of A, B > that we have assumed in eq. ( |134[ ) guaran- 
tees that ifjQ~\x; A, B) is well behaved and hence acceptable as x — ► 0, tt /2a. 
Using this expression for the ground state wave function and eq. ( |103j ) 
one can also obtain explicit expressions for the bound state eigenfunctions 
ip^ (x; A, B) . In particular, in this case, eq. (|103|) takes the form 

ip n (x; {ai}) = (—4 — h At&nax — B cot ax)ip n _i{x; {a 2 }). (141) 
On defining a new variable 

?/ = 1 - 2 sin 2 ax (142) 
and factoring out the ground state state wave function 

$n(v\ {ai}) = ip (y; {ai})R n (y; {«i}) (143) 

with ip being given by eq. (|139|) , we obtain 

R n (y;A.B) = a(l - y 2 )^-R n ^(y, A + a, B + a) 

dy 

+[{A-B)-{a + B + a)y)R n ^{y; A + a,B + a) (144) 



It is then clear [^6|] that R n (y; A, B) is proportional to the Jacobi poly- 
nomial so that the unnormalized bound state energy eigenfuncti ons for this 
potential are 

MV, A B) = {1- yf'\l + y y/ 2 Pt 1/2 ' 8 ~ 1/2 (y). (145) 

The procedure outlined above has been applied to all known SIPs j|6], |125| 
and the energy eigenfunctions ip^'(y) have been obtained in Table 4.1, where 
we also give the variable y for each case. 
Several remarks are in order at this time. 

1. The Poschl- Teller I and II superpotentials as given by eqs. ( |134j ) and 
(| 1 3 3|) respectively have not been included in Table 4.1 since they are 
equivalent to the Scarf I (trigonometric) and generalized Poschl- Teller 
superpotentials 

W\ = —Ataxiax + Bsecax, 

Wi = A coth ar — .Bcosech ar, (146) 
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by appropriate redefinition of the parameters p5 |. For example, one 
can write 



.A + B. w ar. ,-B — A. , ,ar. 
^2 = (—3—) tanh(— ) - (—3—) coth(— ), (147) 



which is just the Poschl- Teller II superpotential of eq. ( |133| ) with re- 
defined parameters. 

Throughout this section we have used the convention ofh = 2m — 1. It 
would naively appear that if we had not put h — 1, then the shape in- 
variant potentials as given in Table 4.1 would all be ^.-dependent. How- 
ever, it is worth noting that in each and every case, the ^-dependence 
is only in the constant multiplying the x-dependent function so that in 
each case we can always redefine the constant multiplying the function 
and obtain an ^-independent potential. For example, corresponding 
to the superpotential given by eq. (|134|) , the ^-dependent potential is 
given by (2m = 1) 



V x (x\ A,B) = W 2 -hW = -(A + B) 2 + A(A + ha) sec 2 ax 

+ B(B + ha)cosec 2 ax. (148) 

On redefining 

A(A + ha) = a; B(B + ha) = b, (149) 

where a, b are ^-independent parameters, we then have a ^-independent 
potential. 

3. In Table 4.1, we have given conditions (like A > 0, B > 0) for the su- 
perpotential (P3D , so that ^o 1} = iVexp(- J x W(y)dy) is an acceptable 
ground state energy eigenfunction. Instead one can also write down 
conditions for ip^ = N exp(J x W(y)dy) to be an acceptable ground 
state energy eigenfunction. 

4. It may be noted that the Coulomb as well as the harmonic oscillator 
potentials in n-dimensions are also shape invariant potentials. 
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5. Does this classification exhaust all shape invariant potentials? It was 



believed that the answer to the question is yes [§5], |134|| but as we shall 
see in the next subsection, the answer to the question is in fact negative. 
However, it appears that this classification has perhaps exhausted all 
SIPs where a 2 and ai are related by translation. 

No new solutions (apart from those in Table 4.1) have been obtained 
so far in the case of multi-step shape invariance and when a 2 and ai 
are related by translation. 

What we have shown here is that shape invariance is a sufficient con- 
dition for exact solvability. But is it also a necessary condition? This 



question has been discussed in detail in ref. [52] where it has been 
shown that the solvable Natanzon potentials |S7| are in general not 
shape invariant. 

Before ending this subsection, we would like to remark that for the SIPs 
(with translation) given in Table 4.1, the reflection and transmission ampli- 
tudes R\(k) and Ti(k) (or phase shift 5i(k) for the three-dimensional case) 
can also be calculated by operator methods. Let us first notice that since for 
all the cases a 2 = a\ + a, hence Ri(k; a±) and Ti(k; a\) are determined for all 
values of a-^ from eqs. ( |104[ ) and ( |105| ) provided they are known in a finite 



strip. For example, let us consider the shape invariant superpotential 

W = ntanhx, (150) 
where n is positive integer (1,2,3,...). The two partner potentials 
Vi(x; n) = n 2 — n{n + l)sech 2 a;, 

V 2 (x; n) = n 2 — n{n — l)sech 2 x, (151) 

are clearly shape invariant with 

Oi = n , a 2 = n — 1 . (152) 

On going from V\ to V 2 to V3 etc., we will finally reach the free particle 
potential which is reflectionless and for which T — 1. Thus we immediately 
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conclude that the series of potentials Vi, V2, ... are all reflectionless and the 
transmission coefficient of the reflectionless potential V\(x;n) is given by 



Ti{k,n) 



(n — ik)(n — 1 — ik)...(l — ik) 
(— n — ik)(—n + 1 — ik)...(—l — ik) 
r(— n — ik)T(n + 1 — ik) 
T(-ik)T(l-ik). 



(153) 



The scattering amplitudes of the Coulomb [I7| and the potential correspond- 
ing to W = A tanh x + i?sech x have also been obtained in this way. 



There is, however, a straightforward method |48[ for calculating the scat- 
tering amplitudes by making use of the n'th state wave functions as given in 
Table 4.1. In order to impose boundary conditions appropriate to the scat- 
tering problem, two modifications of the bound state wave functions have to 
be made: (i) instead of the parameter n labelling the number of nodes, one 
must use the wave number k! so that the asymptotic behaviour is exp(ik'x) 
as x — > 00. (ii) the second solution of the Schrodinger equation must be kept 
(it had been discarded for bound state problems since it diverged asymptot- 
ically). In this way the scattering amplitude of all the SIPs of Table 4.1 have 
been calculated in ref. 48 



4.3.2 Solutions Involving Scaling 

For almost nine years, it was believed that the only shape invariant potentials 
are those given in Table 4.1 and that there were no more shape invariant 
potentials. However, very recently we have been able to discover a huge 
class of new shape invariant potentials 59 . It turns out that for many of 



these new shape invariant potentials, the parameters 02 and a\ are related by 
scaling (a 2 = qa\,0 < q < 1) rather than by translation, a choice motivated 
by the recent interest in g-deformed Lie algebras. We shall see that many 
of these potentials are reflectionless and have an infinite number of bound 
states. So far none of these potentials have been obtained in a closed form 
but are obtained only in a series form. 

Let us consider an expansion of the superpotential of the from 

00 

W(x;a 1 ) = Y / 9j(xK (154) 

j=0 
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and further let 

<i2 = qa>i, < q < 1 . (155) 

This is slightly misleading in that a reparameterization of the form a2 = qai, 
can be recast as a' 2 = a[ + a merely by taking logarithms. However, since the 
choice of parameter is usually an integral part of constructing a SIP, it is in 
practice part of the ansatz. For example, we will construct below potentials 
by expanding in ai, a procedure whose legitimacy and outcome are clearly 
dependent on our choice of parameter and hence reparameterization. We 
shall see that, even though the construction is non-invariant, the resulting 
potentials will still be invariant under redefinition of a\. On using eqs. (|154|) 
and ( |155| ) in the shape invariance condition (|9~6"D, writting R{a{) in the form 

oo 

R( ai ) = Y, RjO-i , ( 156 ) 

3=0 



and equating powers of a± yields [£>9], [58 



2g' (x) = R ; g[{x) + 2d 1 g (x)gi(x) = ndi, (157) 

n— 1 

g' n (x) + 2d n g {x)g n (x) = r n d n - d n ^ gj(x)g n -j(x), (158) 

where 

r n = R n /(l-q n ), d n = (l-q n )/(l + q n ), n = 1,2,3,... (159) 

This set of linear differential equations is easily solvable in succession to give 
a general solution of eq. (|9~6|). Let us first consider the special case ^o(^) = 0, 
which corresponds to Rq = 0. The general solution of eq. ( |158j ) then turns 
out to be 

/n— 1 
dx[r n - Y 9j( x )9n-j(x)], n=l,2, ... (160) 
i=i 

where without loss of generality we have assumed the constants of integration 
to be zero. We thus see that once a set of r n are chosen, then the shape 
invariance condition essentially fixes the g n (x) (and hence W(x;ax)) and 
determines the shape invariant potential. Implicit constraints on this choice 
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are that the resulting ground state wave function be normalizable and the 
spectrum be sensibly ordered which is ensured if R(q n ai) > 0. 

The simplest case is r x > and r n = 0, n > 2. In this case the eq. ( |16U| ) 
takes a particularly simple form 

g n (x) = (5 n x 2n ~ l , (161) 



n-1 



(2^1)5 *!>-> (162) 



where 

f3i = tZ x n, (3 n -- 
and hence 

oo 

W(x; a x ) = J2 PAx 2 ^ 1 = y/aTF(^x) . (163) 
For a 2 = qai, this gives 

W(x;a2) = y/qW(Jqx,a 1 ) , (164) 



which corresponds to the self-similar W of Shabat and Spiridonov [ |60| , |6"Tf . It 
is worth pointing out that these self-similar potentials can be shown to satisfy 
g-supersymmetry p2]| . It may be noted here that instead of choosing r n = 
0, n > 2, if any one r n (say rj) is taken to be nonzero then one again obtains 
self-similar potentials and in these instances the results obtained from 

shape invariance and self-similarity are entirely equivalent and the Shabat- 
Spiridonov self-similarity condition turns out to be a special case of the shape 
invariance condition. 

It must be emphasized here that shape invariance is a much more general 
concept than self-similarity. For example, if we choose more than one r n to 
be nonzero, then SIP are obtained which are not contained within the self- 
similar ansatz. Consider for example, r n = 0, n > 3. Using eq. ( |160| ) one can 
readily calculate all the g n (x), of which the first three are 

g 1 (x) = dirix, g 2 (x) = d 2 r 2 x - ^-d\r\x' i , 

2 2 

9z(x) = ~7;dirid 2 r 2 d^x i + — dlrld 2 d 3 x . (165) 

3 15 

Notice that in this case W(x) contains only odd powers of x. This makes 
the potentials V\ 2 (x) symmetric in x and also guarantees unbroken SUSY. 
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The energy eigenvalues follow immediately from eqs. ( |101[ ) and ( |156| ) and 
are given by (0 < q < 1) 

pd)/ x r (l + 9)(l-9 n ) , r (l + g 2 )(l-g 2 ") , . 

(«i) = r i t; \ + r 2 7; ^ ,n = 0,l,2,.. (166) 

where I\ = G?iriai,r 2 = d 2 r 2 a\ while the unnormalized ground state wave 
function is 

2 4 

a x ) = exp[-|-(r 1 + T 2 ) + ^{d 2 T\ + 24^2 + d A T\) + 0(a; 6 )] (167) 

The wave functions for the excited states can be recursively calculated from 
the relation (|103|) . 

We can also calculate the transmission coefficient of this symmetric po- 
tential (k = k') by using the relation ( |125|) and the fact that for this SIP 
a 2 = qa\. Repeated application of the relation (|125|) gives 

[ik - W{oo, a^jik - Wjoo, a 2 )]...[ik - W{oo, a n )] 

1{ ' ai) [ik + W(oo, a x )][ik + W{oo, a 2 )]...[ik + W(oo, a n )} U ' Un+l) 

(168) 

where 

W(oo, Oj) = - Ef ] . (169) 

Now, as n — > oo,a n+1 = q n a\ — > 0(0 < q < 1) and, since we have taken 
go(x) = 0, one gets W(x; a n+ i) — > 0. This corresponds to a free particle 
for which the reflection coefficient Ri(k;a\) vanishes and the transmission 
coefficient is given by 

j=i [ik + Wfaaj)] 

The above discussion keeping only r 1; r 2 7^ can be readily generalized 
to an arbitrary number of nonzero rj. The energy eigenvalues for this case 
are given by (Tj = djTja{) 

E n\ a i) = 2^ V j (1 _ j) » n = 0,l,2,... (171) 
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All these potentials are also symmetric and reflectionless with T\ as given 
by eq. ( |170| ). The limits q — > and q — > 1 of all these potentials are simple 
and quite interesting. At q = 1, the solution of the shape invariance condition 



(96) is the standard one dimensional harmonic oscillator with W(x) = Rx/2 
while in the limit q — > the solution is the Rosen-Morse superpotential 
corresponding to the one soliton solution given by 

W{x) = v^Rtanh(v^Rx). (172) 

Hence the general solution as obtained above with < q < 1 can be regarded 
as the multi-parameter deformation of the hyperbolic tangent function with q 
acting as the deformation parameter. It is also worth noting that the number 
of bound states increase discontinuously from just one at q = to infinity for 
q > 0. Further, whereas for q = 1 the spectrum is purely discrete, for q even 
slightly less than one, we have the discrete as well as the continuous spectra. 

Finally, let us consider the solution to the shape invariance condition (|96|) 
in the case when R ^ 0. From eq. fll57j) it then follows that go(x) = R^x/2 
rather than being zero. One can again solve the set of linear differential 
equations (|157| ) and ( |158| ) in succession yielding g±(x), g^ (x),.. .Further, the 



spectrum can be immediately obtained by using eqs. ( |101| ) and ( |156| ). For 
example, in the case of an arbitrary number of nonzero Rj (in addition to 
Rq), it is given by 

ffl = ^ + gr / 1+ ff^^ (173) 



which is the spectrum of a g-deformed harmonic oscillator ||136| , |137|| . It is 
worth pointing out here that, unlike the usual q-oscillator where the space 
is noncommutative but the potential is normal (u 2 x 2 ), in our approach the 
space is commutative, but the potential is deformed, giving rise to a multi- 
parameter deformed oscillator spectrum. 

An unfortunate feature of the new SIPs obtained above is that they are 
not explicitly known in terms of elementary functions but only as a Taylor 
series about x = 0. Questions about series convergence naturally arise. Nu- 
merical solutions pose no serious problems. As a consistency check, Barclay 
et al. have checked numerically that the Schrodinger equation solved with 
numerically obtained potentials indeed has the analytical energy eigenvalues 
given above. From numerical calculations one finds that the superpotential 
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and the potential are as shown in Figs. 4.1 and 4.2 |59j corresponding to the 
case when T\ ^ 0,r n = 0,n > 2. These authors, however, see no evidence 
of the oscillations in W and V as reported by Degasperis and Shabat JB5| ]. 
A very unusual new shape invariant potential has also been obtained |p9| 
corresponding to r 1 = l,r 2 = — l,r n = 0,n > 3 (with q = 0.3 and a = 0.75) 
which is shown in Fig. 4.3. In this case, whereas Vi(x) is a double well 
potential, its shape invariant partner potential ^(x) is a single well. 

It is worth pointing out that even though the potentials are not known in 
a closed form in terms of elementary functions, the fact that these are reflec- 
tionless symmetric potentials can be used to constrain them quite strongly. 
This is because, if we regard them as a solution of the KdV equation at 
time t = 0, then being reflectionless, it is well known that such solutions 
as t — > ±oo will break up into an infinite number of solitons of the form 
2/c 4 2 sech 2 /cjX ||138|, |139||. On using the fact that the KdV solitons obey an 



infinite number of conservation laws corresponding to mass, momentum, en- 
ergy one can immediately obtain constraints on the reflectionless SIPs 
obtained above 



4.3.3 Solutions of Multi-Step Shape Invariance 

Having obtained potentials which are multi-parameter deformations of the 
Rosen-Morse potential corresponding to the one soliton solution, an obvious 
question to ask is if one can also obtain deformations of the multi-soliton 
solutions. The answer is yes [^] and as an illustration we now explicitly 
obtain the multi-parameter deformations of the two soliton case by using the 
formalism of two-step shape invariance as developed earlier in this section. 

Let us take the scaling ansatz = qa\ and expand the two superpoten- 
tials W and W in powers of a\ 

oo oo 

W(x] (Xi) = ^2 gj(x)a{; W(x; a x ) = hj(x)a{ (174) 

j=0 j=0 

Further, we write R and R in the form 

oo oo 

R( ai ) = R A\ ^K) = RA ( 175 ) 

3=0 3=0 
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Using these in eqs. ( |107| ) and ( |109| ) and equating powers of a\ yields 
(n=0,l,2,...) 

n n 

#n + E 9j9n-j = E h i h n-j ~ K + R n (176) 
j=0 j=0 

n oo 

K + E h 3 h n-3 = q n E 9i9 n -3 - v n g'n + Rn (177) 

3=0 j=0 

This set of linear equations is easily solved in succession. For example when 
R = R = (and hence go(x) = h (x) = 0) and further R n = R n = 0, n > 3 
one can readily calculate all the g n (x) and h n (x). The first two of each are 

(Ri - Ri) 

92 = ( ^ 2 rff x + 37T^ [(1 - - R ^ - 2(1 + 9)RiRi)Lm 

(Ri - qRi) R 2 x 

h\ = 7 : X- 



2\ 



[l-q) ' {l-q 

X ^2 , /1 „\ /i „2\ D 2 



■[(1 + q)R{ + (1 + g)(l - q 2 )Ri - 2q(l - q)R 1 R 1 ] 



3(l + g)(l-g) 2 

(179) 



The energy eigenvalues which follow from eqs. Q117 ) and ( |118| ) are |5 



4n\ai) = EiRi+RMj^TJ 



2 



1 - q jn ) 



qi) 



7=1 U — 9 J j=l I 1 — 9 J 



For the special case of R 2 = Rz = the spectrum was obtained previously 
by Spiridonov from the considerations of the two-step self-similar potentials 

The limit q — > of the above equation is particularly simple and yields 
the Rosen-Morse potential corresponding to the two soliton case i.e. 

W = 2\[& tanh \[&x, W = V# tanh \[&x (181) 
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provided R = 3R. This procedure can be immediately generalized and one 

can consider shape invar iance with a scaling ansatz in 3,4, p steps and 

thereby obtain multi-parameter deformations of the 3,4, ...p soliton Rosen- 
Morse potential. 



4.3.4 Other Solutions 

So far we have obtained solutions where a 2 and a\ are related either by scaling 
or by translation. Are there shape invariant potential where a 2 and a\ are 
neither related by scaling nor by translation? It turns out that there are 
other possibilities for obtaining new shape invariant potentials. Some of the 
other possibilities are: a 2 = qa? x with p— 2,3,... and a 2 = qa x /(l + pa x ). Let 
us first consider the case when 

a 2 = qa\ (182) 



i.e p = 2. Generalization to arbitrary p is straightforward J59j. On using eqs. 
( |155| ) and ( |156|) one obtains the set of equations 



2m m 

92m( X ) +J29j( X )92m-j(x) = q™ 9j (x)9m-j (x) ~ q™ 9m( X ) + R 2m, (183) 
3=0 3=0 

2m+l 

92m+l( X )+ J2 9j( x )92m+l-j(x) = R 2m +i, (184) 
3=0 

which can be solved in succession and one can readily calculate all the g n (x). 
For example, when only R\ and R 2 are nonzero, the first three g's are 

gi(x) = Rxx, g 2 (x) = (R 2 - qRi)x - ^R^x 3 , 

g 3 {x) = ^-Ri{qRi - R 2 )x 3 + ^-R\x\ (185) 
3 15 

The corresponding spectrum turns out to be (Eq (ai) = 0) 

p n p n 

= ^ + § U^f^ = 1. 2, (186) 

q ,=i q 3=1 

The q — > limit of these equations again correspond to the Rosen-Morse 
potential corresponding to the one soliton solution. One can also consider 
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shape invariance in multi-steps along with this ansatz thereby obtaining de- 
formations of the multi-soliton Rosen-Morse potential. 

One can similarly consider solutions to the shape invariance condition 
(pop in case 

ai = -p- (187) 
1 +pai 

when < q < 1 and pa± << 1 so that one can expand (1 +pa\)~ 1 in powers 
of a\. For example, when only R\ and R2 are nonzero, then one can show 
that the first two nonzero g n are 

giix) = JTT7V ^x) = (R 2 + JTTV) )x- (1 + g)2(1 + g2) y (188) 
and the energy eigenvalue spectrum is = 0) 

Ei^ Rl ± ^ + R 2 ± M-yi (189) 

Generalization to the case when several i?j are nonzero as well as shape 
invariance in multi-steps is straight forward. 

We would like to close this subsection with several comments. 



1. Just as we have obtained q- deformations of the reflectionless Rosen- 
Morse and harmonic oscillator potentials, can one also obtain deforma- 
tions of the other SIPs given in Table 4.1? 

2. Have we exhausted the list of SIP? We now have a significantly ex- 
panded list but it is clear that the possibilities are far from exhausted. 
In fact it appears that there are an unusually large number of shape 
invariant potentials, for all of which the whole spectrum can be ob- 
tained algebraically. How does one classify all these potentials? Do 
these potentials include all solvable potentials |52|, |53[| ? 

3. For those SIP where a 2 and a\ are not related by translation, the spec- 
trum has so far only been obtained algebrically. Can one solve the 
Schrodinger equation for these potentials directly? 
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4. There is a fundamental difference between those shape invariant poten- 
tials for which ci2 and a\ are related by translation and other choices 
(like ci2 = qa,i). In particular, whereas in the former case the poten- 
tials are explicitly known in a closed form in terms of simple functions, 
in the other cases they are only known formally as a Taylor series. 
Secondly, whereas in the latter case, all the SIP obtained so far have 
infinite number of bound states and are either reflectionless (or have 
no scattering), in the former case one has also many SIP with nonzero 
reflection coefficients. 



4.4 Shape Invariance and Noncentral Solvable Poten- 
tials 

We have seen that using the ideas of SUSY and shape invariance, a number 
of potential problems can be solved algebraically. Most of these potentials 
are either one dimensional or are central potentials which are again essen- 
tially one dimensional but on the half line. It may be worthwhile to enquire 
if one can also algebraically solve some noncentral but separable potential 
problems. As has been shown recently ||131|| , the answer to the question is 



yes. It turns out that the problem is algebraically solvable so long as the 
separated problems for each of the coordinates belong to the class of SIP. 
As an illustration, let us discuss noncentral separable potentials in spherical 
polar coordinates. 

In spherical polar coordinates the most general potential for which the 
Schrodinger equation is separable is given by ||140| 



V(r, 9, 0) = V(r) + M + (190) 

r z r 1 sin 9 

where V(r), V(9) and V(4>) are arbitrary functions of their argument. First, 
let us see why the Schrodinger equation with a potential of the form given 
by eq. ( |190| ) is separable in the (r, 9, (p) coordinates. The equation for the 
wave function tp(r,9,<f)) is 



d 2 ip 2dip. 1 ,d 2 ip dip. 1 <fty 



(E-V)ip (191) 
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It is convenient to write ip(r, 6, 4>) as 



^(r,e,<p) = R{r \ H{ ° ) K^). (192) 



Substituting eq. Q192] ) in eq. ( |191| ) and using the standard separation of 
variables procedure, one obtains the following equations for the functions 
K{<t>),H{9) and R{r) : 

— +y(^(^ m 2 ^), (193) 

- + [V(0) + (m 2 - ~) cosec 2 #] ff(0) = l 2 H(8), (194) 

- ^ + + ^-^1 ^( r ) = E ^( r )' ( 195 ) 

where m 2 and / 2 are separation constants. 

As has been shown in ref. [|131f] , the three Schrodinger equations given by 
(|193|) , (|194|) and (|195|) may be solved algebraically by choosing appropriate 



SIPs for V(4>), V(8) and V(r). Details can be found in ref. ||131| 



Generalization of this technique to noncentral but separable potentials in 
other orthogonal curvilinear coordinate systems as well as in other dimensions 
is quite straightforward. Further, as we show below, one could use this trick 
to discover a number of new exactly solvable three-body potentials in one 
dimension. 



4.5 Shape Invariance and 3-Body Solvable Potentials 

Many years ago, in a classic paper, Calogero ||141|| gave the complete solution 
of the Schrodinger equation for three particles in one dimension, interacting 
via two-body harmonic and inverse-square potentials. Later, Wolfes | |1 12| 



used Calogero's method to obtain analytical solutions of the same problem 
in the presence of an added three-body potential of a special form. Attention 
then shifted to the exact solutions of the many-body problem and the general 
question of integrability ||143 , 144 , 145 ]. Recently, there has been renewed 



interest in the one- dimensional many-body problem in connection with the 
physics of spin chains [ |146| , |147| , |148| , [L49[| . Also, there has been a recent 
generalization of Calogero's potential for N-particles to SUSY QM 
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The purpose of this subsection is to show that using the results for SIPs 
derived above, one can discover a number of new 3-body potentials for which 
the 3-body problem in one dimension can be solved exactly ||132|| . There is 
a rule of thumb that if one can solve a 3-body problem then one can also 
solve the corresponding n-body problem. Thus hopefully one can also solve 
the corresponding statistical mechanics problem. 

The important point to note is that three particles in one dimension, 
after the center of mass motion is eliminated, have two degrees of freedom. 
This may therefore be mapped on to a one-body problem in two dimensions. 
Calogero [p.41|| considered the case where the two dimensional potential is 
noncentral but separable in polar coordinates r, 0. From the above discussion, 
it is clear that if the potentials in each of the coordinates r and are chosen 
to be shape invariant, then the whole problem can be solved exactly. 

Calogero's [p.41]] solution of the three-body problem is for the potential 



c 



^/8E( 



X 7 ' 



x 



,) 2 + sE< 

i<j 



X 7 ' 



(196) 



where g > — 1/2 to avoid a collapse of the system, 
a three-body potential 



Wolfes 142 showed that 



w 



/[(xi + x 2 - 2x 3 )~ 2 + (x 2 + x 3 - 2xi)" 2 + (x 3 + x 1 - 2x 2 )- 2 ] (197) 



is also solvable when it is added to Vc, with or without the pairwise cen- 
trifugal term. The last two terms on the right-hand side of eq. (|197|) are 
just cyclic permutations of the first. Henceforth, such terms occuring in 
any potential are referred to as "cyclic terms". In this subsection, we give 
more examples of three-body potentials that can be solved exactly. Our first 
solvable example is the three-body potential of the form 



2r 2 



(xi + x 2 - 2X2 
{xi - x 2 ) 



+ cyclic terms 



(198) 



which is added to the Calogero potential Vc of eq. ( |196j ). In the above 
equation, 



-\(xi - x 2 f + (x 2 - x 3 ) 2 + (x 3 - Xi) 



(199) 



To see why potentials given by eqs. (|196|) to (|198|) are solvable, define the 
Jacobi coordinates 

R= -( Xl +x 2 + x 3 ), (200) 
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x= {Xl ~ X2 \ y= (*i + *2- 2 **). (20 i) 

After elimination of the center of mass part from the Hamiltonian, only 
the x- and y-degrees of freedom remain, which may be mapped into polar 
coordinates 

x = r sin , y = r cos 4> . (202) 

Obviously, the variables r, <p have ranges < r < oo and < <fi < 2n. It is 
straightforward to show that 

(xi — x 2 ) = V2rsin(j), 

(x 2 — x 3 ) = a/2 r sin(0 + 27r/3), 

(X3 — Xi ) = V2r sin(0 + 47r/3). (203) 

It turns out that Vc, Vw as well as V\ are all noncentral but separable po- 
tentials in the polar coordinates r,0. As a result, the Schrodinger equation 
separates cleanly in the radial and angular variables, and the wave function 
can be written as 

Vvz(r,0) = ^W). (204) 



if 

In all three cases, the radial wave function obeys the equation 



d 2 3 2 2 (A 2 -±) 

1 — co 1 r H 

dr 2 8 r 2 



+ -uo 2 r 2 + V l — A> R nl = E nl R nl (r) , (205) 



where Bf is the eigenvalue of the Schrodinger equation in the angular vari- 
able. Eq. ( P05| ) corresponds to a SIP and the eigenvalues E nl and the eigen- 
functions R n i which follow from Table 4.1 are 



E n i = v3/2 u(2n + Bi + 1) , n,/ = 0,l,2... (206) 



Rni = r B > exp[^( V / 372)^r 2 ]L^[i( V / 372)^r 2 ], (207) 



where Bi > 0. To examine the angular part of the eigenfunction Fi((j>), take 
the potential Vc + V\. Then, in the variable 0, the Schrodinger equation is 



d 



' +f E cosec 2 [0 + 2(m + E cot [0 + 2(m - F,(0) 

^ m=l 6 1 m=l 15 J 
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= Bf Fi((p). (208) 

On using the identities 

3 

cosec 2 [0 + 2(m - l)vr/3] = 9cosec 2 30 (209) 



m=l 

3 



cot[0 + 2(m- 1)tt/3] = 3cot 30 (210) 



m=l 



eq. ( [2081 ) reduces to 

{"^ + ^cosec 2 30+ ^Acot 30 j ^ = Bf F^). (211) 
Now this is again a SIP and hence its eigenvalues and eigenf unctions are 

Bf = 9(l + a + 1/2) 2 - ^-fl/{l + a + 1/2) 2 , (212) 

lb 

Fi{4>) = exp(-^Z)(sin 30) [ exp[|A^]P~'~^ cot 30) (213) 

where -P" ,/3 is the Jacobi polynomial and 

a = 1/2(1 + 2#) 1/2 . (214) 

It must be noted that g > —1/2 for meaningful solutions. As expected, we 
recover the results for the Calogero potential Vc in the limit f\ = 0. 

For distinguishable particles, a given value of defines a specific order- 
ing. For < < 7r/3, eq. (|203| ) implies X\ > x 2 > x%, and other ranges 
of correspond to different orderings [fL41| , |142| ]. For singular repulsive 
potentials, crossing is not allowed, and Fi((f>) of eq. (|213| ) is zero outside 
< < 7r/3. Following Calogero, the wave function for the other ranges 
may be constructed. Similarly, for indistiguishable particles, symmetrized or 
antisymmetrized wave functions may be constructed. 

Proceeding in the same way and using the results of Table 4.1, it is easily 
show n that the other exactly solvable potentials are [ 132 

V<i = -ijj 2 y^X x i ~ x j) 2 + 3g [(xi + X2 — 2xs)~ 2 + cyclic terms ] 



_3\/3 fi 

9 



(xi-x 2 ) 

+ cyclic terms 



(xi + x 2 — 2x s 



(215) 
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V 3 = ^ 2 J2( x ~ x i) 2+ 9j2( x i- x j 



i<j 



1 



(xi + x 2 - 2x 3 ) 
(xi - x 2 ) 2 



+ cyclic terms 



(216) 



V4 = -oj 2 y^X x i ~ x j) 2 + 3 g [ (xi + x 2 — x 3 ) 2 + cyclic terms 

8 i<j 



1 h 



3y/2 r 



(xi - x 2 ) 



{x\ + x 2 - 2x 3 ) 2 



+ cyclic terms 



(217) 



Notice that in all these cases one has combined the SIP in the angular variable 
with the harmonic confinement. The three body scattering problem has also 
been studied in these cases after droping the harmonic term |132f| . The 
possibility of replacing the harmonic confinement term with the attractive — 
type interaction has also been considered. Note that in this case one has both 
discrete and continuous spectra. Further, in this case one can obtain exact 
solutions of three-body problems for all the SIPs discussed above (i.e.Vi,... 
V4, Vc, Vw, Vc + Vw) along with the attractive - potential ||132 . 



Operator Transforms — New Solvable Po- 
tentials from Old 



In 1971, Natanzon |56] wrote down (what he thought at that time to be) 
the most general solvable potentials i.e. for which the Schrodinger equation 
can be reduced to either the hypergeometric or confluent hypergeometric 
equation . It turns out that most of these potentials are not shape invari- 
ant. Further, for most of them, the energy eigenvalues and eigenf unctions 
are known implicitly rather than explicitly as in the shape invariant case. 
One might ask if one can obtain these solutions from the explicitly solv- 
able shape invariant ones. One strategy for doing this is to start with a 
Schrodinger equation which is exactly solvable (for example one having a 
SIP) and to see what happens to this equation under a point canonical coor- 
dinate transformation. In order for the Schrodinger equation to be mapped 
into another Schrodinger equation, there are severe restrictions on the nature 
of the coordinate transformation. Coordinate transformations which satisfy 
these restrictions give rise to new solvable problems. When the relationship 
between coordinates is implicit, then the new solutions are only implicitly 
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determined, while if the relationship is explicit then the newly found solvable 
potentials are also shape invariant [ 127 , 128 , 129]. In a more specific special 
application of these ideas, Kostelecky et al. | |13U|| were able to relate, using 
an explicit coordinate transformation, the Coulomb problem in d dimensions 
with the d-dimensional harmonic oscillator. Other explicit applications of the 
coordinate transformation idea are found in the review article of Haymaker 
and Rau |20 |. 

Let us see how this works. We start from the one- dimensional Schrodinger 
equation 

f +[V(x)-E n ])i; n (x) 



dx 2 



0. 



Consider the coordinate transformation from x to z defined by 

m = £, 

dx 



so that 



d d 
dx dz 



(218) 



(219) 



(220) 



The first step in obtaining a new Schrodinger equation is to change coordin 
ates and divide by f 2 so that we have: 



d 2 



dz 2 



f d [V-E n 
fdz + f 2 



4>r< 







(221) 



To eliminate the first derivative term, one next rescales the wave function: 

i> = r 1/2 i>- (222) 
Adding a term e n ip to both sides of the equation yields 

d 2 i) n 



where 



dz< 



V{E n 



+ [V(E n ) + e n ]ip n = e n ip n} 



V ~ E n 
P 



r 



r 
v 



(223) 



(224) 



In order for this to be a legitimate Schrodinger equation, the potential 
V(E n ) + e n must be independent of n. This can be achieved if the quantity 
G defined by 

G = V - E n + e n f (225) 
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is independent of n. How can one satisfy this condition? One way is to have 
P and G to have the same functional dependence on x(z) as the original 
potential V. This further requires that in order for V to b e independent 
of n, the parameters of V must change with n so that the wave function 
corresponding to the n'th energy level of the new Hamiltoni an is related to 
a wave function of the old Hamiltonian with parameters which depen d on 
n. This can be made clear by a simple example. 

Let us consider an exactly solvable problem - the three dimensional har- 
monic oscillator in a given angular momentum state with angular momentum 
(5. The reduced ground state wave function for that angular momentum is 

^ (r) = r /3+V Qr2/2 , (226) 

so that the superpotential is given by 

W(r)=ar-((3+l)/r, (227) 

and Hi is given by: 

#i = ~ + + « V - MP + 3/2). (228) 

By our previous argument we must choose / = ^ to be of the form: 

f 2 = (^) 2 = - 2 +Br 2 + C. (229) 
ar r 2 

The solution of this equation gives z = z(r) which in general is not invertible 
so that one knows r = r(z) only implicitly as discussed before. However for 
special cases one has an invertible function. Let us for simplicity now choose 

f = r, z = r 2 /2 . (230) 

As discussed earlier, the energy eigenvalues of the three dimensional harmonic 
oscillator are give by 

E n = Aan , (231) 
so that the condition we want to satify is 

V - Aan + e n f = G = ^ + Er 2 + F. (232) 
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Equating coefficients we obtain 



D = + 

E = e n + a 2 = 7 , 

F = -Aan - 2a((3 + 3/2) = -2Ze 2 . (233) 

We see that for the quantities D, E, F to be independent of n one needs 
to hav e that a, which describes the strength of the oscillator potential, be 
dependent on n. Explicitly, solving the above three equations and choosing 
(3 = 21 + 1/2, we obtain the relations: 

Ze 2 Z 2 e 4 

a ^ n) = 2(i + n + i) ' e " = 7 - WTTTW ' (234) 

We now choose 7 = a 2 (l,n = 0) so that the ground state energy is zero. 
These energy levels are those of the hydrogen atom. In fact, the new Hamil- 
tonian written in terms of z is now 

H d2 1 l l ± X ) Ze2 1 Z ^ A (235) 
H ~ dz* + z + 4(/ + l)2 [Z6b) 

and the ground state wave function of the hydrogen atom is obtained from 
the ground state wave function of the harmonic oscillator via 

Higher wave functions will have values of a which depend on n so that the 
different wave functions correspond to harmonic oscillator solutions with dif- 
ferent strengths. 

All the exactly solvable shape invariant potentials of Table 4.1 can be 
inter-related by point canonical coordinate transformations ||128| , |129|| . This 
is nicely illustrated in Fig. 5.1 . In general, r cannot be explicitly found in 
terms of z, and one has 



dz/dr = f = ^A/r 2 + Br 2 + C , (237) 

whose solution is: 
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z 



VA + Cr 2 + Br A VA log(r 2 ) 



2 2 



y/A log(2 A + Cr 2 + 2 \f~A VA + Cr 2 + Br 4 ) 



2 



+ 



C log(C + 2£r 2 + 2v A B v M + Cr 2 + £r 4 ) 
1 v'77 



(238) 



This clearly is not invertible in general. If we choose this general coordinate 
transformation, then the potential that one obtains is the particular class 
of Natanzon potentials whose wave functions are confluent hypergeometric 
functions in the variable r and are thus only implicitly known in terms of the 
true coordinate z. In fact, even the expression for the transformed potential 
is only known in terms of r: 



and thus only implicitly in terms of z. Equating coefficients, we get an 
implicit expression for the eigenvalues: 



as well as the state dependence on a and (5 necessary for the new Hamiltonian 
to be energy independent: 



5.1 Natanzon Potentials 

The more general class of Natanzon potentials whose wave functions are 
hypergeometric functions can be obtained by making an operator transfor- 
mation of the generalized Poschl- Teller potential whose Hamiltonian is: 



V(z, D, E) = 1/flD/r 2 + Er 2 + F- f 2 /4 + ff/2] 



(239) 




(240) 




(241) 




(242) 



This corresponds to a superpotential 



W = atanhr — /3coth r 



(243) 
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and a ground state wave function given by: 

•0o = sin]/ r cosh - " r . (244) 

The energy eigenvalues were discussed earlier and are 

E n = (a - (3) 2 - (a - (3 - 2n) 2 . (245) 

The most general transformation of coordinates from r to z which pre- 
serves the Schrodinger equation is described by: 

f = -^-2 tV + C=(^) 2 (246) 

smlTr coslrr dr' 

From this we obtain an explicit expression for z in terms of r. 

r- , t -3A + B-C + (A + B-C) cosh 2r 

z = VA tan x { } 

2 v 7 ^ J -2 A + 2B-C + 2(A + B) cosh 2r + C cosh 2 2r 
- v^B \og{-A + 3B-C + (A + B + C) cosh 2r 



+ 2v^ \j-2 A + 2B-C + 2(A + B) cosh 2r + C cosh 2 2r} 

+ VC hg{A + B + C cosh 2r 

+ VC yf-2 A + 2B-C + 2(A + B) cosh 2r + C cosh 2 2r} 

+ v^B log{2 sinh 2 r} (247) 

However the expression for the transformed potential is only known in terms 

V(=r,;5)- 1 7(7 + 1) i Q f2 i -) (248) 

f 2 sinh 2 r cosh 2 r 4/ 2 2/ 

and thus only implicitly in terms of z. Equating coefficients, we get an 
implicit expression for the eigenvalues: 

[( 7 + 1/2) 2 - Ae n Y' 2 - [(5 - 1/2) 2 - Btn] 1 ' 2 -{a- Ce n ) 1/2 = 2n+l (249) 

as well as the state dependence on a and (5 necessary for the new Hamiltonian 
to be energy independent: 

a n = [( 7 + l/2) 2 -Ae n ] 1 / 2 -l/2, 

(3 n = [{5-l/2) 2 -Be n ] l l 2 + l/2. (250) 



61 



5.2 Generalizations of Ginocchio and Natanzon Poten- 
tials 

In ref. ]52] it was explicitly shown that the Ginocchio and Natanzon po- 
tentials, whose wave functions are hypergeometric functions of an implicitly 
determined variable are not shape invariant and thus one could construct 
new exactly solvable potentials from these using the factorization method. 
It is convenient to use the original approach []S7|, ^B] to define the Ginoc- 
chio and Natanzon potentials and then obtain their generalization to more 
complicated solvable potentials whose wave functions are sums of hyperge- 
ometric functions. These generalized potentials have coordinates which are 
only implicitly known. The operator which allows one to construct the eigen- 
functions of H 2 from those of Hi converts single hypergeometric functions of 
the implicitly known coordinate to sums of hypergeometric functions. This 
process yields totally new solvable potentials. The resulting potentials are 
ratios of polynomials in the transformed coordinate in which the wave func- 
tions are sums of hypergeometric functions. The transformed coordinate is 
only implicitly known in terms of the coordinate system for the Schrodinger 
equation. 

The original method of obtaining the Natanzon potentials was to find a 
coordinate transformation that mapped the Schrodinger equation onto the 
equation for the hypergeometric functions. 

If we denote by r the coordinate appearing in the Schrodinger equation : 

l-^ + v ( r )M r ) = Q ( 251 ) 

where — oo < r < oo and z the coordinate describing the hypergeometric 
function F(a, j3; 7; z) appearing in the differential equation: 

d 2 F dF 
z (l- z )— + [ 1 -(a + p+ l)z\ — - a(3F = 0, (252) 



where < z < 1, then the transformation of coordinates is given by 
dz 2z(l-z) 



dr ' R 1 / 2 



R = az + (ci — Co — a)z + cq (253) 



In terms of the coordinate z the most general potential V(r) which cor- 
respon ds to a Schrodinger equation getting mapped into a hypergeometric 
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function equation is given by: 



V(r) 



[fz(z -l) + h (l-z) + h lZ + l]/R + 
a + (ci -c )(2z- 1) 5A z 2 (l-zf 

z(z-i) Ir* m 



{a + 



(254) 



where 



A = (a - c - ci) 2 - 4c ci. 



We note that z is only implicitly known in terms of r. The potential V(r) 
is a function of six dimensionless parameters /, ho, hi, a, Cq and c±. These six 
parameters will be related to the energy e n and the parameters (a, (3, 7) of 
the hypergeometric function below. 

The class of potentials called the Ginocchio potentials is a subclass of the 
Natanzon class which has only two independent variables nu and A where 



Co = 0, ci = 1/A 4 , a = ci-l/A 2 , h = -3/4, h x = -1, and / = (z/ + l/2) 2 -l. 



The Ginocchio potentials are more easily discussed, however by con- 
sidering them as derived from coordinate transformations which map the 
Schrodinger equation into the differential equation for the Gegenbauer poly- 
nomials. 

The transformation is given by 



V(r) = {-A 2 ^ + l) + i (l-A 2 )[5(l-AV-(7-AV + 2]}(l-y 2 ). (256) 



Although y is not known explictly in terms of r, it is implicitly known via: 



rA 2 = tanlr 1 ^) - (1 - A 2 ) 1/2 tanh- 1 [(l - X 2 ) 1/2 y], A < 1, 
rA 2 = tanh- 1 (y) + (A 2 -l) 1/2 tan- 1 [(A 2 -l) 1/2 y], A > 1. (257) 



and changing the variable x to an angle x = cos 6, one obtains from the 
Schrodinger equation, the equation for the Gegenbauer polynomials: 




(255) 



1 
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+ („ + n + 1/2) 2 - ^-4^]Z^ n = (258) 
du z sin 

where 

v _ \i-y 2 \ l/2 

Z ~ A 

e n = -/i 2 A 4 . (259) 
The corresponding energy eigenvalues and eigenfuctions are 

yU„A 2 = [A> + 1/2) 2 + (1 - A 2 )(n + 1/2) 2 ] 1 / 2 - (n+ 1/2), (260) 

W» = (1 - i/ 2 )^ /2 [^(i/)]- (2 ^ +1)/4 ^ +1 / 2 )(x). (261) 
Thus the superpotential is : 

W(r) = = i(l - \ 2 )y(y 2 - 1) + /m/A 2 (262) 

For the Natanzon potential ( |254] ) we have that the energy eigenvalues are 
given by: 

(2n + l) = (1 - ae n + /) 1/2 - (1 - c e n + h ) 1/2 

(l-C^ + htfP^On-Pn-Sn. (263) 

and the corresponding (unnormalized) energy eigenfunctions are given by : 

i/j n = i? 1/4 ^" /2 (l - z) s "' 2 F(-n, a n -n;l + (3 n ; z). (264) 

It is easy to determine the superpotential from the ground state wave 
function: 

W(r) = [*o*-(l + &)(l-*)]/i2 1/2 

+ [(c 1 -c -a)z + 2c ](l-z)/(2R 3/2 ). (265) 

Once we have the superpotential and the explicit expression for the wave 
functions we can determine all the wave functions of the partner potential 
using the operator: 
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a d \ dz d TJT , . 2z(l - z) d Tjr . . 

A ^Tr + W(r) = TrTz + W(Z) = + W(Z) 

for the case of the Natanzon potential and the operator: 

Al = Tr + W{r) = tij + W{V) = {l ~ y2)[1 + (A ' ~ l)y2] Ty + W{V) 

for the case of the Ginocchio potentials. In general, these operators take a 
single hypergeometric function into a sum of two hyper geometric functions 
which cannot be reexpressed as a single hypergeometric function except in 
degenerate cases where one obtains the shape invariant special cases discussed 
earlier. Explicitly, one has for the unnormalized eigenfunctions of if 2, the 
partner to the origina 1 Natanzon Hamiltonian: 

A = i?- 1/ V"/ 2 (l- 2 ) 5 "/ 2 [(l-z)[/3„-/5o-^(5„-5 )]i ? (-^,a„-n;l+/5„;z) 

- 22(1 - z)n{n + 1 + p n + S n )F(-n + 1, a n - n + 1; 2 + f3 n ; z)} (266) 
and for the partner potential V2 = W 2 + W one has 

V 2 = e + p + S )(J3o + 5 + 2)z(z - 1) + ((3 2 - 1)(1 - z) + (5 2 - l)z + 1]/R 

+[a - [ Ci (3Aj + 5 ) + c ((3 + 35 )a(l -fa- 5 )]/[z(z - 1)] 

-(2z - l)[( Cl - c )(/3o + S + 1) + a(P - 5 )]/[z(z - 1)] 

+ 7A/(4R)}z 2 (l - z 2 )/R 2 . (267) 

By using the hierarchy of Hamiltonians one can construct in the usual manner 
more and more complicated potentials which are ratios of higher and higher 
order polynomials in z (as well as R) which are isospectral to the original 
Natanzon potential except for the usual missing states. The wave functions 
will be sums of hypergeometric functions. The arguments are identical for the 
Ginocchio class with the hypergeometric functions now being restricted to 
being Gegenbauer polynomials. The potentials one obtains can have multiple 



local minima and several of these are displayed in [52]. In [52 1 it was also 



shown that the ser ies of Hamiltonians and SUSY charges form the graded 
Lie algebra si (1/1) ® SU{2). However this algebra did not lead to any new 
insights. We shall also see in Sec. 13 that the series of Hamiltonians form a 
parasupersymmetry of order p if the original Hamiltonian has p bound states. 
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6 Supersymmetric WKB Approximation 



The semiclassical WKB method [151] is one of the most useful approxima- 



tions for computing the energy eigenvalues of the Schrodinger equation. It 
has a wider range of applicability than standard perturbation theory which is 
restricted to perturbing potentials with small coupling constants. The pur- 
pose of this section is to describe and give applications of the supersymmetric 
WKB method (henceforth called SWKB) |79| , |80H which has been inspired by 
supersymmetric quantum mechanics. 



6.1 SWKB Quantization Condition for Unbroken Su- 
per symmetry 

As we have seen in previous sections, for quantum mechanical problems, the 
main implication of SUSY is that it relates the energy eigenvalues, eigenfunc- 
tions and phase shifts of two supersymmetric partner potentials V\(x) and 
V^x). Combining the ideas of SUSY with the lowest order WKB method, 
Comtet, Bandrauk and Campbell |79|] obtained the lowest order SWKB quan- 
tization condition in case SUSY is unbroken and showed that it yields energy 
eigenvalues which are not only exact for large quantum numbers n (as any 
WKB approximation scheme should in the semiclassical limit) but which are 
also exact for the ground state (n = 0). We shall now show this in detail. 

In lowest order, the WKB quantization condition for the one dimensional 
potential V(x) is ||151f] 

[ X2 j2m[E n - V(x)] dx = (n + l/2)7wr, n = 0,1,2,... , (268) 

where x\ and x 2 are the classical turning points defined by E n = V(x\) = 
V(x 2 ). For the potential V\(x) corresponding to the superpotential W(x), 
the quantization condition (|268| ) takes the form 

[ X2 J2m[Ei 1] - W 2 {x) + -JL=W'(x)] dx = (n + l/2)fm. (269) 
J xi V y2m 

Let us assume that the superpotential W(x) is formally 0{%). Then, the 
W term is clearly 0(h). Therefore, expanding the left hand side in powers 
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of h gives 

t ^2m[E^ - W 2 (x)] dx + % f" -, W '} X) dX + ... = (n + l/2)h7r, (270) 

J a 2Ja JE { n ) -W 2 (x) 



where a and b are the turning points defined by E^ = W 2 (a) = W 2 (b). The 
0(h) term in eq. (|270|) can be integrated easily to yield 

W(x) 



h 

— sin 

2 



(271) 



In the case of unbroken SUSY, the superpotential W(x) has opposite signs 
at the two turning points, that is 



- W(a) = W(b) = . (272) 

For this case, the 0(h) term in Q271|) exactly gives hTr/2, so that to leading 
order in h the SWKB quantization condition when SUSY is unbroken is 



fraiHof 

/ ^2m[E { n ] - W 2 (x)] dx = nhrr, n = 0,1,2,... (273) 

J a 

Proceeding in the same way, the SWKB quantization condition for the po- 
tential V 2 (x) turns out to be 

[ b \/2m[Ei 2) - W 2 (x)} dx=(n + l)hn, n = 0,1,2,... (274) 

J a 

Some remarks are in order at this stage. 

(i) For n = 0, the turning points a and b in eq. (|273|) are coincident and 



Eq = 0. Hence SWKB is exact by construction for the ground state energy 
of the Hamiltonian H 1 = (—h 2 /2m)d 2 /dx 2 + V\(x). 



(ii) On comparing eqs. ( |273| ) and ( |274j ), it follows that the lowest order 
SWKB quantization condition preserves the SUSY level degeneracy i.e. the 
approximate energy eigenvalues computed from the SWKB quantization con- 
ditions for V\(x) and V 2 (x) satisfy the exact degeneracy relation -E^+x = E^\ 

(iii) Since the lowest order SWKB approximation is not only exact as 
expected for large n, but is also exact by construction for n = 0, hence, 
unlike the ordinary WKB approach, the SWKB eigenvalues are constrained 
to be accurate at both ends of allowed values of n at least when the spectrum 
is purely discrete. One can thus reasonably expect better results than the 
WKB scheme. 
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6.2 Exactness of the SWKB Condition for Shape In- 
variant Potentials 

How good is the SWKB quantization condition [eq. ( |273| )]? To study this 
question, the obvious first attempts consisted of obtaining the SWKB bound 
state spectra of several analytically solvable potentials like Coulomb, har- 
monic oscillator, Morse, etc. |79], |80| . In fact, it was soon shown that the 
lowest order SWKB condition gives the exact eigenvalues for all SIPs [fl5j|! 
The proof of this statement follows from the facts that the SWKB con- 
dition preserves the level degeneracy and a vanishing ground state energy 
eigenvalue. For the hierarchy of Hamiltonians discussed in Sec. 3, the 
SWKB quantization condition takes the form 



\ 



2m[E { n ) - R( a k) ~ W 2 (a s ; x)] dx = ntm . (275) 

k=l 



Now, since the SWKB quantization condition is exact for the ground state 
energy when SUSY is unbroken, hence 

E { S) = ± R(a k ) (276) 

k=l 

as given by eq. (|275|) , must be exact for Hamiltonian H^ s \ One can now go 
back in sequential manner from to to and and use the 

fact that the SWKB method preserves the level degeneracy -E^+x = . On 
using this relation n times , we find that for all SIPs, the lowest order SWKB 
condition gives the exact energy eigenvalues |fl5|| . This is a very substantial 
improvement over the usual WKB formula eq. ( |268| ) which is not exact 
for most SIPs. Of course, one can artificially restore exactness by ad hoc 



Langer-like corrections [152]. However, such modifications are unmotivated 
and have different forms for different potentials. Besides, even with such 
corrections, the higher order WKB contributions are non-zero for most of 



these potentials| |152| , |153| . 



What about higher order SWKB contributions? Since the lowest order 
SWKB energies are exact for shape invariant potentials, it would be nice to 
check that higher order corrections vanish order by order in h. By starting 
from the higher order WKB formalism, one can readily develop the higher 
order SWKB formalism p3l. It has been explicitly checked for all known 
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SIPs that up to 0(h 6 ) there are indeed no corrections.. This result can be 
extended to all orders in % [jS4 , BjJ. Conditions on the superpotential which 
ensure exactness of the lowest order SWKB condition have been given in ref. 



so 



It has been proved above that the lowest order SWKB approximation 
reproduces the exact bound state spectrum of any SIP. This statement has 
indeed been explicitly checked for all SIPs known until last year i.e. solutions 
of the shape invariance condition involving a translation of parameters 02 = 
Oi+ constant. However, it has recently been shown [|154| that the above 



statement is not true for the newly discovered class [[58], [K], |135|| of SIPs 



discussed in Sec. 4.2.2, for which the parameters 02 and a\ are related by 
scaling a<i = qa\. What is the special feature of these new potentials that 
interferes with the proof that eq. ( |273|) is exact for SIPs ? To understand this, 
let us look again at the derivation of the lowest order SWKB quantization 
condition. In the derivation, W 2 is taken as 0(%) while hW is 0(h) and 
hence one can expand the integrand on the left hand side in powers of h. This 
assumption is justified for all the standard SIPs [[14], |55| since W 2 is indeed 
of 0(h°) while hW is indeed of 0(h). One might object to this procedure 
since the resulting potential V\ is then ft- dependent. However, in all cases, 
this ^-dependence can be absorbed into some dimensionful parameters in the 
problem. For example, consider 

W(x) = Atanhx , V x (x) = A 2 - A(A + h)sech 2 x. (277) 

Taking A such that A(A + h) is independent of h gives the desired %- 
independent potential (the additive constant is irrelevant and so can contain 
h). Such a move may appear to be of limited value since one cannot apply 
SWKB directly to a superpotential W which is now Ti-dependent. However, 
because A is a free parameter, one can continue the SWKB results obtained 
for A (and hence W) independent of h over to this superpotential and so 
obtain an SWKB approximation for a ^-independent potential. 

What about the new potentials ? In the simplest of these cases, the 
only free parameter in the problem (apart from q) is the combination R\di, 
on which W depends as W(x, Ri^) = \/RiaiF(\J Ridix/h). Incorporating 
different dependences on h in R\a\ will give different ones in W, V\ and E n , 
but F is a sufficiently complicated function that there is no way of eliminating 
h from W 2 . This is a direct consequence of the scaling reparameterisation 
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a2 = qa\ not involving %. If W 2 (x;ai) were independent of H, so would 
W 2 (x; 02) be and in taking the lowest order of the shape invariance condition 
one would get W 2 (x;ai) = W 2 (x;a 2 ), which corresponds to the harmonic 
oscillator. Furthermore, with a 2 = qa\, W 2 and hW are now of a similar 
order in h. The basic distinction between them involved in deriving eq. 
( |273| ) is thus no longer valid and we are prevented from deriving the SWKB 
condition for these new potentials. 

We thus see that the SWKB quantization condition is not the correct 
lowest order formula in the case of the new SIPs and hence it is not re- 
ally surprising that eq. ( |273| ) does not give the exact eigenvalues for these 
potentials. In other words, it remains true that the lowest order SWKB 
quantization condition is exact for SIPs (if the SUSY is unbroken), but only 
in those cases for which the formula is applicable in the first place. It is 
thus still the case that the SIPs given in refs. p4"| , are the only known 
ones for which the lowest order SWKB formula is exact and the higher order 
corrections are all zero. 



6.3 Comparison of the SWKB and WKB Approaches 

Let us now compare the merits of the two schemes [WKB and SWKB]. For 
potentials for which the ground state wave function (and hence the super- 
potential W) is not known, clearly the WKB approach is preferable, since 
one cannot directly make use of the SWKB quantization condition eq. ( |273| ). 
On the other hand, we have already seen that for shape invariant potentials, 
SWKB is clearly superior. An obvious interesting question is to compare 
WKB and SWKB for potentials which are not shape invariant but for whom 
the ground state wave function is known. One choice which readily springs 
to mind is the Ginocchio potential given by [|57| 

V(x) = (1 - y 2 ) j-A 2 ^ + 1) + ^^[2 - (7 - AV + 5(1 - A 2 )?/ 4 ] j 

(278) 

where y is related to the independent variable x by 

^ = (l- 2 / 2 )[l-(l-AV]- (279) 
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Here the parameters v and A measure the depth and shape of the potential 



respectively. The corresponding superpotential is |52| 

W(x) = (l-\ 2 )y(y 2 -l)/2 + ^\ 2 y 
where /z n is given by 



(280) 



/i n A 2 = ^\\\v + 1/2) 2 + (1 - A 2 )(n + 1/2) 2 ] - (n + 1/2) (281) 
and the bound state energies are 



n 2 \ 4 



n 



0,1,2, 



(282) 



For the special case A = 1, one has the Rosen-Morse potential, which is shape 
invariant. The spectra of the Ginocchio potential using both the WKB and 
SWKB quantization conditions have been computed [[TjJ. The results are 
shown in Table 6.1. In general, neither semiclassical method gives the exact 
energy spectrum. The only exception is the shape invariant limit A = 1, 
in which case the SWKB results are exact, as expected. Also, for n — 0, 1 
the SWKB values are consistently better, but there is no clear cut indica- 
tion that SWKB results are always better. This example, as well as other 
potentials studied in ref. |$6, [ 



BjJ , support the conjecture that 
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shape invariance is perhaps a necessary condition for the lowest order SWKB 
approximation to yield the exact spectrum |M| . 

So far, we have concentrated our attention on the energy eigenvalues. 
For completeness, let us note that several authors have also obtained the 
wave functions in the SWKB approximation [79. 91, 92. M. As in the WKB 



method, the SWKB wave functions diverge at the turning points a and b. 
These divergences can be regularized either by the uniform approximation 



511, |155|| or by appropriate retention of higher orders in % ||92"|, |93|. 



6.4 SWKB Quantization Condition for Broken Super- 
symmetry 

The derivation of the lowest order SWKB quantization condition for the case 
of unbroken SUSY is given in Sec. |D] [eqs. Q268j ) to ( [273|) 1. For the case of 
broken SUSY, the same derivation applies until one examines the 0(h) term 
in eq. (|271|) . Here, for broken SUSY, one has 



W{a) = W(b) = VE& (283) 
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and the 0(h) term in ( 271Q exactly vanishes. So, to leading order in h the 
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SWKB quantization condition for broken SUSY is []94 , 

f ^2m[E { n ] - W 2 (x)} dx = (n + l/2)hn, n = 0,1,2, 

J a 



(284) 



As before, it is easy to obtain the quantization condition which includes 
higher orders in h [36] and to test how well the lowest order broken SWKB 
condition works for various specific examples. As for the case of unbroken 
SUSY, it is found that exact spectra are obtained for shape invariant po- 
tentials with broken SUSY PS] . For potentials which are not analytically 
solvable, the results using eq. ( |284j ) are usually better than standard WKB 
computations. Further discussion can be found in ref. |9T> . 



7 Isospectral Hamiltonians 

In this section, we will describe how one can start from any given one- 
dimensional potential V\(x) with n bound states, and use supersymmetric 
quantum mechanics to construct an n-parameter family of strictly isospec- 
tral potentials Vi(Ai, A 2 , . . . , A n ; x) i.e., potentials with eigenvalues, reflection 
and transmission coefficients identical to those for V±(x). The fact that such 
families exist has been known for a long time from the inverse scattering ap- 
proach |67| , but the Gelfand-Levitan approach to finding them is technically 
much more complicated than the supersymmetry approach described here. 
Indeed, the advent of SUSY QM has produced a revival of interest in the de- 
termination of isospectral potentials [^, |69|, |156| , |157| , |158| fl66| , [70], [72], [75 



In Sec. |7.1| we describe how a one parameter isospectral family is obtained 
by first deleting and then re-inserting the ground state of V\(x) using the 
Darboux procedure [Fj3], [7Dfl . The generalization to obtain an n-parameter 



family is described in Sec. |77^ ]72). When applied to a reflectionless potential 
(Sec. |7T3| ), the n-parameter families provide surprisingly simple expressions 
for the pure multi-soliton solutions [73] of the Korteweg-de Vries (KdV) and 



other nonlinear evolution equations [159, 160, 138, 161, 162 



7.1 One Parameter Family of Isospectral Potentials 

In this subsection, we describe two approaches of obtaining the one-parameter 
family Vi(Ai; x) of potentials isospectral to a given potential V\(x). 
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The first approach follows from asking the following question: Suppose 
V2(x) is the SUSY partner potential of the original potential Vi(x), and let 
W(x) be the superpotential such that V 2 {x) = W 2 +W and V 1 (x) = W 2 -W. 
Then, given V 2 (x), is the original potential V\(x) unique i.e., for a given 
V 2 (x), what are the various possible superpotentials W(x) and corresponding 
potentials Vi(x) = W 2 — W'l Let us assume @] that there exists a more 
general superpotential which satisfies 

V 2 (x) = W 2 (x) + W\x). (285) 

Clearly, W = W is one of the solutions to eq. ( |285|) . To find the most general 
solution, let us set 

W(x) = W(x) + (j){x) (286) 
in eq. ( 285 ). We find then that y(x) = <p^ 1 (x) satisfies the Bernoulli equation 

y'(x) = 1 + 2Wy (287) 

whose solution is 

^ = *(«) = ^b^(*)+AJ. (288) 

Here x 

X x {x) = I ipl(x')dx'. (289) 



Ai is a constant of integration and ip\{x) is the normalized ground state wave 
function of V\(x) = W 2 (x) — W'(x). Thus the most general W(x) satisfying 
eq. ( |285| ) is given by 

W{x) =W(x) + 4~ ln[Xi (x) + Ai] , (290) 
ax 

so that the one parameter family of potentials 

d 2 

VAx) = W 2 (x) - W'ix) = VAx) - 2— lnfJ^x) + AJ (291) 

ax 2 

have the same SUSY partner V^x). 

In the second approach, we delete the ground state ipi at energy E\ for 
the potential V\{x). This generates the SUSY partner potential V 2 (x) = 
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V\ — 2-n^ln'0i, which is almost isospectral to Vi(x) i.e., it has the same 
eigenvalues as V\(x) except for the bound state at energy E\. The next step 
is to reinstate a bound state at energy E\. 

Although the potential Vi does not have an eigenenergy Ei, the function 
satifies the Schrodinger equation with potential V2 and energy E\. The 
other linearly independent solution is /f^ ip\ (x')dx' /ipi. Therefore, the most 
general solution of the Schrodinger equation for the potential V 2 at energy 
E 1 is 



$x(Ai) = (Xx + AOM • (292) 

Now, starting with a potential V2, we can again use the standard SUSY 
(Darboux) procedure to add a state at E\ by using the general solution 
$i(Ai), 

U 1 (A 1 ) = U 2 -2^1n$ 1 (A 1 ) . (293) 

The function l/$i(Ai) is the normalizable ground state wave function of 
Vi(Ai), provided that Ai does not lie in the interval — 1 < Ai < 0. Therefore, 
we find a one-parameter family of potentials Vi(Ai) isospectral to V\ for 
Ai > or Ai < — 1. 



Vi(Ai) = V 1 -2^ ln^x^iCAO) 

ax A 

= V 1 -2^1n(J 1 + A 1 ) . (294) 
The corresponding ground state wave functions are 

&(Ai;ar) = l/$i(Ai) . (295) 

Note that this family contains the original potential V\. This corresponds to 
the choices Ai — > ±00. 

To elucidate this discussion, it may be worthwhile to explicitly construct 
the one-parameter family of strictly isospectral potentials corresponding to 
the one dimensional harmonic oscillator 1701. In this case 



W(x) = %-x (296) 
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so that 



Vi(x) = ^a*-% . (297) 



The normalized ground state eigenfunction of V\(x) is 



Mx) = {^) exp(-^74) (298) 

Using eq. (|289|) it is now straightforward to compute the corresponding X\ (x) . 
We get 



2 f°° 2 

Tx(x) = 1 - ^ erfc ( ] ; erfc(x) = -j= J e~ l dt . (299) 




Using eqs. ( |294| ) and ( |295| ), one obtains the one parameter family of 
isospectral potentials and the corresponding wave functions. In Figs. 7.1 
and 7.2, we have plotted some of the potentials and wave functions for the 
case uj = 2. We see that as Ai decreases from oo to 0, V\ starts developing 
a minimum which shifts towards x = — oo. Note that as Ai finally becomes 
zero this attractive potential well is lost and we lose a bound state. The 
remaining potential is called the Pursey potential Vp{x) [36]. The general 
formula for Vp(x) is obtained by putting Ai = in eq. ( f294| ). An analogous 
situation occurs in the limit Ai = —1, the remaining potential being the 
Abraham-Moses potential [ |65|j . 

7.2 Generalization to n-Parameter Isospectral Family 

The second approach discussed in the previous subsection can be generalized 
by first deleting all n bound states of the original potential Vi(x) and then 
reinstating them one at a time. Since one parameter is generated every time 
an eigenstate is reinstated, the final result is a n-parameter isospectral family 



W% - Recall that deleting the eigenenergy E\ gave the potential V 2 (x). The 
ground state ip 2 for the potential V 2 is located at energy E 2 . The procedure 
can be repeated "upward", producing potentials V 3 , V4, . . . with ground states 
^3,^4; • • • at energies E 3 ,E 4 , . . . , until the top potential U n+1 (x) holds no 
bound state (see Fig. 7.3, which corresponds to n = 2). 
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In order to produce a two-parameter family of isospectral potentials, we 
go from V\ to V 2 to V3 by successively deleting the two lowest states of V\ 
and then we re-add the two states at E 2 and E x by SUSY transformations. 
The most general solutions of the Schrodinger equation for the potential V3 
are given by $ 2 (A 2 ) = (X 2 + A 2 )/V2 at energy E 2 , and A 2 &i(\i) at energy 
Ei (see Fig. 7.3). The quantities X i are defined by 

Xi(x) = [ X ^(x')dx' . (300) 
Here the SUSY operator Ai relates solutions of the potentials Vi and 

A i = ^-(]n^ i )' . (301) 
Then, as before, we find an isospectral one-parameter family U 2 (A 2 ), 

U 2 (A 2 ) = V 2 - 2^ ln(J 2 + A 2 ) . (302) 

The solutions of the Schrodinger equation for potentials U 3 and U 2 (A 2 ) are 
related by a new SUSY operator 

A\{\ 2 ) = -A + (i n $ 2 (A 2 ))' . (303) 
Therefore, the solution $i(Ai,A 2 ) at Ex for U 2 (A 2 ) is 

$ 1 (A 1 ,A 2 ) = 4(A 2 )A 2 $i(Ai) . (304) 

The normalizable function l/$i(Ai,A 2 ) is the ground state at E\ of a new 
potential, which results in a two-parameter family of isospectral systems 
U 1 (A 1 ,A 2 ), 



Vi(Ai,A 2 ) = U 1 -2^1n(^ 2 $ 2 (A 2 )<fi(A 1 ,A 2 )) 

= K-2^1n(^i(T 2 + A 2 )$ 1 (A 1 ,A 2 )) , (305) 
for Aj > or \ < — 1. A useful alternative expression is 
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^(Ai, A 2 ) = -V 2 (X 2 ) + 2($' 1 (A 1 , A 2 )/$i(Ai, A 2 )) 2 + 2E 1 . (306) 

The above procedure is best illustrated by the pyramid structure in Fig. 7.3. 
It can be generalized to an n-parameter family of isospectral potentials for an 
initial system with n bound states. The formulas for an n-parameter family 
are 

$;(A;) = (li + Ai)M , i=l,---,n; (307) 
A = ^-(ln^)' I (308) 
it (A,, • ■ - , A n ) = ~ + (In $,(A 4 , ■ ■ • , An))' ; (309) 



^i(Ai, Aj + i, • • • , A n ) 
= A+i(Ai+i, Aj +2 , • • • , \ n )Al +2 (\ i+ 2, A i+3 , • • • , A n ) • • • Ajj(A n ) 
x A n A n _ 1 ---A m $ l (A i ) ; (310) 



d 2 

^(Ai, • • • , A n ) = V X - 2— ln(VlV2 • ■ • ^n*n(A„) • ■ • $l(A 1; • • • , A„)) . (311) 

etar 

The above equations |72[ summarize the main results of this section. 

7.3 n-Soliton Solutions of the KdV Equation 

As an application of isospectral potential families, we consider reflectionless 
potentials of the form 

V x = -n(n + l)sech 2 a; , (312) 

where n is an integer, since these potentials are of special physical inter- 
est. V\ holds n bound states, and we may form a n-parameter family 
of isospectral potentials. We start with the simplest case n = 1. We 
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have V\ = — 2sech 2 x, E\ = — 1 and ijii — 2 asechx. The corresponding 1- 
parameter family is 

Vi(Ai) = -2sech 2 (x+iln(l + ^)) . (313) 

2 Ai 

Clearly, varying the parameter Ai corresponds to translations of V\(x). As 
Ai approaches the limits + (Pursey limit) and — 1~ (Abraham-Moses limit), 
the minimum of the potential moves to — oo and +oo respectively. 

For the case n — 2, V\ — — 6sech 2 x and there are two bound states at 
Ei = —4 and E 2 = —1. The SUSY partner potential is V 2 = — 2sech 2 x. 
The ground state wave functions of Vi and V 2 are ip% = ^sech 2 x and ip 2 = 
-j^sechx. Also, X\ — |(3 tanh x — tanh 3 x + 2) and X 2 = |(tanh x + 1). After 
some algebraic work, we obtain the 2-parameter family 



Vi(Ai,A 2 ) 



12 



[3 + 4 cosh(2x - 25 2 ) + cosh(4a; - 2S{)] 
[cosh(3x — 62 — Si) + 3 cosh(a; + 5 2 — Si)]' 



Si = --Hi 



1 

A~ 



1,2 



As we let Ai — > — 1, a well with one bound state at Ei will move in the +x 
direction leaving behind a shallow well with one bound state at E 2 . The 
movement of the shallow well is essentially controlled by the parameter X 2 . 
Thus, we have the freedom to move either one of the wells. 

It is tedious but straightforward to obtain the result for arbitrary n and 
get Vi(Ai, A2, • ■ • , A n , x). It is well known that one-parameter (t) families of 
isospectral potentials can also be obtained as solutions of a certain class of 



nonlinear evolution equations ||159 
form (q = 0, 1, 2, • ■ ■) 



160| , |138| , |162|| . These equations have the 



-u t = {L u )" u x (314) 
where the operator L u is defined by 

Luf(x) = f xx - Auf + 2u x / dyf(y) (315) 

and u is chosen to vanish at infinity. [For q = we simply get — Ut = u x , 
while for q = 1 we obtain the well studied Korteweg-de Vries (KdV) equa- 
tion]. These equations are also known to possess pure (i.e., reflectionless) 



78 



multisoliton solutions. It is possible to show that by suitably choosing the 
parameters Aj as functions of t in the n-parameter SUSY isospectral family of 
a symmetric reflectionless potential holding n bound states, we can obtain an 
explicit analytic formula for the pure n-soliton solution of each of the above 
evolution equations f75fl . These expressions for the multisoliton solutions of 
eq. ( |314j ) are much simpler than any previously obtained using other proce- 
dures. Nevertheless, rather than displaying the explicit algebraic expressions 
here, we shall simply illustrate the 3-soliton solution of the KdV equation. 
The potentials shown in Fig. 7.4 are all isospectral and reflectionless holding 
bound states at E\ = -25/16, E 2 = —1,-^3 = -16/25 . As t increases, note 
the clear emergence of three independent solitons. 

In this section, we have found n-parameter isospectral families by repeat- 
edly using the supersymmetric Darboux procedure for removing and inserting 
bound states. However, as briefly mentioned in Sec. 7.1, there are two other 
closely related, well established procedures for deleting and adding bound 
states. These are the Abraham-Moses procedure |65| and the Pursey proce- 



dure []66|]. If these alternative procedures are used, one gets new potential 
families all having the same bound state energies but different reflection and 



transmission coefficients. Details can be found in reference \ffl 



8 Path Integrals and Supersymmetry 

In this section, we will describe the Lagrangian formulation of SUSY QM 
and discuss three related path integrals: one for the generating functional 
of correlation functions, one for the Witten index - a topological quantity 
which determines whether SUSY is broken, and one for a related "classical" 
stochastic differential equation, namely the Langevin equation. We will also 
briefly discuss the superspace formalism for SUSY QM. 

Starting from the matrix SUSY Hamiltonian which is 1/2 of our previous 
H [eq. (|60|)1 for convenience: 

H = l -p 2 + \w 2 {x)I - ^}W'{x), 

we obtain the Lagrangian 

L = \-x 2 + i^d t ip - \w 2 {x) + \[ip,^}W'(x). (316) 
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It is most useful to consider the generating functional of correlation functions 
in Euclidean space. We rotate t — > it and obtain for the Euclidean path 
integral: 



Z[j, n, rj*) = J [dx] [#] [#*] exp[-S E + J jx + #* + 77V], (317) 

where 

S E = jTdr {^xl + l -W 2 {x) - r[d T - W'(x)]4) 
and if) and if)* are now elements of a Grassman algebra: 

{p,tl>} = {*i>,ii>} = {P,P} = 0, (318) 

and 

dx 

Xj- — 
dr 

. The Euclidean action is invariant under the following SUSY transformations 
which mix bosonic and fermionic degrees of freedom: 



5x = e*ip + ip*e, 
Sifj* = - e * (d T x + W(x)) , 

Sifj = - e (-d T x + W(x)), (319) 

where e and e* are two infinitesimal anticommuting parameters. These trans- 
formations correspond to iV = 2 supersymmetry. 

The path integral over the fermions can now be explicitly performed using 
a cutoff lattice which is periodic in the the coordinate x but antiperiodic in 
the fermionic degrees of freedom at r = and r = T. Namely we evaluate 
the fermionic path integral: 

[dif)\[dif)*\e£ W\dr-W{*M ( 320 ) 

by calculating the determinant of the operator [d T — W'(x)] using eigenvectors 
which are antiperiodic. 



We have, following Gildener and Patrascioiu ||164 |, that 

det[d T -W'(x)] = l[X m , 
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where 
so that 

Mr) = C m exp[/ dr'[\ m + W']. (321) 
J o 

Imposing the antiperiodic boundary conditions: 

i> m (T) = -^ m (0) 

yields: 

A m = ^?^-if <W' W . (322) 

Regulating the determinant by dividing by the determinant for the case where 
the potential is zero we obtain: 

, r d T -W'(x). , /- T , W'(x) , nnn . 

det — — - — = cosh / dr (323) 

a T Jo 2 

Rewriting the cos/i as a sum of two exponentials we find, as expected that Z 
is the sum of the partition functions for the two pieces of the supersymmetric 
Hamiltonian. Namely when the external sources are zero: 

Tie~ HlT + Tre- H ' 2T = Z^ + Z+. (324) 

For the case when SUSY is unbroken, only the ground state of H 1 con- 
tributes as T — > oo We also have: 

Z ± = J [dx\ssxp[-S%\ (325) 



where 



-L 



T (x 2 



W 2 (x) W'(x)' 



A related path integral is obtained for the noise averaged correlation 
functions coming from a classical stochastic equation, the Langevin equation. 
If we have the stochastic differential equation 

x = W(x(t))+ti(t) (326) 

where 77 (r) is a random stirring forcing obeying Gaussian statistics, then the 
correlation functions of x are exactly the same as the correlation functions 
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obtained from the Euclidean quantum mechanics related to the Hamiltonian 
Hi. To see this we realize that Gaussian noise is described by a probability 
functional: 

P[r,\ = NeM~ fdr^] (327) 
2 Jo t 



normalized so that: 



/ 



/ 



/ 



Dr]P[r]] = 1, 

Dr)P[r)]f(T) = 0, 

Dr ] P[r ] ]f(r)f(r') = F 5(r-r'). 
The correlation functions averaged over the noise are: 

< x(n)x(T 2 )... >= J Dr]P[r]}x(Ti)x(T2)... (328) 

where we have in mind first solving the Langevin equation explicitly for 
x(i](t)) and then averaging over the noise as discussed in Sec. 2.2. An- 
other way to calculate the correlation function is to change variables in the 
functional integral from r\ to x. 

< x(n)x(r 2 )... >= J D[x] P[r)} Det \ |^ | x(n)x(r 2 )... (329) 

This involves calculating the functional determinant, 

subject to the boundary condition that the Green's function obey causality, 
so one has retarded boundary conditions. One has 

Det A = exp / dt Tr In ( [4~ - W'(x(t))]6(t - r') ] dr. (331) 
dx J \ dr J 

When there are no interactions the retarded boundary conditions yield 

G (t-t') = 9(t-t'). (332) 
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Expanding ln(l — GqW) one finds because of the retarded boundary condi- 
tions that only the first term in the expansion contributes so that 



Det 



dx{r') 




drW'(x)}. 



(333) 



Choosing F 



H so that 




= iVexp[-— / dr(x 2 + W 2 (x))\ 



(334) 



we find that the generating functional for the correlation functions is exactly 
the generating functional for the correlation functions for Euclidean quantum 
mechanics corresponding to the Hamiltonian Hi. 



Thus we see that we can determine the correlation functions of x for the 
Hamiltonian Hi by either evaluating the path integral or solving the Langevin 
equation and averaging over Gaussian noise. 

An equation related to the Langevin equation is the Fokker-Planck equa- 
tion, which defines the classical probability function P c for the equal time 
correlation functions of Hi. Defining the noise average: 



Z[j] = N J D[x}ex P [--j o 



dr (x 2 + W 2 (x) - W\x) - 2j(r)x(r))\ (335) 



P c (z) =< 5(z - x(t) >, 



„= J Dr]S(z - x{t))P[n] 



(336) 



one obviously has: 




One can show [|165|1 that P c obeys the Fokker-Planck equation: 




(337) 



S3 



For an equilibrium distribution to exist at long times t one requires that 



P(z,t)^P(z) 

and 

J P{z)dz = 1. 

Setting ^ = in the Fokker-Planck equation, we obtain 

P(z) = Ne- 2 ^ w{y)dy = ^(zf . (338) 

Thus at long times only the ground state wave function contributes (we are in 
Euclidean space) and the probability function is just the usual ground state 
wave function squared. We see from this that when SUSY is broken, one 
cannot define an equilibrium distribution for the classical stochastic system. 

A third path integral for SUSY QM is related to the Witten index. As 
we discussed before, one can introduce a "fermion" number operator via 

n F = = . (339) 

Thus 

(-l) F =[^^]=a 3 . (340) 

The Witten index is given by A = Tr(— 1) F . As we discussed earlier, the 
Witten index needs to be regulated and the regulated index is defined as: 

A 03) = Tr(-l) F e-^ = Tr(e^ Hl - e^ 2 ) . (341) 

In Sec 2.2 we showed how to determine A(/3) using heat kernel methods and 
how it was useful in discussing non-perturbative breaking of SUSY. Here 
we will show that the Witten index can also be obtained using the path 
integral representation of the generating functional of SUSY QM where the 
fermion determinant is now evaluated using periodic boundary conditions to 
incorporate the factor (— 1) F . It is easy to verify a posteriori that this is the 
case. Consider the path integral: 

A(P) = J [dx] [d9] [dV*}eJo ^(*.*.*')* > (342) 
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where 



L E = ~a£ + ^W 2 - ^*[d T - W'(x)\9. 



To incorporate the (— 1) F in the trace, one changes the boundary conditions 
for evaluating the fermion determinant at r = 0, /3 to periodic ones: 



x(0) = x(P) 



*(0) = 



The path integral over the fermions can again be explicitly performed using 
a cutoff lattice which is periodic in the fermionic degrees of freedom at r = 
and t = f3. We now impose these boundary conditions on the determinant 
of the operator [d T — using eigenvectors which are periodic. 

We again have 

det[d T -W'(x)] = Y[\ m . 



Imposing periodic boundary conditions: 

i(2m)7T 



A, 



(3 (3 



(343) 



Regulating the determinant by dividing by the determinant for the case where 
the potential is zero we obtain: 



det 



d T - W'(x) 



sinh 



P , W(x) 
dr 



(344) 



Again rewriting the sinh as a sum of two exponentials we find, as expected 
that we obtain the regulated Witten index: 



A((3) = Z_-Z + = Tie-^ - Tre"^ 2 



(345) 



8.1 Superspace Formulation of Supersymmetric Quan- 
tum Mechanics 

One can think of SUSY QM as a degenerate case of supersymmetric field 
theory in d — 1 in the superspace formalism of Salam and Strathdee [|167| 
(This idea is originally found in unpublished lecture notes of S. Raby ||166|| ). 



superfields are defined on the space (x n ; 9 a ) where x is the space coordinate 
and 9 a are anticommuting spinors. In the degenerate case d — 1 the field 
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is replaced by x(t) so that the only coordinate is time. The anticommuting 
variables are and 0* where 

{0,0*} = {0,0} = [0,t]=O. 

Consider the following SUSY transformation: 

t' = t -i(0*e-e*0), 0' = + e, 0*' = 0* + e*. (346) 

If we assume that finite SUSY transformations can be parametrized by 

then from 

5A = i[e*Q* + Qe,A] (347) 
we infer that the operators Q and Q* are given by: 

Q = id e - 0*d t , Q* = -id e * - 0d t . (348) 

Now these charges obey the familiar SUSY QM algebra: 

{Q, Q*} = 2id t = 2H, [Q, H] = 0. (349) 

The Lagrangian in superspace is determined as follows. A superfield made 
up of x and and 0* can at most be a bilinear in the Grassman variables: 

0( X) e, 0*) = x(t) + i9^(t) - iip*0* + 00*D(t). (350) 

Under a SUSY transformation, the following derivatives are invariant: 

D e = d e - i0*d t 

or in component form: 

D 0( f) = itj}- 0*D - i0*x + 0*0tp, (351) 

and 

Dq* = dg* — iOdt 

or in component form: 

[D e( f)]* = -iijj* -0D + i0x + 0*0ip*. (352) 
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The most general invariant action is: 

S = J dtd8*d6 (-\D e </>\ 2 - /(</>)) (353) 

Again the expansion in terms of the Grassman variables causes a Taylor 
expansion of / to truncate at the second derivative level. Integrating over 
the Grassman degrees of freedom using the usual path integral rules for 
Grassman variables: 



j ode = J e*de* = 1, J de = J de* = o, 



one obtains 



S = J dt(^x 2 + r[dt-f"(x)}^ + ^D 2 + Df'(x)). (354) 

Eliminating the constraint variable D = —f'(x) = W(x) we obtain our 
previous result for the action (now in Minkowski space): 

S = Jdt Qx 2 + V* [d t - W'{x))ip - ^W 2 ^j (355) 



A more complete discussion of this can be found in ref. [ I I 



9 Perturbative Methods for Calculating En- 
ergy Spectra and Wave Functions 

The framework of supersymmetric quantum mechanics has been very useful 
in generating several new perturbative methods for calculating the energy 
spectra and wave functions for one dimensional potentials. Four such meth- 
ods are described in this section. 

In Sees. 9.1 and 9.2, we discuss two approximation methods (the vari- 
ational method and the 5— expansion) for determining the wave functions 
and energ y eigenvalues of the anharmonic oscillator making use of SUSY 
QM. Sec. 9.3 contains a description of a SUSY QM calculation of the energy 
splitting and rate of tunneling in a double well potential. The result is a 
rapidly converging series which is substantially better than the usual WKB 
tunneling formula. Finally, in Sec. 9.4, we describe how the large N expan- 
sion (N = number of spatial dimensions) used in quantum mechanics can be 
further improved by incorporating SUSY. 
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9.1 Variational Approach 

The anharmonic oscillator potential V(x) = gx 4 is not exactly solvable. To 
determine the superpotential one has to first subtract the ground state energy 
E and solve the Riccati equation for W(x): 

Vi(x) = gx 4 ~E = W 2 - W, (356) 

Once the ground state energy and the superpotential is known to some order 
of accuracy, one can then determine the partner potential and its ground 
state wave function approximately. Then, using the SUSY operator 

one can construct the first excited state of the anharmonic oscillator in the 
usual manner. Using the hierarchy of Hamiltonians discussed in Sec. 3, 
one can construct from the approximate ground state wave functions of the 
hierarchy and the approximate superpotentials W n all the excited states of 
the anharmonic oscillator approximately. 

First let us see how this works using a simple variational approach. For 
the original potential, we can determine the optimal Gaussian wave function 
quite easily. Assuming a trial wave function of the form 

^o=(^) 1/4 e-^ 2 , (357) 

71 



we obtain 



(In this subsection, we are taking m — 1 in order to make contact with pub- 
lished numerical results). Minimizing the expectation value of the Hamiltonia 
n with respect to the parameter (5 yields 

£ = (^) 4/ V /3 , /5 = (^) 1/ V /3 . 

This is rather good for this crude approximation since the exact ground state 
energy of the anharmonic oscillator determined numerically is Eq = .668<? 1//3 
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whereas (|) 4//3 = .681 . The approximate potential W resulting from this 



variation calculation is 



W(x) 



dlogj)Q 
dx 



2(3x, 



which leads to a Gaussian approximation to the potential 

V 1G = 4(3x 2 - 2(3. 
The approximate supersymmetric partner potential is now 

V 2G = A(3x 2 + 2f3. 



(359) 



(360) 



(361) 



Since V 2G differs from V XG by a constant, the approximate ground state wave 
function for V 2 is given by eq. ([357]). The approximate ground state energy 
of the second potential is now 



< H 2 >= v < 0\^L + V 2G \0 > v = H, + 4/3. 



(362) 



Thus we have approximately that the energy difference between the ground 
state and first excited state of the anharmonic oscillator is 

E l -E =A(3 = 4( 3 -)^g 1 /\ 
The approximate (unnormalized) first excited state wave function is 



(!) = d _ 2(3x ^ ^ _ A(3xe -^ 
dx 



(363) 



This method can be used to find all the excited state wave functions and 
energy levels of the anharmonic oscillator by using the methods discussed in 



ref. 1163 



Let us now look at a more general class of trial wave functions. If we 
choose for the trial ground state wave functions of the hierarchy of Hamilto- 
nians, 



4 fc) = iV fc exp[-l(gP 



N k = 2 X /7J7X(1 



-1/2 



(364) 
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we obtain much better agreement for the low lying eigenvalues and eigen- 
functions. It is convenient in this case to first scale the Hamiltonian for the 
anharmonic oscillator, 

(365, 

by letting x — *> x/g 1 ^ 6 and H — > g 1 ^H. Then we find the ground state 
energy of the anharmonic oscillator and the variational parameters p\ and 
ri\ by forming the functional 

1 d 2 

^o(pi,nO = (Vol - 2dx^ + x ^ o) - (366) 
Thus we first determine p\ and n\ by requiring 

gf = °> S = °- (367) 
The equation for the energy functional for the anharmonic oscillator is 

% ""'»-^-r(ir + ' , 'f(i)' (368) 

Minimizing this expression, we obtain the following variational result: 

E = 0.66933, n x = 1.18346, p x = 0.666721 . (369) 

This ground state energy is to be compared with a numerical evaluation 
which yie Ids 0.667986 . Since the trial wave function for all ground states is 
given by eq. Q364] ), the variational superpotential for all k is 

W kv = n k \x\ 2n *- l {p k )- nk . (370) 

Since we are interested in the energy differences E n — E n -\ of the anharmonic 
oscillator, we consider the variational Hamiltonian 

H vk +i = \A kv A\ v (371) 

which approximately determines these energy differences. We obtain the 
approximate energy splittings by minimizing the energy functional 

SE k (p k ,n k ) = \{^ k+1) \ - ^ + W 2 k + WU\4 vk+1) > ■ (372) 
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Performing the integrals one obtains the simple recursion relation: 




) 




(373) 



One can perform the minimization in p analytically leaving one minimization 
to perform numerically 

The results for the variational parameters and for the energy differences 
are presented in Table 9.1 for the first three energy eigenvalues and compared 
with a numerical calculation, based on a shooting method. 



9.2 5 Expansion Method 

In this method [|7£| we consider the anharmonic oscillator as an analytic 
continuation from the harmonic oscillator in the paramater controlling the 
anharmonicity. That is we consider simultaneously potentials of the form 



where M is a scale parameter, 5 measures the anharmonicity, and C is the 
ground state energy of the anharmonic oscilator. C is subtracted as usual 
from the potential so it can be factorized. The standard anharmonic oscil- 
lator corresponds to S = 1 and M = (2(/) 1 / 3 . To approximately determine 
W{x) from V\(x) we assume that both W(x) and V\(x) have a Taylor series 
expansion in 5. Thus we write: 



where E n corresponds to the order Taylor expansion of the dependence of 
the ground state energy on the parameter 5. We assume 




(374) 




(375) 



oo 




(376) 



n=0 
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and insert these expressions in eq. ([374]) and match terms order by order. 
At lowest order in S the problem reduces to the supersymmetric harmonic 
oscillator. We have: 

W 2 -W^ = M 2 x 2 - 2E , (377) 

whose solution is 

W (x) = Mx, E = (378) 

To next order we have the differential equation: 

- 2W 1 W = -M 2 x 2 ln(Mx 2 ) + 2E 1 (379) 



dx 

which is to be solved with the boundary condition W n (0) = 0. The order S 
contribution to the energy eigenvalue E\ is determined by requiring that the 
ground state wave function be square integrable. Solving for W\ we obtain 

Wi(x) = -e Mx2 f dye~ My2 \M 2 y 2 \n(My 2 ) - 2E X ] . (380) 
Jo 

To first order in 5 the ground state wave function is now : 

Mx)=e- Mx \l-5 ! X dyW x {y)) 



Imposing the condition that ipo vanishes at infinity, we obtain: 

Ex = ^M^(3/2) , i){x) = r'(x)/r(x). (381) 

Writing M = (2(7) 1//3 , we find that the first two terms in the 5 expansion for 
the ground state energy are 

E = \{2gy*[l + ~V(3/2)5]. (382) 

At 5 = 1, we get 

E = M15g 1/3 . (383) 

A more accurate determination of the ground state energy can be obtained 
by calculating up to order 5 2 and then analytically continuing in S using Pade 
approximants. This is discussed in ref. ||169| . 
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9.3 Supersymmetry and Double Well Potentials 

Supersymmetric quantum mechanics has been profitably used to obtain a 
novel perturbation expansion for the probability of tunneling in a double 



well potential Since double wells are widely used in many areas of 

physics and chemistry, this expansion has found many applications ranging 
from condensed matter physics to the computation of chemical reaction rates 



TO [Tfll , [T72| , [1731, p74| , TO |17§ ,|77], |l7l, 0, H HI- HOI HH • In what 



follows, we shall restrict our attention to symmetric double wells, although 



an extension to asymmetric double wells is relatively straightforward ||181 |. 

Usually, in most applications the quantity of interest is the energy differ- 
ence t = Ei — E between the lowest two eigenstates, and corresponds to the 
tunneling rate through the double-well barrier. The quantity t is often small 
and difficult to calculate numerically, especially when the potential barrier 
between the two wells is large. Here, we show how SUSY facilitates the 
evaluation of t. Indeed, using the supersymmetric partner potential V^x), 
we obtain a systematic, highly convergent perturbation expansion for the 
energy difference t. The leading term is more accurate than the standard 
WKB tunneling formula, and the magnitude of the nonleading terms gives a 
reliable handle on the accuracy of the result. 

First, we briefly review the standard approach for determining t in the 
case of a symmetric, one- dimensional double- well potential, Vi(x), whose 
minima are located x = ±Xo- We define the depth, D, of V\{x) by D = 
Vi(0) — Vi(xq). An example of such a potential is shown in Fig. 9.1. For 
sufficiently deep wells, the double-well structure produces closely spaced pairs 
of energy levels lying below Vi(0). The number of such pairs, n, can be 
crudely estimated from the standard WKB bound-state formula applied to 
V\(x) for x > 0: 



nir 



o 



: [T4(0) -Vxix^dx, (384) 



where x c is the classical turning point corresponding to energy Vi(0) and we 
have chosen units where h = 2m = 1. We shall call a double- well potential 
"shallow" if it can hold at most one pair of bound states, i.e., n < 1. In 
contrast, a "deep" potential refers to n > 2. 

The energy splitting t of the lowest-lying pair of states can be obtained 



by a standard argument ||126|| . Let x( x ) be the normalized eigenfunction for 
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a particle moving in a single well whose structure is the same as the right- 
hand well of V\{x) (i.e., x > 0). If the probability of barrier penetration is 
small, the lowest two eigenfunctions of the double- well potential Vi(x) are 
well approximated by 

4%) = ix(x)±x(-x)]/V2 . (385) 
By integration of Schrodinger's equation for the above eigenfunctions, it can 



be shown that ||126 ] 



t = E 1 -E = 4 X (0)x'(0) , (386) 

where the prime denotes differentiation with respect to x. This result is 
accurate for "deep" potentials, but becomes progressively worse as the depth 
decreases. Use of WKB wave functions in eq. ( |386| ) yields the standard 
result: 

%kb = {^"(xo)] 1 / 2 /^} exp (-2 [°[Vi{x) - V^xo^dx^j . (387) 



The same result can also be obtained via instanton techniques | |182 



Using the supersymmetric formulation of quantum mechanics for a given 
Hamiltonian, Hi = —d 2 /dx 2 + V\{x), and its zero-energy ground state wave 
function ipo(x), we know that the supersymmetric partner potential V 2 (x) is 
given by 

V 2 (x) = V 1 (x)-2(d/dx)(^' /^ ) 

= -V l {x)+2^'J^f . (388) 

Alternatively, in terms of the superpotential W(x) given by W(x) = 
—ip' /ipo we can write 

V 2 ,x{x) = W 2 (x) ± dW/dx . (389) 

From the discussion of unbroken SUSY in previous sections, we know that 
the energy spectra of the potentials V 2 and V\ are identical, except for the 



ground state of V\ which is missing from the spectrum of V 2 |10| . Hence, 
for the double- well problem, we see that if V\(x) is "shallow" (i.e., only the 
lowest two states are paired), then the spectrum of V 2 is well separated. In 
this case, V 2 is relatively structureless and simpler than V\. Previous papers 



31, 32, R3j have implicitly treated just the case of shallow potentials, and, 
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not surprisingly, have found that the use of SUSY simplifies the evaluation 
of the energy difference t. In contrast, let us now consider the case of a 
deep double well as shown in Fig. 9.1 . Here, the spectrum of V 2 has a 
single unpaired ground state followed by paired excited states. In order to 
produce this spectrum, V 2 has a double-well structure together with a sharp 
"<5— function like" dip at x — 0. This central dip produces the unpaired 
ground state, and becomes sharper as the potential V\{x) becomes deeper. 

As a concrete example, we consider the class of potentials whose ground 
state wave function is the sum of two Gaussians, centered around ±xo, 

i/) (x) ~ e - {x - xo)2 + e - {x+xo)2 . (390) 

The variables x and Xq have been chosen to be dimensionless. The corre- 
sponding superpotential W(x), and the two supersymmetric partner poten- 
tials Vi(x) and V^x), are given respectively by 

W(x) = 2[x - x Q tanh(2xx )], (391) 



V 2 ,i(x) = 4[x - x tanh(2xx )] 2 ± 2[1 - 2£ 2 1 sech 2 (2xx )] . (392) 

The minima of V\(x) are located near ±xo and the well depth (in the 
limit of large x ) is D ~ 4xq. We illustrate the potentials V\(x) and V 2 (x) in 
Fig. 9.2 for the two choices xq = 1.0 and xq = 2.5. We see that in the limit 
of large xq, for both Vi(x) and V^x), the wells become widely separated and 
deep and that V^z) develops a strong central dip. 

The asymptotic behavior of the energy splitting, t, in the limit Xq — > 00 
can be calculated from eq. (|386|) , with xi x ) given by one of the (normalized) 
Gaussians in eq. (|390| ). We find that 

t -> 8a:o(2/7r) 1/2 e- 2:c o . (393) 
The same result can be obtained by observing that V\(x) — > 4(| x J — ^o) 



as Xq — > 00. This potential has a well known ||183|| analytic solution, which 
involves solving the parabolic cylindrical differential equation. After carefully 
handling the boundary conditions, one obtains the separation of the lowest 
two energy levels to be 8xq(2/7i) 1 ^ 2 exp(— 2xq), in agreement with eq. (|393|) . 
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We now turn to the evaluation of t via the ground state energy of the 
super symmetric partner potential V^x). In general, since V^(x) is not ana- 
lytically solvable, we must solve an approximate problem and calculate the 
corrections perturbatively. The use of SUSY, coupled with the observation 
that the magnitude of t is in general small, allows us to construct a suit- 
able unperturbed problem. Consider the Schrodinger equation for V^x) and 
E = 0. From supersymmetry [eq. ( |388| )1 we see immediately that 1/"0q is a 
solution. Since t is small, we expect this solution to be an excellent approx- 
imation to the correct eigenfunction for small values of x. However, 1/i/jq 
is not normalizable and hence is not acceptable as a starting point for per- 
turbation theory. One possibility is to regularize the behavior artificially at 
large | x | pi| . This procedure is cumbersome and results in perturbation 



corrections to the leading term which are substantial. Instead, we choose 
for our unperturbed problem the second linearly independent solution of the 
Schrodinger equation given by ||184 



1 r°° 

<j,(x) = — / ^(x')dx', x > , (394) 

and <j)(x) = 4>(—x) for x < 0. Clearly, <f>(x) is well behaved at x — ±oo and 
closely approximates 1/Vo at small x; thus we expect it to be an excellent 
approximation of the exact ground state wave function of V^(x) for all values 
of x. The derivative of <f)(x) is continuous except at the origin, where, unlike 
the exact solution, it has a discontinuity | + —<$> |o-= — 2^0 (0) ■ Hence cf>(x) 
is actually a zero-energy solution of the Schrodinger equation for a potential 
Vo(x) given by 

^o(x) = V r 2 (x)-4^ 2 (0)5(x) , (395) 

where we have assumed that ipo(x) is normalized. We calculate the perturba- 
tive corrections to the ground state energy using AV = +4-0Q (0) 6 (x) as the 
perturbation. Note that the coefficient multiplying the 5-function is quite 
small so that we expect our perturbation series to converge rapidly. 

For the case of a symmetric potential such as V^x), the perturbative 
corrections to the energy arising from AV can be most simply calcualted 
by use of the logarithmic perturbation-theory [|185|| formulation of the usual 
Rayleigh- Schrodinger series. The first and second order corrections to the 
unperturbed energy E = are 
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2£(0) 



POO 

£(ar) = / <f) 2 {x')dx, £ (2) = -2 

J X 



POO 
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(p(x') 



2 



dx' 



(396) 

For our example, we numerically evaluate these corrections in order to 
obtain an estimate of t. The results are shown in Fig. 9.3 for values of 
Xq < 2. Estimates of t correct to first, second, and third order calculated 
from logarithmic perturbation theory are compared with the exact result for 
V2, obtained by the Runge-Kutta method. The asymptotic behavior of t 
given by eq. ( |393| ) is also shown. This asymptotic form can also be recovered 
from eq. (|396|) by a suitable approximation of the integrand in the large-Xo 
limit. Even for values of xq < 1/V2, in which case V\(x) does not exhibit a 
double-well structure, the approximation technique is surprisingly good. The 
third-order perturbative result and the exact result are indistinguishable for 
all values of xq. 

In conclusion, we have demonstrated how SUSY can be used to calculate 
t, the energy splitting for a double-well potential. Rather than calculating 
this splitting as a difference between the lowest- lying two states of Vi(x), 
one can instead develop a perturbation series for the ground state energy t 
of the partner potential Vi^x). By choosing as an unperturbed problem the 
potential whose solution is the normalizable zero-energy solution of V2(x), 
we obtain a very simple 5— function perturbation which produces a rapidly 
convergent series for t [Q. The procedure is quite general and is applicable 



to any arbitrary double-well potential, including asymmetric ones [181]. The 



numerical results are very accurate for both deep and shallow potentials. 



9.4 Supersymmetry and Large-iV Expansions 



The large- N method, where N is the number of spatial dimensions, is a pow- 
erful technique for analytically determining the eigenstates of the Schrodinger 
equation, even for potentials which have no small coupling constant and hence 
not amenable to treatment by standard perturbation theory [|186|, |187|, |188 , 
Pq , [19li |192t p3[ , Pi p5| , p00| . A slightly modified, physically motivated 
approach, called the "shifted large-iV method" H, [193, HI HH 

incorpo- 
rates exactly known analytic results into 1/N expansions, greatly enhancing 
their accuracy, simplicity and range of applicability. In this subsection, we 
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will descibe how the rate of convergence of shifted 1/N expansions can be 
still further improved by using the ideas of SUSY QM [[77] . 



The basic idea in obtaining a 1/N expansion in quantum mechanics con- 
sists of solving the Schrodinger equation in N spatial dimensions, assuming 
TV to be large, and taking 1/N as an "artificially created" expansion parame- 
ter for doing standard perturbation theory. At the end of the calculation, one 
sets N = 3 to get results for problems of physical interest in three dimensions. 

For an arbitrary spherically symmetric potential V(r) in N dimensions, 
the radial Schrodinger equation contains the effective potential 

V eS (r) = V(r) + i*- 1 ^-^ 2 ; k = N + 21 . (397) 

It is important to note that N and I always appear together in the combi- 
nation k = N + 21. This means that the eigenstates, which could in principle 
have depended on the three quantities N,l,n, in fact only depend on k and 
n, where n is the radial quantum number which can take values 0,1,2,... One 
now makes a systematic expansion of eigenstates in the parameter 1/k, where 
k = k — a. Of course, for very large values of N, the two choices k and k 
are equivalent. However, for N = 3 dimensions, a properly chosen shift a 
produces great improvements in accuracy and simplicity. At small values of 
r, the n = wave function ipo(r) has the behavior r^' 1 ^ 2 . If one sets 

^ (r) = r (fc - 1)/2 $ (r) (398) 

where $o( r ) is finite at the origin, then eq. ( |398| ) readily gives the supersym- 
metric partner potential of V e g(r) to be 

, , (k + l)(k- l)h 2 h 2 d 2 , . . , , 

V 2 (r) = V(r) + ± Q n 2 ; T^ ln ®° r 399 

Vzfr) and V e s(r) have the same energy values [except for the ground 
state]. However, large- iV expansions with the partner potential V^(r) are 
considerably better since the angular momentum barrier in eq. ( |399| ) is given 
by (k' — l)(k' — 3)h 2 /8mr 2 , where k' = k + 2. So, effectively, one is working 
in two extra spatial dimensions! Thus, for example, in order to calculate the 
energy of the state with quantum numbers k,n of V c s(r) one can equally 
well use k' = k + 2, n — 1 with V^r). To demonstrate this procedure, let us 
give an explicit example. Using the usual choice of units h = 2m = 1, the 
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s-wave Hulthen effective potential in three dimensions and its ground state 
wave function are: 



V£{r) = -i^—r r + ^~~~r~^~ , Mr) ~ (1 - e~ dr ) , (400) 



1 - e~ 5r 4 



where the parameter 5 is restricted to be less than 2. The supersymmetric 
partner potential is 

2S 2 p~ Sr 

v " = v " + jr^pw ■ < 401 > 

As r tends to zero, goes like 2r~ 2 , which as mentioned above, corre- 
sponds to the angular momentum barrier (k' — l){k' — 3)h /8mr 2 for k' = 5 
(N = 5,1 = 0). Let us compute the energy of the first excited state of 
Vgjj(r). For the choice 5 = 0.05, the exact answer is known to be 0.748125 
201||. The results upto leading, second and third order using a shifted 1/N 



expansion for Vg are 0.747713, 0.748127 and 0.748125 . The corresponding 
values using the supersymmetric partner potential Vf 1 are all 0.748125 ! It is 
clear that although excellent results are obtained with the use of the shifted 
1/N expansion for the original potential Vjg(r) in three dimensions, even 
faster convergence is obtained by using the supersymmetric partner poten- 
tial, since we are now effectively working in five dimensions instead of three. 
Thus, SUSY has played an important role in making a very good expansion 
even better | 77| . In fact, for many applications, considerable analytic simpli- 



fication occurs since it is sufficient to just use the leading term in the shifted 



1/N expansion for V^r). Other examples can be found in ref. [77 . 



10 Pauli Equation and Supersymmetry 

So far, we have discussed the concept of SUSY for Schrodinger Hamiltonians 
in one-dimension and for central potentials in higher dimensions which are 
essentially again one-dimensional problems. In this section we shall show that 
the concept of SUSY can also be applied to some geniune two-dimensional 
problems. In particular we show that the Pauli Hamiltonian which deals with 
the problem of a charged particle in external magnetic field can always be put 
in SUSY form provided the gyromagnetic ratio is equal to two. It is worth 
emphasizing here that not only the uniform magnetic field problem (i.e. the 
famous Landau level problem) but the nonuniform magnetic field problems 
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can also be put in SUSY form. Using the concepts of SUSY and shape 
invariance we show that some of the nonuniform magnetic field problems can 
be solved analytically. 

The Pauli Hamiltonian for the motion of a charged particle in external 
magnetic field in two dimensions is given by (h = 2m = e = 1) 

H={ Px + A x f + {p y + A y ) 2 + |(V x A) z a z (402) 
It is easily seen that this H along with the supercharges Q 1 and Q 2 defined 

by©! 

Q 1 = -^[-(p y + A y )a x + (p x + A x )a y ] 

Q 2 = ^[(Px + A x )a x + ( Py + A y )a y ] (403) 

satisfy the N = 1 supersymmetry algebra provided the gyromagnetic ratio g 
is two 

{Q a , Q p } = H5 a ^ [H, Q a ] = 0; a, (3 = 1, 2 (404) 

It is interesting to note that SUSY fixes the value of g. The above Hamil- 
tonian has an additional 0(2) ® 0(2) symmetry coming from <j z and an 0(2) 
rotation in the A 1 , A 2 plane (A 1 = p x + A x ; A 2 = p y + A y ). Let us first 
consider the problem in an asymmetric gauge i.e. 

A y (x,y) = 0,A x (x,y) = W(y) (405) 

where W(y) is an arbitrary function of y. In this case the Pauli Hamiltonian 
takes the form 

H = (p x + W(y)) 2 +p 2 y - W\y)a z (406) 

Since this H does not depend on x, hence the eigenfunction ip can be factor- 
ized as 

i>(x,y) = e lk ^(y) (407) 

where k is the eigenvalue of the operator p x (— oo < k < oo). The Schrodinger 
equation for ip(y) then takes the form 

iJl + (W(y) + k) 2 - W'{y)a]iP{y) = E^{y) (408) 



dy 
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where a{= ±1) is the eigenvalue of the operator a z . Thus we have reduced 
the problem to that of SUSY in one dimension with superpotential W(y) + k 
where W(y) must be independent of k. This constraint on W(y) strongly 
restricts the allowed forms of shape invariant W(y) for which the spectrum 
can be written down algebraically. In particular from Table 4.1 we find that 
the only allowed forms are (i) W(y) = ui c y + c\ (ii) W(y) = atemhy + c\ 
(iii) W(y) = at&ny + C\ (iv) W(y) — C\ — c 2 exp (— y) for which W(y) can 
b e written in terms of simple functions and for which the spectrum can be 
written down algebraically |202| , 101 . In particular when 



W{y) = to c y + c a (409) 

which corresponds to the uniform magnetic field, then the energy eigenvalues 
known as Landau levels are given by 

E n = {2n + 1 + a)u c ; n = 0,1,2... (410) 

Note that the ground state and all excited states are infinite-fold degenerate 
since E n does not depend on k which assumes a continuous sequence of value 
(— oo < k < oo). 

The magnetic field corresponding to the other choices of W are (ii) B = 
—asech 2 y (iii)£> = — asec 2 y(^ L < y < |) (iv) B = +c 2 exp (— y) and as 
mentioned above, all these problems can be solved algebrically by using the 
results of Sec. 4. 

Let us now consider the same problem in the symmetric gauge. We choose 

A x = u c yf(p), A y = -u c xf(p) (411) 

where p 2 = x 2 + y 2 and uj c is a constant. The corresponding magnetic field 
B z is then given by 

B z (x, y) = d x A y ~ d y A x = -2uJ(j>) - u c pf{p) (412) 
In this case the Hamiltonian ( [402| ) can be shown to take the form 

H = "(^2 + J*) + ~ 2UJ J L Z ~ ( 2C °J + UcPf')<Tz (413) 
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where L z is the z-component of the orbital angular momentum operator. 
Clearly the corresponding Schrodinger problem can be solved in the cylindri- 
cal coordinates p, 0. In this case, the eigenfunction ip(p, 0) can be factorized 
as 

^(p,0)= J R(p)e im Vv / P (414) 

where m — 0, ±1, ±2, .. is the eigenvalue of L z . In this case the Schrodinger 
equation for R(p) takes the form 

d 2 m 2 — 1 I A 

+ ulp 2 f - 2oo Jm - (2ojJ + uj c pf(p))a + ^ 1 ]R(p) = ER(p) 

P ' (415) 

where a(= ±1) is the eigenvalue of the operator o z . On comparing with 
Table 4.1, it is easily checked that there is only one shape invariant potential 
(f(p) = 1) for which the spectrum can be written down algebraically. This 
case again corresponds to the famous Landau level problem i.e. it corresponds 
to the motion of a charged particle in the x — y plane and subjected to a 
uniform magnetic field (in the symmetric gauge) in the z-direction. The 
energy eigenvalues are 

E n = 2(n + m+ | m \)u c ;n = 0, 1,2... (416) 

so that all the states are again infinite-fold degenerate. It is worth noting 
that the nonuniform magnetic fields can also give this equi-spaced spectrum 



[ |203|| . However, they do so only for one particular value of m while for other 
values of m, the spectrum is in general not equi-spaced. 

The fact that in this example there are infinite number of degenerate 
ground states with zero energy can be understood from the Aharonov-Casher 
theorem ||204|| which states that if the total flux defined by $ = / B z dxdy = 



n + e(0 < e < 1) then there are precisely n — 1 zero energy states. Note that 
in our case $ is infinite. 



11 Supersymmetry and the Dirac Equation 

There have been many applications of SUSY QM in the context of the Dirac 
equation. In view of the limitations of space, we shall concentrate on only a 
few of these applications ||100|| . In particular, we discuss the supersymmetric 



structure of the Dirac Hamiltonian (Hd) and show how the methods used 
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to obtain analytical solutions of the Schrodinger equation can be extended 
to the Dirac case. First of all, we consider the Dirac equation in 1+1 di- 
mensions with Lorentz scalar potential <p(x). We show that whenever the 
one-dimensional Schrodinger equation is analytically solvable for a potential 
V(x), then there always exists a corresponding Dirac scalar potential prob- 
lem which is also analytically solvable ||101|| . It turns out that, on the one 



hand, 4>{x) is essentially the superpotential of the Schrodinger problem and 
on the other hand, it can be looked upon as the kink solution of a scalar 
field theory in 1+1 dimensions. Next we discuss the celebrated problem of 
the Dirac particle in a Coulomb field [ |103|| and show that its eigenvalues and 
eigenfunctions can be simply obtained by using the concept of SUSY and 
shape invariance as developed in Sees. 2 and 4. We also discuss the problem 
of the Dirac equation in an external magnetic field in two dimensions and 
show that there is always a supersymmetry in the problem in the massless 
case. We also classify a number of magnetic field problems whose solutions 
can be algebraically obtained by using the concepts of SUSY and shape in- 
variance ||101|| . In addition, we show that the Euclidean Dirac operator in 
four dimensions, in the background of gauge fields, can always be cast in 
the language of SUSY QM. Finally, we discuss the path integral formulation 
of the fermion propagator in an external field and show how the previously 
known results for the constant external field can be very easily obtained by 



using the ideas of SUSY plj . 



11.1 Dirac Equation With Lorentz Scalar Potential 

The Dirac Lagrangian in 1+1 dimensions with a Lorentz scalar potential 
4>{x) is given by 

C = vfyfdyft - injxp. (417) 

The scalar potential 4>(x) can be looked upon as the static, finite energy, kink 
solution corresponding to the scalar field Lagrangian 

= \d^{x)d»<P{x) - V{4>). (418) 
Such models have proved quite useful in the context of the phenomenon 



of fermion number fractionalization 1 205 , 206, 207 1 which has been seen in 



certain polymers like polyacetylene. Further, a variant of this model is also 
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relevant in the context of supersymmetric field theories in 1+1 dimensions 



208| , |209|| . Note that the coupling constant in eqs. (|417|) and ( |41S| ) has been 
absorbed in and V(<f>) respectively. 

The Dirac equation following from eq. ( [417D is 



irfd^x, t) - <j>{x)ij){x, t) = 0. (419) 

Let 

tl>(x,t) = exp(-iut)ip(x) (420) 
so that the Dirac equation reduces to 

7 °^(x) + t^^l - <f>(xWx) = 0. (421) 
dx 

We choose 

v=^=(! i).v=^=(j _° i ).^)=C:g) <-) 

so that we have the coupled equations 

Aipi(x) = uj%1) 2 (x) 

A^ 2 (x) = oj^(x), (423) 

where 

A = ^- + </>{x), A* = -4- + <P(x). (424) 
dx dx 

We can now easily decouple these equations. We get 

A*Aiffi = lu 2 iPu AA^ 2 = to 2 ip 2 . (425) 

On comparing with the formalism of Sec. 2, we see that there is a supersym- 
metry in the problem and <p(x) is just the superpotential of the Schrodinger 
formalism. Further ipi and ip2 are the eigenfunctions of the Hamiltonians 
H- = A + A and H + = AA^ respectively with the corresponding potentials 
being V T (x) = <j) 2 (x) =F <j)'(x). The spectrum of the two Hamiltonians is thus 
degenerate except that H_(H + ) has an extra state at zero energy so long 
as (p(x — > ±oo) have opposite signs and <p(x — > +oo) > 0(< 0). Using the 
results of the Sees. 2 and 4 we then conclude that for every SIP given in 
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Table 4.1 there exists an analytically solvable Dirac problem with the corre- 
sponding scalar potential (j)(x) being the superpotential of the Schrodinger 
problem. In particular, using the reflectionless superpotential given by 

W(x) = ntanhx (426) 

one can immediately construct perfectly transparent Dirac potentials with 
n bound states ||210|| . Further, using the results for the SIP with scaling 



ansatz (a 2 = qa\) |58], p9|l , one can also construct perfectly transparent Dirac 



potentials with an infinite number of bound states. 

11.2 Supersymmetry and the Dirac Particle in a Coulomb 
Field 

The Dirac equation for a charged particle in an electromagnetic field is given 
by (e = % = c — 1) 

[iY(d^ + iA^) -m]i(> = (427) 
For a central field i.e. A = and Aq(x, t) = V(r), this equation can be 



written as ||211 



i^- = Ei\) = (&p + (3m + V)ip (428) 
where 

with a 1 being the Pauli matrices. For central fields, this Dirac equation can 
be separated in spherical coordinates and finally for such purposes as the 
computation of the energy levels one only needs to concentrate on the radial 
equations which are given by [pi 1 



k n 

G'(r) + — -(a 1 -V)F = 
r 

kF 

F'(r) (a 2 + V)G = (430) 

r 

where 

cti = m + E, a 2 = m — E (431) 
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and Gk is the "large" component in the non-relativistic limit. Of course 
the radial functions G k and F k must be multiplied by the appropriate two 
component angular eigenfunctions to make up the full four- component solu- 



tions of the Dirac equation ||211|| . These coupled equations are in general not 
analytically solvable; one of the few exceptions being the case of the Dirac 
particle in a Coulomb field for which 

V(r) = --,7 = Ze 2 (432) 
r 

We now show that the Coulomb problem can be solved algebraically by 
using the ideas of SUSY and shape invariance. To that end, we first note 
that in the case of the Coulomb potential, the coupled equations ( |430| j can 
be written in a matrix form as 

G' k (r)\ , 1 fk - 7 \ (G k \ _ ( ai \ (G k 



where k is an eigenvalue of the operator — (a-L + 1) with the allowed values 
k = ±1, ±2, ±3..., and satisfies | k \— J+\- Following Sukumar |103| , we now 
notice that the matrix multiplying 1/r can be diagonalized by multiplying it 
by a matrix D from the left and D~ x from the right where 

k + s —7 



«=U 7 * + .;■ s = V* 2 -^. (434) 

On multiplying eq. (|433|) from the left by the matrix D and introducing the 
new variable p = Er leads to the pair of equations 

A'G = + («5) 



where 



■E S' 



G f)= d ( G f ) (13G) 



and 

d s 7 1+ d s 7 

A=- + ±Ai = -- + -. 437 

dp p s dp p s 



106 



Thus we can easily decouple the equations for F and G thereby obtaining 

H-F = A+AF = £~)F, 

H + G = AA + G = 4-5c. (438) 

S 

We thus see that there is a supersymmetry in the problem and H± are shape 
invariant supersymmetric partner potentials since 

2 2 

H + (p; s, 7 ) = s + 1, 7 ) + ^ - T^TTvr ( 439 ) 

On comparing with the formalism of Sec. 4 it is then clear that in this case 

7 2 7 2 

02 = s + 1, ai = s, R(a 2 ) 77 77 (440) 

so that the energy eigenvalues of H_ are given by 

^-# s ^=g fl w = ^-^- (441) 

Thus the Coulomb bound state energy eigenvalues E n are given by 

Tfl 

E n = 3 -,n = 0,1,2, ... (442) 

1 _|_ 7 2 1 1/2 ' ' ' ' V ' 

L 1 (s+n) 2 J 

It should be noted that every eigenvalue of H_ is also an eigenvalue of H + 
except for the ground state of if_ which satisfies 

AF = =^ F (p) = p s exp{--fp/s) (443) 

Using the formalism for the SIP as developed in Sec. 4, one can also alge- 
braically obtain all the eigenfunctions of F and G. 

Notice that the spectrum as given by eq. ( [442j ) only depends on | k | 
leading to a doublet of states corresponding to k —\ k | and k = — | k | for 
all positive n. However, for n — 0, only the negative value of is allowed 
and hence this is a singlet state. 
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11.3 SUSY and the Dirac Particle in a Magnetic Field 

Let us again consider the Dirac equation in an electromagnetic field as given 
by eq. (|427|) but now consider the other case when the vector potential is 



nonzero but the scalar potential is zero i.e. A = 0, A ^ 0. Now as shown 
by Feynman and Gell-Mann |[212|| and Brown | 213 |, the solution of the four 



component Dirac equation in the presence of an external electromagnetic 
field can be generated from the solution of a two component relativistically 
invariant equation. In particular, if ip obeys the two component equation 

[(P + A) 2 + m 2 + a-(B + iE)]if> = (E + A ) 2 ip (AAA) 

then the four component spinors that are solutions of the massive Dirac 
equation are generated from the two component ip via 

^-[(aiP + ^+E-Ao-m)^)- (445) 

Thus, in order to solve the Dirac equation, it is sufficient to solve the much 
simpler two-component eq.( [444j ) and then generate the corresponding Dirac 
solutions by the use of eq. ( |445| ) . In the special case when the scalar po- 
tential Aq (and hence E) vanishes, the two-component equation then has 
the canonical form of the Pauli equation describing the motion of a charged 
particle in an external magnetic field. If further, m = and the motion is 



confined to two dimensions, then the Pauli eq. ( [444] ) exactly reduces to the 
eq. ( |402j ) of the last section. Further, since 

(H D ) 2 = [a-(P + A)} 2 = H Pauli (446) 

hence, there is a supersymmetry in the massless Dirac problem in external 
magnetic fields in two dimensions since Hp, Qi and Q 2 (see eq. ( |403|) satisfy 
the SUSY algebra as given by eq. ( |404j ) . Clearly, this supersymmetry will 
also be there in two Euclidean dimensions. We can now immediately borrow 
all the results of the last section. In particular, it follows that if the total flux 
$(= / B z dxdy)=n + e(0 < e < 1) then there are precisely n — 1 zero modes 
of the massless Dirac equation in two dimensions in the background of the 
external magnetic field B (B = B z ) ||204| , 214|| . Further, in view of eqs. ( [44^ 



and ( 445|) we can immediately write down the exact solution of the massless 



Dirac equation in an external magnetic field in two dimensions in all the four 
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situations discussed in the last section when the gauge potential depended 
on only one coordinate (say y). Further using the results of that section, 
one can also algebraically obtain the exact solution of the Dirac equation in 
an uniform magnetic field in the symmetric gauge when the gauge potential 
depends on both x and y. 

Even though there is no SUSY, exact solutions of the Pauli and hence 
the Dirac equation are also possible in the massive case. On comparing 
the equations as given by (|444| ) (with A = 0, E = 0) and ( |402|) it is clear 
that the exact solutions in the massive case are simply obtained from the 
massless case by replacing E by E — m 2 . Summarizing, we conclude that 
the exact solutions of the massive (as well as the massless) Dirac equation 
in external magnetic field in two dimensions can be obtained algebraically 
in case the magnetic field B(= Bz) has any one of the following four forms 
(i) B = constant (ii) B = —a sech 2 y (iii) B = — a sec 2 y(— n/2 < y < ti/2) 
(iv) B = — C2exp(— y). Further, in the uniform magnetic field case, the 
solution can be obtained either in the asymmetric or in the symmetric gauge 



20^, [1011. 



11.4 SUSY and the Euclidean Dirac Operator 

In the last few years there has been a renewed interest in understanding 
the zero modes and the complete energy spectrum of the square of the Dirac 
operator for a Euclidean massless fermionic theory interacting with the back- 
ground gauge fields [^, |215| , |216|| . Let us first show that there is always 



a supersymmetry in the problem of the Euclidean massless Dirac operator in 
the background of the gauge fields. On defining 

Q±=^(l±75)7^ (447) 

where 

D^ = d^ + (448) 
one finds that the Dirac operator can be written as 

= Q+ + Q_ (449) 
Let us take the following representation of the Euclidean 7 matrices 

o)^=(° o)^=(i -1) ( 45 °) 
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where a 1 are the usual Pauli matrices. It is then easily seen that the operators 
Q+, Q- and H = —( r y^D fJi ) 2 satisfy the usual N = 1 SUSY algebra 

H = -(^D,) 2 = {Q+,Q_}, [H,Q+] = = [#,Q_] (451) 

This supersymmetry is popularly known as chiral supersymmetry. In view 
of the above representation of the 7 matrices, we find 

iD + a-3 \ _ ( L 

and hence 



V>* = i- w : +s .n =U» 0) < 452 > 



tt I r. X2 r. r. ( &-{B + E) \ / 



(453) 



where we have used the convention 



Ffj, v = d^A v — d u A^, Bi = -€ijkFj k , E t = F 0i (454) 

Using the techniques of Sees. 2 and 4, one can algebraically obtain the 
eigenvalues and the eigenfunctions of H for few special cases ||101|| . Further, 
once we have an eigenfunction of H with a nonzero eigenvalue E±, then one 
can easily obtain the eigenfunctions <f>± of 7 /J -t> /i with eigenvalues ±(E 1 ) 1 ^ 2 
by the construction 

<f )± = ( 1 ,D,±(E 1 ) 1 ^ (455) 

Finally, it is worth pointing out that in case the gauge potential has 
the special form 

4* = f^d vX (456) 
where the constant 4 by 4 matrix has the properties: 

f jw fvixi f iivfva $fia (457) 

then there exists another breakup of the Dirac operator 7 M -D^ which leads to 
the so called complex supersymmetry [ |101|| . 
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11.5 Path Integral Formulation of the Fermion Prop- 
agator 

We have seen that there is always a supersymmetry associated with the 
Euclidean Dirac operator in four as well as two dimensions in the background 
of the gauge fields. This SUSY was first successfully exploited in the context 
of the study of the chiral anomalies ||217|| . In a nut shell, the solvability of the 
Dirac operator is related to the ability to integrate the path integral. The 
quantities one wishes to calculate are the Green's functions 

G(x,x;t) =< x | e~ Hr \ x' >\ x=x > (458) 

and 

G 5 (x, x- r)=<x\ Tr( l5 e- HT ) | x' >\ x=x , (459) 

By using Schwinger's proper time formalism we can determine S(x, x; A) 
from G. Similarly, the index Z5 which in the limit r = is related to the 
chiral anomaly, is just the spatial integral over G 5 . 

Here we would like to show that by introducing the fermionic degrees 
of freedom, the path integral for the Dirac operator, which is initially a 
path ordered integral gets reduced to an ordinary path integral. This trick 
was first introduced by Rajeev ||218|| who was interested in reformulating 



quantum electrodynamics as a supersymmetric theory of loops. Once we 
introduce the fermionic variables then G(G§) is determined by choosing the 
antiperiodic (periodic) boundary conditions for the fermions. By integrating 
over the fermionic degrees of freedom, a purely bosonic path integral is then 
obtained. For the case of a constant external field strength F^, the bosonic 
path integral is a Gaussian, which allows one to trivially obtain G and G$. 

Let us consider the square of the Dirac operator i.e. (7^-D^) 2 . It can also 
be written as 

H =(P + A) 2 -^F tlu (460) 

where 

^ = ^,71- (461) 

If we do not introduce auxiliary fermions then the related matrix valued 
Lagrangian would be 

L(x, x) = -x^Xf, - iA^x^ + \^ v F»v (462) 
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and we would obtain for Feynman's path integral representation 



< x | e \ x >= / Px A1 (r)Pexp[- / dr'L(x,x)] 



(463) 



where P denotes path ordering. Now, taking analogy from SUSY quantum 
mechanics, it has been noticed [[218| , |101|| that one can introduce the Grass- 
man variables ip^ via 



Then H can be written as 



H = (P + A) 2 + ~^F^ U 



1 i 



and hence the Lagrangian now becomes 

which is invariant under the SUSY transformations 

5 X/t = -ieip^Stp 
We now obtain for the path integral 



(464) 



(465) 



(466) 



(467) 



G(x, x\ r) 



J Vx„(t)Vip„exp[- J T dr'L ss (x,x)), (468) 



where we impose antiperiodic boundary conditions on the fermion at and r. 
Since the fermionic path integral is quadratic, we can perform the functional 
integral over the fermionic degrees of freedom exactly for arbitrary F^. 
The result of the fermionic path integral is 

Det ' /2 <;f- + f ->, (469) 

where the prime denotes the omission of the zero mode. To evaluate the 
determinant one puts F^ v in skew diagonal form: 



' 







" 


—X\ 




















X2 








-%2 


. 



(470) 
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and imposes either periodic or antiperiodic boundary conditions. One obtains 
for the fermion determinant using antiperiodic boundary conditions relevant 
forG 

JJcosh-^ Xidr. (471) 

and for the fermion determinant with periodic boundary conditions relevant 
for the determination of the index: 

(-rnt-shmj f Xidr']. (472) 

Thus for any arbitrary external field one can explicitly perform the fermionic 
path integral and is left with the purely bosonic path integral which, for ex- 
ample, could by performed numerically by Monte Carlo techniques. For the 
particular case of a constant external field one can further explicitly perform 
the resulting bosonic path integral. When 



A„ = -Ax v , (473) 



with F^v being a constant matrix, one obtains for the remaining bosonic path 
integral: 

TT- Xi (474) 
\ L 2 sinh 2f 1 ; 

G and G§ are easy to calculate in two and four Euclidean dimensions for 
the constant external field case ||101|| using this method. In two dimensions, 



where only F i = B exists one has that {x\ = B) 101 



G = Btanh(5r); G 5 = B = -.e^F^ (475) 



It is worth remarking here that the same result has been obtained with 
more difficulty by Akhoury and Comtet Hl6| . In four dimensions one has 
a simpler derivation of the result of Schwinger ||219|| rewritten in Euclidean 
space. One finds that 

Xi = [F + (F 2 - G 2 ) 1/2 ] 1/2 , x 2 = [F - (F 2 - G 2 ) 1/2 ] 1/2 , (476) 

where 

F = -F F 
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so that 



From this we obtain that: 



G = -F F* 



X\ Xi = \f, u F; u . (477) 



2 cosh^ 21 



G(x,x;r)=F, u F; u Y[^i 



G 5 (x,x;r) = -F flu F; u . (478) 



8=1 ^u-"" 2 

1 

T 

Once one has obtained the Green's function it is easy to reconstruct the 
effective action which allows one to obtain the rate of pair production from 
a strong external electric field. The effective action is given by 

S eff = J dxL{x), L(x) = J™ dss- l e- m2s G{s) (479) 

12 Singular Superpotentials 

So far, we have only considered nonsingular superpotentials W(x) which give 
rise to the supersymmetric partner potentials V%(x) and V^x). This choice 
of the superpotential was based on the ground state wave function ipo(x) of 
V\(x) by the relation W(x) = —i/;' (x) / i/; (x) . In this section, we describe a 
general procedure for constructing all possible superpotentials which yield a 
given potential V(x) upto an additive constant [|114|| . This general procedure 
is based on any arbitrary solution 4>(x) of the Schrodinger equation for V(x), 
rather than just the ground state wave function. We shall see that singular 
superpotentials are given by = —<p'/<p and the singularities are located 
at the zeros of excited wave functions. Such singularities produce interesting 
properties for the potentials Vu^x) and V^^x) obtained via the Riccati 
equations Vi^)(x) = W 2 (x) - W'(x) and V 2 ^)(x) = W 2 (x) + W'(x). Singular 
superpotentials are responsible for negative energy eigenstates for Vi^(x) 
and often give rise to a breakdown of the degeneracy of energy levels for 
Vi(<t>)( x ) an d V2( ( j,)(x) ||110| , |112| , |1 13| , |220| , |114|| . Another interesting physical 



phenomenon occurs if one considers a singular superpotential resulting from 



a solution <j)(x) for an energy E in the classical energy continuum. ||114|| The 
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isospectral family of V\^) (x) is then found to have bound states (normalised 
eigenfunctions) at energy E. Thus SUSY QM provides a systematic pro- 
cedure for generating potentials possessing the purely quantum mechanical 
phenomenon of bound states in the continuum ||116| , |lTI 



12.1 General Formalism, Negative Energy States and 
Breakdown of the Degeneracy Theorem 

Given any nonsingular potential V(x) with eigenfunctions ip n {x) and eigen- 
values E n (n = 0,1,2,...), let us now enquire how one can find the most 
general superpotential W(x) which will give V(x) upto an additive constant 



114j| . To answer this question consider the Schrodinger equation for V(x): 

- 0" + V{x)(j> = e<j) (480) 

where e is a constant energy to be chosen later. For convenience, and without 
loss of generality, we will always choose a solution <p(x) of eq. ( [48 Oj ) which 
vanishes at x — — oo. Note that whenever e corresponds to one of the eigen- 
values E n , the solution <f>(x) is the eigenfunction ip n (x). If one defines the 
quantity = —4>'/4> and takes it to be the superpotential, then clearly the 
partner potentials generated by W<f, are 

V m = Wl + , V m = Wl-W' 4> = < ^ = V(x) - e , (481) 

where we have used eq. ( |480| ) for the last step. The eigenvalues of Vi(^) are 
therefore given by 

E nW =E n -e. (482) 

One usually takes e to be the ground state energy E and to be the ground 
state wave function ipo(x), which makes £Wo) = and gives the familiar case 
of unbroken SUSY. With this choice, the superpotential Wq{x) = —ip' /ijjQ 
is nonsingular, since ipo{x) is normalizable and has no zeros. The partner 
potential V^) has no eigenstate at zero energy since A q i/jq(x) = [d/dx + 
Wq(x)]%I)q(x) = 0; however, the remaining eigenvalues of V2U) are degenerate 
with those of Viuy 

Let us now consider what happens for other choices of e, both below and 
above the ground state energy Eq. For e < Eq, the solution <p(x) has no 
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nodes, and has the same sign for the entire range — oo < x < +00. The 
corresponding superpotential W^{x) is nonsingular. Hence the eigenvalue 
spectra of and V^) are completely degenerate and the energy eigenval- 
ues are given by eq. ( |482|) . In particular, Ew^) = E — e is positive. Here, 



has the same sign at z = ±00, and we have the well-studied case of broken 



SUSY 10 



If the constant e is chosen in the range E < e < Ei, then the solution 
4>(x) of eq. ( |482|) will have only one node, say at the point x = x s . Near the 
point x = x s , the node of 4>(x) makes the superpotential W(x) singular 

(j) = a {x-x s ), <p' = a, W{x) = -{x - x s )'\ (483) 

and the partner potentials have the behavior [using eq. ( [48 1[ ) ] 



V m = 0, V m = 2(x-x s y 2 . (484) 



It is well-known [|221| , [222| , |223|1 that for any singular potential V(x) = 
X(x — x s )~ 2 , the behavior of the wave function at x = x s is governed by the 
value of A. For values A > |, the potential V(x) has a "strong" singularity 
which forces ip(x s ) = and the range of x is effectively broken into two 
disjoint pieces x < x s and x > x s with no communication between them. For 
— 7 < A < |, the potential has a singularity of "intermediate" strength. It is 
not strong enough to make if)(x s ) vanish and in fact the two regions x < x s 
and x > x s do communicate and one has the full range — 00 < x < +00. 
Here, in principle, one can have singular wave functions which require self- 



adjoint extensions ||221|| ; in what follows, we will deal with regular solutions. 
Finally, for A < — 7, the Hamiltonian is unbounded from below. 

From eq. ( |484j ) , it is clear that V^) has a "strong" singularity which 
makes i/)(x s ) = 0. Thus, the problem of finding the eigenstates of Vzu) is 
really two separate problems; one in the range — 00 < x s and the other in 
the range x s < x < +00. Clearly, this is very different from the range of 
Viu) which is the whole real axis —00 < x < +00. Hence, in general, the 
degeneracy theorem obtained from SUSY for the spectra of V^) and \4(</>) 
is not valid. The above discussion can be readily extended to all values of 
the constant e - the only difference being the number of poles present. A 
summary of results is given in Table 12.1. 

So far we have considered superpotentials (both nonsingular and singular) 
which give rise to a nonsingular potential Vi(^). However, if the potential 
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V(x) itself has singularities, then so must all superpotentials which produce 
it. Consider the case of a simple pole singularity at x = xq. Then, near 
x = x , the singular superpotential W(x) = g(x — x ) _1 gives the following 
behavior to the corresponding partner potentials: 

g(g + i) 9(9-1) 

Then clearly, for g > | or g < — |, both Vi(^) and V^) have "strong" 
singularities and both have two disjoint regions x < xo and x < Xq. Here, 
one expects degenerate energy levels (corresponding to broken or unbroken 
SUSY). Similarly, for — | < g < +|, both Vi(^) and V^) have "intermediate" 
strength singularities and the whole region — oo < x < +00 is valid for 
both. Here again, one obtains degeneracy. However, for the two regions 
— I < g < — ~ and | < g < | only one of the partner potentials has a 
"strong" singularity, whereas the other has a singularity of "intermediate" 
strength. Therefore, the two potentials have different Hilbert spaces and, 
in general, degeneracy is not there. The above discussion is borne out by 



the work of Jevicki and Rodrigues ||110|| who consider a superpotential of the 
form W(x) = g/x — x. 

The general discussion of singular superpotentials, negative energy states 
and breakdown of the degeneracy theorem is best illustrated with some spe- 
cific examples. 

Consider the harmonic oscillator potential V(x) = x 2 . The energy eigen- 
values are E n = 2n and the first few eigenstates are 

Using the solution = ipo(x), one gets the usual nonsingular superpotential 
Wo = —ipot/ipa = x, leading to the following partner potentials and eigenval- 
ues: 

U 1(0) =x 2 -1, U 2(0) =x 2 + 1, £$, } = 2n, E® Q) = 2n + 2, n = 0,1,2,... 

(487) 

This is standard unbroken SUSY. The harmonic oscillator potential can also 
be obtained if one starts from the solution <ft = ipi (x) . One gets the singular 
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superpotential W\ = —ipi/ipi — x — 1/x. It has a pole at x = and is a 
special case of the form discussed in Ref. ||110|| . The partner potentials are 



V m = x 2 - 3, V 2(1) = x 2 + 4 - 1, (488) 

Ob 



and both are exactly solvable [ 54|, [54|. The two eigenvalue spectra are 



sW ) = 2n-2, ^=4/1 + 4, 71 = 0,1,2,... (489) 

There occurs a negative energy state in Vim at ^o(i) = — ^. This is expected 

since we chose to be the first excited state and fi'jl = 0, which pushes 
the ground state to a negative energy. Proceeding along the same lines, we 
can get the harmonic oscillator potential using yet another superpotential. 
Taking (f> = ip2( x ), one gets W 2 = —ip' 2 /ip2 = x — 4x/(2x 2 — 1), which has 
poles at x = ±1/a/2. The corresponding potentials are 

9 „ 8(2x 2 + l) . . 

V m =x 2 -5, V m =x 2 -3+ ( ^ 2 _ 1)2 - (490) 

The eigenvalue spectrum of Vi( 2 ) is E^Ls = 2n — 4 indicating two negative 
energy states at —4 and —2, which was anticipated. The potential V 2 (2) is 
not analytically solvable. The "strong" singularities at x — ±l/v2 break 
the x-axis into three disjoint regions — oo < x < — l/y/2, — l/v2 < x < 
l/\/2, l/\/2 < x < +oo and the degeneracy theorem breaks down. 

In summary, the potentials Vi( ) and V^o) give a realization of SUSY QM 
with a non-degenerate zero energy ground state and pairing of excited states. 
For the case of Vim and V 2 (i) the eigenvalues are given by eq. (|489| ) . There 



is partial degeneracy of the spectrum [see Fig. 12.1] due to the fact that the 
potentials are symmetric. The states which are missing in V^n) are the even 
parity states, since the 2/x 2 barrier requires the vanishing of wave functions 
at x = 0. Finally, for the partner potentials Vi( 2 ) and V 2 ( 2 ), degeneracy is 
completely absent. 

Our second example is the Morse potential V(x) = A 2 + B 2 e~ 2ax — B(a + 
2A)e~ ax . The superpotential based on the ground state is W (x) = A — 
Be~ ax . For concreteness, we shall take A = 4, a = 1 and B = 1. The 
corresponding energy eigenvalues are E^h. = 16 — (4 — n) 2 ; there are four 
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bound states with eigenvalues 0, 7, 12 and 15 and a continuum starting above 
16. The two lowest eigenfunctions are |54 



Mx) ~ e~ {ix+e ~ x \ ^i(x) ~ e~ i3x+e ~ x) (7 - 2 e - ;r ). (491) 

The supersymmetric partner potentials constructed from the nonsingular su- 
perpotential Wq{x) = 4 — e~ x are 

V m = 16 - 9e~* + e~ 2x , V m = 16 - 7e~ x + e" 2x . (492) 

The spectrum of tfyo) is identical to that of Vi(o), except that there is no 
state at zero energy. An alternative singular superpotential W\{x) which 
also yields the Morse potential is 

= _ £ = 21-15e- + 2e- 
u ; 7-2 e -* V ' 

There is a simple pole at x s = — ln(3.5). The partner potentials are 

n n T 2t Tr 441 - 567e- x + 28le 2x - 56e~ 3x + Ae~ 4x 
V m =9-9e- x + e- 2x , V m = (7 _ 2e _ x)2 

(494) 

By construction, the potential V^n) is the Morse potential and V^m has a 
singularly at x s = — ln(3.5). Expanding V%(i) about the singular point gives 
^2(1) ~ 2/(x — x s ) 2 , which requires the wave function to vanish at x — x s . 
This effectively breaks the potential and the real axis into two parts: Vfyi) 
(left) for — oo < x < x s and V 2 (i) (right) for x s < x < +oo. The potentials 
and their eigenstates are plotted in Figure 12.2. The potential Vi(i) has 
energy levels located at -7, 0, 5, 8 and a continuum above 9. As expected 
from Table 12.1, there is a negative energy state at -7. We have calculated the 
energy levels of Ifyi) (right) and Vfyi) (left) numerically. They are E}J (right) 

= 6.08, 8.65 with a continuum above 9 and E$ (left) = 22.96, 50.69, 82.16, 
116.8,.... The potential VWi) (left) has an infinite number of bound states 
and as is obvious from the spectra, the degeneracy between V 2 (i) and is 
completely broken by the "strong" singularity in W\(x). Note that since the 
Morse potential is asymmetric, no partial degeneracy remains. This is unlike 
our first example (harmonic oscillator) which had symmetric potentials and 
VWi) and Vi(i) had a partial degeneracy. 
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Our last example is the reflectionless potential V(x) = A 2 — A(A + 
l)sech 2 x. The motivation for considering this example is to make contact 
with singular superpotentials considered by Casahorran and Nam [fL12[| . The 
ground state wave function is given by = (sechx) A . The nonsingular su- 
perpotential Wq = Atanhx comes from Using the property of shape 
invariance |44| one readily obtains ip\(x) = tanhx(sechx) A_1 as the wave 
function for the first excited state. If one takes <f> = ipi, the resulting su- 
perpotential is W\ = — V'i/V'i = — 2Acosech 2x + (A — l)cothx, which 



agrees with eq. (3.14) of ref. ||112|| . The energy eigenvalues for Vi(o) are 



^n(o) = A 2 — (A — n) 2 and one has unbroken SUSY. The potentials generated 
by W x {x) are V m = (A - l) 2 - A(A + l)sech 2 x, V 2{1) = (A - l) 2 - A(A - 
l)sech 2 x + 2cosech 2 a;. These potentials are shape invariant |54], and their 
eigenvalues are = (A-l) 2 -(A-n) 2 , E$ l} = (A-l) 2 -(A-2n-3) 2 , in- 
dicating partial degeneracy, as expected for symmetric potentials. Similarly, 
our previous experience indicates that the partner potentials constructed 
from W2(x) will have no degeneracy. Thus, we have not only reproduced 
some families of singular superpotentials previously considered in the litera- 
ture, but given the general method to construct new ones. 

In conclusion, we have shown how excited state wave functions can be 
used to construct singular superpotentials in SUSY QM. Our method pro- 
vides a complete and unified picture of the origin of negative energy states 
and the presence or absence of degeneracy. Although our technique is per- 
fectly general, our examples were taken from the class of shape invariant 
potentials, since these are analytically solvable. 

12.2 Bound States in the Continuum 



In 1929, Von Neumann and Wigner ||116|| realized that it was possible to 



construct potentials which have quantum mechanical bound states embed- 
ded in the classical energy continuum (BICs). Further developments, by 
many authors | |117| |224| , |225| , |226|| have produced more examples and a bet- 



ter understanding of the kind of potential that can have such bound states, 
although there is not as yet a fully systematic approach. These authors 
have also suggested possible applications to atoms and molecules. Capasso 
et al. |[227|| have recently reported direct evidence for BICs by constructing 



suitable potentials using semiconductor heterostructures grown by molecular 
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beam epitaxy. Finally, it is interesting to note that BICs have found their 



way into a recently written text 118 



In this subsection, we show how one can start from a potential with 
a continuum of energy eigenstates, and use the methods of SUSY QM to 
generate families of potentials with bound states in the continuum [BICs] 
[ |115|| . Basically, one is using the technique of generating isospectral potentials 
(discussed in Sec. 7) but this time starting from states in the continuum. The 
method preserves the spectrum of the original potential except it adds these 
discrete BICs at selected energies. As illustrative examples, we compute and 
graph potentials which have bound states in the continuum starting from a 
null potential representing a free particle and the Coulomb potential. 

(a) One Parameter Family of BICs 

Consider any spherically symmetric potential V(r) which vanishes as r — > 
oo. The radial s-wave Schrodinger equation for the reduced wave function 
u(r) (in units where H = 2m = 1) is 

- u" + V(r) u{r) = E u{r), (495) 

where we have scaled the energy and radial variables such that all quantities 
are dimensionless. Eq. ([495 ) has a classical continuum of positive energy 



solutions which are clearly not normalizable. 

As we have seen in Sec. 7, the Darboux |64| procedure for deleting and 



then reinstating the ground state uo(r) of a potential V(r), generates a fam- 
ily of potentials V(r; A) which have the same eigenvalues as V(r). These 
isospectral potentials are labeled by a real parameter A in the ranges A > 
or A < —1. The isospectral potential V(r; A) is given in terms of the orig- 
inal potential V(r) and the original ground state wave function uo{r) by 

V(r; A) = V{r) - 2[ln(7 + A)]" = V{r) - ±^ + (496) 

i + A (I + X) 2 



where 



J Q (r) = ^ u 2 (r')dr'. (497) 



Except in this section, uq was taken to be the nodeless, normalizable 
ground state wave function of the starting potential V(r). However, it is easy 
to generalize the above equations to the case where uq{t) is any solution of 
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eq. ( f495| ) with arbitrary energy Eq. If uq(t) has nodes, this leads to singular 
superpotentials and to singularities in the partner potential V^r). However, 
when the original state at E is re-inserted, the resulting family of potentials 



V(r; A) is free of singularities ||114|| . Our results are best summarized in the 
following statement: 

Theorem: Let Uq(t) and U\(r) be any two nonsingular solutions of the 
Schrodinger equation for the potential V(r) corresponding to arbitrarily se- 
lected energies E and Ei respectively. Construct a new potential V(r; A) as 
prescribed by eq. (|496|) . Then, the two functions 



u (r;X) = ^{, (498) 

and 

u x {r] A) = {Eh - E ) Ul + u W{u , ui), (499) 

[where W denotes the Wronskian, W(uo,ui) = u u[ — uiu' Q ] are solutions of 
the Schrodinger equation for the new potential V(r; A) corresponding to the 
same energies E and E\. 

While the new potential in eq. ( |496p and the new wave functions in eq. 



( [49 8|) were originally inspired by SUSY QM, the easiest proof of the above 
theorem is by direct substitution. One simply computes —u" + V(r; A)«i 
(i=0,l), with the wave functions Ui given in the theorem. After straightfor- 
ward but tedious algebraic manipulations, one gets EiUi, thus establishing 
the theorem. The algebra is considerably simplified by using the following 
identity for the Wronskian of two solutions of the Schrodinger equation: 

4-W{uq, u x ) = {E - E 1 )u u 1 . (500) 
ar 

Let us now take Uq to be a scattering solution at a positive energy Eq = k 2 
of a potential V(r) which vanishes at r=oo. Taking uq(v = 0) = satisfies 
one of the required boundary conditions, but clearly u oscillates as r — » oo 
and has an amplitude which does not decrease. Consequently, the integral 
Io(r) in eq. (f497| ) now grows like r at large r and uq is now square integrable 
for A > 0, while the original wave function uq was not. Negative values of A 
are no longer allowed. Therefore, we see that all the potentials V(r; A) have 
a BIC with energy Eq. Note from eq. (|498|) that uq has the same zeros as 
the original uq. At zeros of Uq, V(r; A) and V(r) are equal. All the other 
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oscillatory solutions of the Schrodinger equation with V(r) get transformed 
into oscillatory solutions to the new Schrodinger equation with V(r; A) with 
the same energy. In particular, note that Ui(r; A) remains a non-normalizable 
scattering solution of the corresponding Schrodinger equation. 

We note that the new potential V(r; A) in eq. (|496|) and the BIC at energy 



Eq are formed using the corresponding wave function Uq{t). Any other state, 
say U\(r), is transformed into a solution of the new Schrodinger equation 



by the operation given in eq. (|499 ) which involves both uq and U\. The 
central column of Table 12.2 gives a convenient overview of the relationship 
of the potentials V and V and the solutions of the corresponding Schrodinger 
equations. 

We now give two examples to explicitly illustrate how one applies the 
above procedure to obtain potentials possessing one BIC. 

Example 1: Free Particle on the Half Line. Consider a free particle 
on the half line (V = for < r < oo). We choose Uq = sin kr, the spherical 
wave solution, corresponding to energy Eq = k 2 > 0, which vanishes at r = 0. 
The integral Iq given in eq. ( |497| ) becomes 

I = [2kr - trin(2Jfer)]/(4Jfe). (501) 

We observe that / ^ >"/2 as r ^ oo. 

The potential family V, defined in eq.( [496| ) becomes 



32^sm^_8^sm(2^) 
El En 



with 

D (r; A) = 2Jfer - sin(2Jfer) + 4JfeA. (503) 
V has a BIC at energy Eq = k 2 with wave function 

tto(A) = Ak sinkr/E . (504) 

For special values of the parameters k and A, the potential V and its 
BIC wave functions are shown in Figs. 12.3a and 12.3b. The original null 
potential has now become an oscillatory potential which asymptotically has 
a 1/r envelope. The new wave function at Eq = k 2 also has an additional 
damping factor of 1/r which makes it square integrable. As uq appears in the 
numerator of V, eq. ( |496| ), every node of uq is associated with a node of V 
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but not every node of V produces a node of uq. The value of the eigenenergy 
Eq is clearly above the asymptotic value, zero, of the potential. Evidently, 
the many oscillations of this potential, none of them able to hold a bound 
state, conspire in such a way as to keep the particle trapped. 

The parameter A which appears in the denominator function D Q (r; A) 
plays the role of a damping distance; its magnitude indicates the value of r 
at which the monotonically growing integral Iq becomes a significant damping 
factor, both for the new potential and for the new wave function. This is 
illustrated graphically in Figs. 12.3a and 12.3b which are drawn for very 
different values of A. [Note that the wave functions shown in the figures are 
not normalized]. The parameter A must be restricted to values greater than 
zero in order to avoid infinities in V and in the wave functions. In the limit 
A — > oo, V becomes identical to V. 

Example 2: Coulomb Potential. Here V — Z/r, and the unbound, 
reduced / = wave function satisfies the Schrodinger equation eq. ( [495| ) , which 
can be written in standard form 

Uq + (1 - 2^)u = (505) 

with p = \f~E~r and f\ = Z/2\/~E. 

The solutions involve confluent hypergeometric functions which in the 
asymptotic limit approach sine waves phase-shifted by a logarithmic term. 
Useful expressions for these solutions in the regions near and far from the 
origin are available in the literature ||228 , |229|| . Stillinger and Herrick ||117 



following the method of Von Neumann and Wigner [ 11611 , have constructed 
BIC potentials and wave functions for the case of the repulsive Coulomb 
potential. Here we use our theorem to construct a one-parameter family of 
isospectral potentials containing a BIC The procedure is the same for both 
positive and negative Z; the only difference being in the sign of rj. The formal 
expressions for the BIC potentials and wave functions have been given above, 
eqs. (|496|) and ( |498| ), in terms of u$. 

The positive energy solution of eq. (|505|) can be written in the usual form 



228| , |229| , [1171 as the real function 

u (p) = C (fj) e-^M(l - if}, 2, 2ip), (506) 

where 

C (f l ) = (e-^ 2 )\T(l + tf l ) | (507) 
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and M(a, b; z) is Rummer's function. Using tabulated expressions for the 
Coulomb wave functions [229] and doing the integral for Jo numerically, we 
have obtained the BIC wave functions for representative values of A. The 
corresponding one-parameter family of potentials obtained by the SUSY pro- 
cedure is given in eq. ( f496| ) with V = Z/r. 

The results are displayed in Fig. 12.4. Fig. 12.4a shows the BIC partner 
to the attractive Coulomb potential for A = 1, k — 1, and Z = -2. Fig. 12.4b 
shows the (unnormalized) wave function of the bound state in the continuum 
for this potential at Eq = k 2 . For comparison the original Coulomb potential 
and wave function are also shown dotted. It is seen that the potential which 
holds a bound state of positive energy shows an oscillatory behavior about 
the Coulomb potential, Vq, as is also evident from the form of eq. ( |49(j| ) for 
V. Since the oscillating component vanishes whenever w vanishes, we have 
V = V at each node of Uq. Compared to the original, unnormalizable wave 
function, the BIC wave function in both cases shows a damped behavior due 
to the denominator function. This is also seen in the figures. 

A similar behavior is also expected for other spherically symmetric po- 
tentials with a continuous spectrum of positive eigenvalues. For arbitrary 
one-dimensional potentials, where the range extends from — oo to +oo, the 
situation is not so clear cut. Our method works for the Morse potential which 
is steeply rising on the negative x-axis with correspondingly damped wave 
functions. It also works for the case of a particle in a constant electric field 
for similar reasons. For potentials, such as V(x) = —Vq sech 2 x, the integral 
J in eq. ([49 7|) is not convergent if the starting point is chosen at — oo, and it 
gets negative contributions if the starting point is selected at finite x- values. 
This leads to a vanishing denominator function in the expressions for some 
wave functions which makes them unacceptable. 

(b) Two Parameter Family of Potentials 

In the previous section (a), we have seen how to generate a one-parameter 
family of potentials with one BIC. We now show how this procedure can be 
extended to construct two-parameter families which contain two BICs. 

In constructing the new wave functions for the one-parameter family, eq. 
( [49 6|) , we observe that the denominator function given in eq. ([498|) was all 
that was needed to create the BIC, while the operation in eq. (|499|) ensured 
that the wave functions for all the other states, there represented by u\, are 
a solution to the new potential. Note again, there is nothing special about 
the ordering of the two energy values nor the relative magnitude of Eq and 
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Ei, therefore we can repeat this procedure by applying the theorem to the 
wave functions and the potential of the one-parameter family, but this time 
we transform the state at E\ into a BIC. The state at E , which already is 
a BIC, is transformed in the step of eq. ( |499| ), suitably modified, to become 
a solution to the new potential. In this way we obtain the two parameter 
family of potentials 



V(r;X,X 1 ) = V-2[ln(h + X 1 )]' , = V 



h + Ai (A + \\ 



with the solutions of the corresponding Schrodinger equation 

M = {Eq - Ei)Uq + Ui W(ui, M ), 



h + Xi 



and 



h= f ~u\(r')dr'. 
Jo 



(508) 

(509) 
(510) 

(511) 



The precise relationship of the new potential and its wave functions, which 
are now both BICs, is illustrated in the last column of Table 12.2. 

While the compact form of Eqs. ( |508| - |510| ) explicitly shows the method 
of construction, it is useful to observe that the integral l\ can be conveniently 
re-cast into a simpler form which contains integrals of the form 



u i (r')dr', 



involving the original wave functions only. Making use of eq. 

we get 



(512) 
for u%, 



h 



(E 1 -E )\j + T ^ + 2(E 1 -E ) T ^-U 



(Jo + A) 2 

The second term is integrated by parts as 



117 



o (Jo + A)- 



W 2 {r')dr' 



-W 2 
/o + A 



+ 



(Jo + A) 



r 2WW 

o (To + A) 



dr'. (513) 



dr'. 



(514) 
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We now use eq. ( |500| ) for the derivative of a Wronskian of two solutions of the 
Schrodinger equation to rewrite the second term and observe, that it exactly 
cancels the last term in eq. ( |513| ). We therefore have 

-W 2 (r) 

h{r) = -^- + {E 1 -E f l x {r). (515) 

Here we have made use of the fact that our boundary conditions imply that 
W(0) = 0. 

As an example, we evaluate the two-parameter potential 



V = V 



( W 2 (r) 
In (J + A) pi - E ) 2 h - fAi + Al ] 



(516) 



The argument of the logarithm can be rewritten as 

{Ex - E ) 2 I h - W 2 (r) + AAi + X(Et - E ) 2 h + XJ . (517) 

We happen to have transformed first the state at energy Eq into a BIC and 
then, in the second step, the state at E\, which introduced the parameters 
A and Ai. Let us now consider applying our procedure in the reverse order, 
that is let us first transform the state at energy E\ into a BIC and then the 
state at energy E , producing the parameters /x and /ii. For this situation, 
the argument corresponding to eq. ( |517| ) is 

(E 1 - E ) 2 I h - W 2 (r) + m + ^(E, - E ) 2 I + ^h- (518) 

Clearly, one expects symmetry. This is guaranteed if the parameters are 
related by 

fi = X(Ei - E ) 2 , /n = Ai/(Si - ^o) 2 . (519) 

This also leads to the same two-parameter wave functions. We also note that 
transforming any state twice by eq. (|498|) does not create a second denomi- 
nator or anything else new, but simply changes the value of the parameter A 
as shown in ref . [7t| . Finally, relation (|509| ) ensures that all other eigenstates 
will be solutions to the new potentials. 

Shown in Fig. 12.5 is a potential with two BICs at energies Eq = 1, 
E\ = 4. Clearly, the above procedure can be readily extended to obtain 
multi-parameter families with multiple BICs at arbitrarily selected energies. 
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Our discussion of BICs has been restricted to effectively one dimensional 
problems, and as stated before there is now some experimental evidence for 
the existence of BICs under appropriately chosen conditions ||227| . Recently, 
there has also been a computation which claims the existence of BICs in QED 
in three dimensions [[230|| . These authors have also speculated that the BIC 
energies they compute in e + e~ scattering in QED are in fact in reasonable 
agreement with unusual peaks observed in recent heavy ion experiments ||231J] . 



13 Parasupersymmetric Quantum Mechanics 
and Beyond 

In the last few years, exotic quantum statistics have been widely discussed 
in the literature. For example, in two space dimensions, one can have a one- 
parameter family of statistics interpolating between Bose and Fermi statistics 
[ j232|| . On the other hand, in three and higher space dimensions parafermi and 
parabose statistics ||233| , |234| , |235| , |236|| are the natural extensions of the usual 



Fermi and Bose statistics. In particular, whereas Fermi and Bose statistics 
describe the two one-dimensional representations of the permutation group, 
parafermi and parabose statistics describe the higher dimensional represen- 
tations of the same group. In view of the fact that the SUSY has provided 
us with an elegant symmetry between fermions and bosons, it is natural to 
enquire if there exists a generalization which includes the above exotic statis- 
tics. Such a question was raised many years ago in the context of parastring 
models ||237|| , but the specific symmetry algebra was still the usual SUSY 



one. 

In this section, we study the possibility of having a symmetry between 
bosons and parafermions. We shall construct parasupersymmetric quantum 
mechanics (PSQM) of a boson and a parafermion of order p(= 1,2,3...). It 
turns out that whereas in the usual SUSY QM, the symmetry generators obey 
structure relations that involve bilinear products, in PSQM of order p, the 
structure relations involve products of (p + 1) parasupersymmetry (PARA- 
SUSY) charges. Various consequences of this algebra are also discussed. 

It is worth adding here that historically the PSQM of order 2 was in- 
troduced first ||106|| and its various consequences were discussed ||238| , |239 



240| , [24 1| , |242|| . Initially it was felt that the generalization to order p was not 
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possible in the sense that the PSQM of order p cannot be characterized with 



one universal algebraic relation [243]. However, later on, one was indeed able 



to construct such a PSQM of order p JTU7 



Very recently, new forms of quantum statistics called orthofermi and or- 
thobose statistics, have been constructed fE44j| . It is then natural to construct 



orthosupersymmetric quantum mechanics (OSQM) where there is symmetry 
between a boson and an orthofermion of order p Rl09|| . Unlike the PSQM 
of order p, one finds that the structure relations now involve only bilinear 
products of the symmetry generators. Various consequences of SUSY QM, 
PSQM and OSQM are also discussed in this section. In particular, it is worth 
pointing out that whereas in SUSY QM and OSQM of order p, the energy 
eigenvalues are necessarily nonnegative, in PSQM of order p they need not 
be so. 



13.1 Parasupersymmetric Quantum Mechanics 

In order to motivate the algebra of PSQM, let us recall that in SUSY QM 
the symmetry between a boson and a fermion is characterized by the algebra 

Q 2 = = Q t2 , [H, Q] = 0, QQt + QtQ = 2H. (520) 

Note that there is an extra factor of 2 on the right hand side compared to the 
algebra given in Sec. 2. This results from assuming m — 1 in this section, 
in conformity with the notation followed by the various authors in this field. 
The SUSY QM algebra is easily motivated by recalling that the fermionic 
operators a, a) satisfy the algebra 

a 2 = = a t2 ,{a,a t } = l. (521) 

A useful representation of a and at is in terms of the 2x2 matrices 

1\ + /0 







1 1 o ■ < 522 » 



Let us now consider the parafermi operators b,tf of order p (=1,2,... ) 



which are known to satisfy the algebra |[236|| 

(b) p+l = = (6 f ) p+1 , (523) 
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[[tf,b],b] = -2b,[[b\b],tf] = 2tf. (524) 

On letting 

J + = &t,j_ = &,j 3 = I[&t,&], (525) 

it immediately follows from eqs. ( |523j ) and (|524|) that the operators J± and 
J 3 satisfy the SU(2) algebra 

[J + ,J_] = 2J 3 ,[J 3 ,J ± ]=±J ± . (526) 

Let us now choose J 3 to represent the third component of the spin | 
representation of the SU(2) group as given by 

J 3 = ^(|,|-l,...-| + l,-|). (527) 

It is now easily seen ||107|| that the operators b and W can be represented 



by the following [p + 1) x (p + 1) matrices [a, (3 = 1, 2..., (p + 1)] 

(fyap = Cp6 a ,i3+i; (b j ) a/3 = Cp8 a+Xt p (528) 

where 



C p = ^(3(p-(3 + l) = C p ^ +1 . (529) 

It is easily checked that the operators b and indeed satisfy the algebra 
as given by eqs. ( |523j ) and ( |524p . One can now ask as to what multilinear 
relation is satisfied by b and b^ apart from the one given by eq. (|523|) ? It 
turns out that the nontrivial relation is 

ftPfet + 6 p ~ 1 6 t 6 + ... + bb^bP- 1 + b ] b p = -p(p + l)(p + 2W- 1 , (530) 

6 

where one has (p + 1) terms on the left hand side. As expected, for p = 1 
this reduces to the bilinear relation for the fermionic operators given in eq. 

This relation between b and strongly suggests that one may have an 
analogous multilinear relation in the algebra of PSQM of order p. To that 
purpose, let us choose the PARASUSY charges Qi and Q\ as (p+ 1) x (p+ 1) 
matrices as given by 

(Qi)*p = (P~ iW p )S a , p+1 ; {Q\) a p = {P + iW a )6 a+ltP (531) 
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where a, (3 — l,2,...,(p + l), so that Q\ and Q\ automatically satisfy 

Qi +1 = = (Q\) p+1 - (532) 
Further, it is easily shown that the Hamiltonian (h — m — 1) 

(H) a/ 3 = H a 5 a p, (533) 

where (r = 1, 2, ...,p) 

Hr = ^ + \{W*-W' r ) + l -C T 

H P+ i = y + + + \ C v ( 534 ) 

commutes with the PARASUSY charges Qi and provided (s = 2,3...,p) 

W 2 S _ X + W' s _ x + C s _! = W 2 S - W' s + C s . (535) 

Here C\,Cii ■■■■,C P are arbitrary constants with the dimension of energy. It 
turns out that the nontrivial relation between Q%, Q\ and H is given by ||107 



QiQl + QVqIQx + ■■■ + QiQiQV 1 + Q\Q p i = 2 P Cfi- 1 H, (536) 

and Hermitian-conjugated relations (which we shall not write explicitly) pro- 
vided 

d + C 2 + ... + C P = 0. (537) 

An example is in order at this stage to illustrate the structure of PSQM 
of orde r p. If one chooses 

W 1 = W 2 = ... = W p = ujx, (538) 



then it follows from eq. (|535| ) that in this case (r = 1, 2, p) 

C r+1 -C r = 2u, (539) 
and the Hamiltonian ( |533[ ) takes a very simple form given by 

H = — + -uj 2 x 2 - J 3 u, (540) 
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where J3 is as given by eq. ( |527| ). This H describes the motion of a particle 
with spin p/2 in an oscillator potential and a uniform magnetic field. The 
spectrum of this Hamiltonian is 



E n>m = (n + - - m)uj, (541) 

where n = 0, 1,2, ... and m = p/2, p/2 — 1, —p/2 so that the ground state 
energy of the system is negative unlike in the usual SUSY QM. It is also clear 
from here that whereas the ground state is nondegenerate, the first excited 
state is two-fold degenerate, etc. and finally the p'th and higher excited 
states are (p + l)-fold degenerate. Of course for p = 1 we recover the well 
known results of SUSY QM. 

Several comments are in order at this stage: 



1. For arbitrary W r too one can show that the spectrum of the PSQM 
Hamiltonian is (p + l)-fold degenerate at least starting from the p'th 
and higher excited states. The nature of the ground and the first (p— 1) 
excited states would however depend on the specific form of W r . 

2. With the PSQM Hamiltonian as given by eqs. ( |533| ) and (533), one 
can associate p supersymmetries characterized by the corresponding p 
superpotentials W r {r = 1,2, ...,p) and the corresponding p SUSY QM 
Hamiltonians are given by 

where H r is given by eq. ( |534j ). 

3. Apart from Qi, there are (p — 1) other conserved, independent PARA- 
SUSY charges Q 2 , Q3, Q p defined by (s = 2, 3, ...p) 

Q s = {P-iWp)5 a , p+l if p^s, 

= -{P-iWp)8 a , p+1 if p = s, (543) 



all of which commute with the Hamiltonian ( |534| ) and also satisfy the 
PSQM algebra as given by eqs. @3_D and (|3j) P7|| . 
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4. There is an interesting application of PSQM in nonrelativistic quantum 
mechanics. As we have seen in Sec. 3, given any Hamiltonian Hi with 
p bound states with energies Ei, E 2 , E p and the corresponding eigen- 
functions ?pi,?p 2 , Vv> one can always generate p other Hamiltonians 
H 2 , H 3 , H(p +1 \ with the same spectrum as Hi except that 1,2, ...,p 
levels respectively are missing from them. Further in that case there 
are p supersymmetries with the corresponding Hamiltonians being pre- 
cisely given by eq. ( |542| ) . Besides, the p constants Ci, C 2 , C p are 
related to the energy eigenvalues by (r = 1, 2, ...,p) 

C r = -[E p + E p+ i + ... + E r+ i -{p- l)E r + E r „i + ... + Ei}. (544) 
p 

Thus instead of associating the symmetry algebra s/(l/l) ® SU(2) or 
U{1) ® SU{2) as we did in Sec. 3, one can also associate parasuper- 
symmetry of order p to the hierarchy of Hamiltonians H\, H 2 , H p 



5. What is the most general solution of the relation ( 535 ) which must be 
satisfied in order to have the PSQM of order pi Treating W r ± W r _i 
as the two variables, eq. ( |535| ) can be reduced to a simple nonlinear 



equation which can be immediately solved ||106|| . This is however not 
very useful as it gives us W r + W r ^i for a given W r — W T -\. Ideally, 
we would like to know the most general solution for W 2 for a given Wi 
and then using this W 2 , one would recursively obtain W3, W4, W p . 
Unfortunately, this problem is still unsolved. Of course shape invariant 
potentials satisfy eq. ( |535| ) but shape invariance is clearly not necessary. 
As discussed above, given any H with p bound states, one can always 
construct Wi, W 2 , ....W p which will satisfy eq. (|535|) . 



One unsatisfactory feature of the PSQM of order p is that except for p — 1, 
H cannot be directly expressed in terms of the PARASUSY charges Qi and 
Q\. This is because, in eq. ( |53(j| ), H is multiplied by Q\ whose inverse does 
not exist. Another unsatisfactory feature is that unlike in SUSY QM, the 
energy eigenvalues are not necessarily nonnegative and there is no connection 
between the nonzero (zero) ground state energy and the broken (unbroken) 
PARASUSY. It turns out that in case all the constants Ci, C 2 ....C P are chosen 
to the zero, then one can take care of both of these unsatisfactory features 
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|108|| . In particular, in that case H as given by eq. ( |534| ) can be expressed in 
terms of any one of the p PARASUSY charges Q s by (s = 1, 2, ...,p) 



H = -KQtQs - QsQt) 2 + QtQsQsQ 



)tlV2 



(545) 



It is immediately clear from here that all the energy eigenvalues are nec- 
essarily nonnegative and that the ground state energy being (> 0) corre- 
sponds to unbroken (broken) PARASUSY. Further, one can also show that 
in this case all the excited states are always (p + l)-fold degenerate while the 
nature of the ground state will depend on the specific form of W s 
the limit of p — > oo it has been shown that the relation ( |545|) reduces to 



10S . In 



H 



\q s q 



(546) 



which can be termed as the PSQM of infinite order whose algebra corresponds 
to that of Greenberg's infinite statistics |[245| . 



There is an interesting application of this version of PSQM to the strictly 
isospectral Hamiltonians discussed in Sec. 7. In particular, one can show that 
any potential V\ = Wf — W[ with at least one bound state, forms PSQM of 
order p along with its SUSY partner potential V2 = Wf + W[ and with the 
strictly isospectral potential families V^x, A2), ^2(2;, A3), T4(^, A p ) where 
A 2 , A 3 , X p are arbitrary parameters which are either > or < — 1 [ |2U3| |. 

Eq. (|545| ) shows an explicit expression for the Hamiltonian H which 
involves just one charge Q s but an overall square root. There is an alternative 
expression for H where it is not necessary to take the square root, but which 
involves all the p PARASUSY charges 



10S : 



2H = Q r Ql + QtQ r + - Y^iQlQs + QlQr ~ 2QlQ r 

^ s=l 



(547) 



where, r, s = 1, 2, ...p and s 7^ r . Further, H and Q r also satisfy the simpler 
relation 

Q r QtQ r = 2Q r H. (548) 
In case one chooses W a to be of form (a = l,2...,p) 



A + a - 1 



x 



(549) 
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then one has a model for conformally invariant PSQM of order p. In this case 
one can show that the dilatation operator D and the conformal operator K 
defined by 

D = -^(xP + Px),K = ^x 2 (550) 

satisfy relations analogous to those given in eqs. (|545|) , ( |547| ) and (|536|) in 
terms of the PARASUSY charges Q a and parasuperconformal charges S a 
defined by [a = 1,2, ...,p;i,j = 1,2, (p + l);r = 2,3, 

(S'i)ii = ( 551 ) 



(SV)ij = -x6 i+1) j if i^r + 1, 

= if % — r + 1. (552) 

13.2 Orthosupersymmetric Quantum Mechanics 

Recently, Mishra and Rajasekaran [[244j| have introduced new forms of quan- 
tum statistics called orthofermi and orthobose statistics. Orthofermi statis- 
tics contain a new exclusion principle which is more stringent than the Pauli 
exclusion principle: an orbital state shall not contain more than one particle, 
whatever be the spin direction. The wave function is thus antisymmetric in 
spatial indices alone with the order of the spin indices frozen. In an analo- 
gous way, one can also define orthobose statistics. All these properties follow 
provided the corresponding creation and annihilation operators and C 
satisfy 

v 

Ck a C mf3 ± 5 a p ^2 gamrna Cky = 5k m 5al3i (553) 
7=1 

CkaCmf3 ± C ma Ck/3 = 0, (554) 

where the upper and the lower signs lead to the orthofermi and the ortho- 
bose cases respectively and the Latin indices k,m, ... and the Greek indices 
a,/3,7, ... correspond to space and spin indices respectively. 

For constructing the quantum mechanics of a boson and an orthofermion, 
we ignore the spatial indices in eqs. (|553|) and (|554|) and obtain 



C a C\ + 6 aP £ C\C, = 5 a p, (555) 
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C a C = 0. (556) 

Eq. ( |555Q implies that 

C,C{ = C 2 C\ = ... = C V C\. (557) 

Following the discussion of the last subsection it is easy to see that a 
useful representation of these operators is in terms of the (p + 1) x (p + 1) 
matrices defined by (a, (3 — 1, 2, ...,p; r, s = 1, 2, (p + 1)) 

(C a ) rs = 5 r: i5 Sja+ i, (C^) rs = £ Sl l<5 r ,a+l (558) 

Let us now try to write down the algebra for the OSQM of order p where 
there is a symmetry between a boson and an orthofermion of order p [ITE . 
On comparing the algebra for the fermionic and the orthofermionic opera- 
tors as given by eqs. ( |521|) , (|555|) and (|556|) and remembering the SUSY 
QM algebra as given by eq. (|520| ) it is easy to convince oneself that the p 
ORTHOSUSY charges Q a , Q} a and the Hamiltonian must satisfy the algebra 
(a,/3 = l,2,...,p) 

QaQl + <W E Q\Ql = MapH, (559) 
7=1 

Q a Qp = 0, [H, Q a ) = 0. (560) 

Note that for p = 1 we recover the usual SUSY QM algebra. It is easily 
checked that if we choose the p ORTHOSUSY charges Q a as (p+ 1) x {p+ 1) 
matrices as given by (r = 1, 2, {p + 1)) 

(Q a ) rs = (P - iW a )5 r)1 5 Sja+1 ; (Qi) rs = (P + iW a )S s>1 5 r>a+1 (561) 
and the Hamiltonian H as 

(H) rs = H r 5 rs , (562) 



where (r = 1, 2, ...,p) 



Hi = ^ + \{Wl + W[) 
Hr + i = ^ + \{Wl-W[) (563) 
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then the OSQM algebra as given by eqs. ( |559| ) and ( |560| ) is indeed satisfied 
provided 

Wl + W' a = Wj + W'p. (564) 
Note that this condition directly follows from the OSQM relations 

QiQ\ = Q2QI- = Q P Ql (565) 

which follow from the OSQM relation ( |559| ). Thus, unlike PSQM, the con- 
stants Ci,C2, ... are not allowed in OSQM. The various consequences of the 
OSQM have been discussed in detail in [|109|| . At this point, it may be worth- 
while to make a relative comparison of SUSY QM with PSQM and OSQM. 



1. First of all, the close similarity in structure between OSQM and PSQM 
must be noted. For order p, both are based on (p + 1) x (p + 1) matrices 
and the structure of H in the two cases is very similar. The chief 
difference between the two is the absence of the constants C\, C2..., C p 
in the former while in the latter they may or may not be zero. 

2. Whereas the ground state energy is zero (> 0) in case the SUSY or the 
ORTHOSUSY is unbroken (spontaneously broken), in general, there is 
no such restriction in the PSQM case and the energy eigenvalues can 
even be negative. However, in the special case when Ci,C 2 , ■■■■,C P are 
all zero then the PSQM has similar prediction as the other two. 

3. Whereas in the SUSY QM, all the excited states are necessarily two- 
fold degenerate, in OSQM of order p they are necessarily (p + l)-fold 
degenerate. On the other hand, in the PSQM of order p, all the levels 
starting from the p'th excited state (and above) are necessarily (p+ 1)- 
fold degenerate except when all the constants are zero in which case all 
the excited states are necessarily (p + l)-fold degenerate. 



14 Omitted Topics 

So much work has been done in the area of SUSY QM in the last 12 years 
that it is almost impossible to cover all the topics in such a review. We have 
therefore decided to give a brief description of some of the omitted topics. 
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For each topic, a few references are provided, so that interested readers can 
go back and trace other references and get a good idea of the developments. 

(i) Supermathematics | |246| . 

(ii) SUSY in Atomic Physics H || 



In a series of papers, Kostelecky et al. |41[ have discussed the relationship 
between the physical spectra of different atoms and ions using SUSY QM . In 
particular, they have suggested that the helium and hydrogen spectra come 
from SUSY partner potentials. This connection has been commented upon 
by Rau ||. 

(iii)SUSY in Condensed Matter and Statistical Physics [[L63], [248], |249 
250|, gglj g5^| 



Ideas from SUSY have been used to give insight into random magnetic 
fields in Ising-like models, polymers, electron localization in disordered media 
and ferromagnets. These topics are covered in a set of review talks found in 
the Proceedings of the Conference on Supersymmetry held at CNLS in 1983 
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(iv)Index Theorem and SUSY ||253| , 25 , [2? 

The Atiyah-Singer index theorem can be related to understanding the 
index of the Dirac operator on suitably defined spaces. From our previous 
discussion of SUSY and the Dirac equation, it is not surprising that the index 
of the Dirac operator is related to the Witten index of the related SUSY QM. 
Using techniques similar to calculating the fermion propagator in an external 
field, it has been possible to give a proof of the Atiyah-Singer index theorem 
using the supersymmetric representation of the index theorem. This work 
is best described in lectures of Alvarez-Gaume given at the Bonn Summer 
School p7 |. The question of axial anamolies are best phrased using these 
methods. 



(v) Factorization Method and Solvable Potentials [PH, [13, ET\ 



The method of factorization which can be traced as far back as Bernoulli(1702) 
and Cauchy(1827) can be shown to be exactly equivalent to solving poten- 
tials by SUSY and shape invariance with translation. The reader interested 
in the history as well as an excellent presentation of the factorization method 
and its connection with SUSY is referred to the work of Stahlhofen 



51 



(vi) Group Theory Method and Solvable Potentials |[254 



In the doctoral thesis of Jainshi Wu, it is proven that there is a one to 
one correspondence between using the differential realizations of the S0(2, 1) 
potential group and the factorization method of Infeld and Hull [113]. Thus 
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all the solvable potentials found by Infeld and Hull can be also obtained 
using group theory methods. By extending the symmetry group to 5*0(3, 1) 
and SU (3, 1) the Yale group was also able to study several three dimensional 
scattering problems. 



(vii) Quasi-Solvable Potentials [255 , [256|| 



These potentials, for which a finite number of states can be determined 
analytically, are intermediate between non-solvable and analytically solvable 
potentials discussed in this review. Quasi-solvable potentials have been dis- 
cussed by Shifman and Turbiner [[25 5| , [25611 and using the techniques of SUSY 
QM a large number of new quasi-solvable problems have been discovered 

(viii)SUSY Breaking and Instantons JTJ, gO) |TT], 



One of the least understood problems in supersymmetric field theories is 
that of the origin of SUSY breaking. One suggestion for SUSY breaking is 
that it is of dynamical origin and that instantons are responsible for it. As a 
testing ground of these ideas, the role of instantons in the dynamical breaking 
of SUSY has been studied extensively in various quantum mechanical models. 

(ix) Propagators for SUSY Partner Potentials |)^, [99| 

Since the propagator for a system with a given potential is determined 
from the Hamiltonian, and the various Hamiltonians in a heirarchy of SIP 
are related by 

Q(a s )H^ = H^Q{a s ), (566) 

one can derive recursion relations for the propagators of the heirarchy. If 
one member the heirarchy has a known propagator (such as the harmonic 
oscillator or the free particle) then it is easy to use these recursion relations 
to derive the propagators for the heirarchy. This also allows one to calculate 
path integrals for SIP. 

(x) SUSY and N-body Problem |6(J, |5(J 



A novel algebraic structure is found for 1/r 2 family of many body prob- 
lems which provide a novel link between SUSY and quantum integrability 
p60|| . Also of considerable interest is a supersymmetric generalization of a N- 
particle quantum mechanical model with combined harmonic and repulsive 
forces 1150|| . 



(xi) Coherent states for SIP |[261 



Recently a Lie algebraic treatment of the SIP has been given and using 
it coherent states have been constructed for these potentials. 
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(xii) SUSY and New Soliton Solutions |[262| 

It is well known that the stability equation for the kink solutions in scalar 
field theories in 1+1 dimensions is a Schrodinger-like equation. It has been 
shown that new scalar field theories with kink solutions can be obtained by 
considering the isospectral deformation of this problem. 

(xiii) SWKB and Tunneling p£| || 



Expression for transmission coefficient T through potential barriers has 
been obtained within the SWKB approximation and it has been shown that 
the analytic continuation of T for the inverted SIP (with translation) leads 
to the exact bound state spectrum. 

Time-Dependent Pauli Equation ||264 



xiv 



Kostelecky et al. [ |264) | have succeeded in factorizing the time dependent 
Pauli Equation and utilizing the SUSY to solve the Pauli Equation for a time 
dependent spatially uniform magnetic induction. 

(xv) SUSY and Nuclear Physics |M7], |157|, |265 . 



Relations between the spectra of even-even and neighbouring even-odd 
nuclei has been obtained by using the concepts of SUSY ||247|| . Also, Baye 
[ p. 5711 has shown that the deep and shallow nucleus-nucleus potentials which 
have been successfully used in the past are in fact SUSY partner potentials. 
It has also been suggested that the relationships seen between the energy 
levels of adjacent superdeformed nuclei can be understood in terms of SUSY 



P65 

Acknowledgments It is a pleasure to thank all of our many collaborators 
without whose efforts this work would not have been accomplished. We would 
especially like to thank Joe Ginocchio for carefully reading the manuscript. 
We would like to thank Los Alamos National Laboratory for its hospitality 
during the writ ing of the manuscript. We gratefully acknowledge the U.S. 
Department of Energy for partial support. 



140 



References 

[1] Y.A. Gel'fand and E.P. Likhtman, JETP Lett. 13 (1971) 323. 

[2] P. Ramond, Phys. Rev. D3 (1971) 2415. 

[3] A. Neveu and J. Schwarz, Nucl. Phys. B31 (1971) 86. 

[4] D. Volkov and V. Akulov, Phys. Lett. B46 (1973) 109. 

[5] J. Wess and B. Zumino, Nucl. Phys. B70 (1974) 39; ibid. B78 (1974) 
I. 

[6] For more details and other references in this field see, for example, M.F. 
Sohnius, Phys. Rep. 128 (1985) 39. 

[7] S. Coleman and J. Mandula, Phys. Rev. 159 (1967) 1251. 

[8] S.Ferrara, D. Freedman and P. van Nieuwenhuizen, Phys. Rev. D13 
(1976) 3214. 

[9] S. Deser and B. Zumino, Phys. Lett. B62 (1976) 335. 
[10] E. Witten, Nucl. Phys. B188 (1981) 513. 
[11] F. Cooper and B. Freedman, Ann. Phys. 146 (1983) 262. 
[12] E. Witten, Nucl. Phys. B202 (1982) 253. 

[13] C. Bender, F. Cooper and A. Das, Phys. Rev. D28 (1983) 1473. 

[14] L. Infeld and T.E. Hull, Rev. Mod. Phys. 23 (1951) 21. 

[15] See for example, Supersymmetry in Physics, editors: V.A. Kostelecky 
and D.K. Campbell, North Holland (1985). 

[16] R. Akhoury and A. Comtet, Nucl. Phys. B245 (1984) 213. 

[17] A. Niemi and L.C.R. Wijerwardhana, Phys. Lett. B138 (1984) 389. 

[18] A. Khare and J. Maharana, Phys. Lett. B145 (1984) 77. 

[19] S. Cecotti and L. Girardello, Ann. Phys. 145 (1983) 81. 

141 



[20] D. Boyanovsky and R. Blankenbecler, Phys. Rev. D30 (1984) 1821. 

[21] B. Freedman and F. Cooper, Physica D15 (1985) 138. 

[22] L. Girardello, C. Imbimbo and S. Mukhi, Phys. Lett. B132 (1983) 69. 

[23] A. Kihlberg, P. Salomonson and B.S. Skagerstam, Zeit. Phys. C28 
(1985) 203. 

[24] R.K. Bhaduri, A. Khare and M.V.N. Murthy, Mod. Phys. Lett. A4 
(1989) 2701. 

[25] L. Alvarez-Gaume, Jour. Phys. A16 (1983) 4177. 

[26] L. Alvarez-Gaume and E. Witten, Nucl. Phys. B234 (1984) 269. 

[27] L. Alvarez-Gaume, Bonn Summer School Lectures (1984), unpublished. 

[28] D. Friedan and P. Windey, Physica D15 (1985) 71. 

[29] G. Parisi and N. Sourlas, Nucl. Phys. B206 (1982) 321. 

[30] P. Salomonson and J. van Holten, Nucl. Phys. B196 (1982) 509. 

[31] M. Bernstein and L. Brown, Phys. Rev. Lett. 52 (1984) 1933. 

[32] P. Kumar, M. Ruiz-Altaba and B.S. Thomas, Phys. Rev. Lett. 57 (1986) 
2749. 

[33] F. Marchesoni, P. Sodano and M. Zannetti, Phys. Rev. Lett. 61 (1988) 
1143. 

[34] W.-Y. Keung, E. Kovacs and U. Sukhatme, Phys. Rev. Lett. 60 (1988) 
41. 

[35] M. de Crumbrugghe and V. Rittenberg, Ann. Phys. 151 (1983) 99. 

[36] L.F. Urrutia and E. Hernandez, Phys. Rev. Lett. 51 (1983) 755. 

[37] H. Ui, Prog. Theor. Phys. 72 (1984) 192, 813; H. Ui and A. Takeda, 
ibid. 72 (1984) 266. 



142 



[38] A. Andrianov, N. Borisov and M. Ioffe, Phys. Lett. A105 (1984) 19, 
Phys. Lett. B181 (1986) 141; A. Andrianov, N. Borisov, M. Eides and 
M. Ioffe, Phys. Lett. A109 (1985) 143. 

[39] A. Khare and J. Maharana, Nucl. Phys. B244 (1984) 409. 

[40] A.B. Balantekin, Ann. Phys. 164 (1985) 277. 

[41] A. Kostelecky and M.M. Nieto, Phys. Rev. Lett. 53 (1984) 2285, Phys. 
Rev. A32 (1985) 1293,3243, Phys. Rev. Lett. 56 (1986) 96. 

[42] A. R. P. Rau, Phys. Rev. Lett. 56 (1986) 95. 

[43] C. Blockley and G. Stedman, Euro. Jour. Phys. 6 (1985) 218. 

[44] L. Gendenshtein, JETP Lett. 38 (1983) 356. 

[45] R. Dutt, A. Khare and U. Sukhatme, Phys. Lett. B181 (1986) 295. 

[46] J. Dabrowska, A. Khare and U. Sukhatme, Jour. Phys. A21 (1988) 
L195. 

[47] F. Cooper, J.N. Ginocchio and A. Wipf, Phys. Lett, A129 (1988) 145. 

[48] A. Khare and U. Sukhatme, Jour. Phys. A21 (1988) L501. 

[49] E. Schrodinger, Proc. Roy. Irish Acad. A46 (1940) 9. 

[50] W.T. Reid, Riccati Differential Equations, Math, in Sci. and Eng., Aca- 
demic Press, N.Y. (1972). 

[51] A. Stahlhofen, Duke Univ. preprint (1989). 

[52] F. Cooper, J.N. Ginocchio and A. Khare, Phys. Rev. D36 (1987) 2458. 
[53] C. Chuan, Jour. Phys. A 24 (1991) L1165. 

[54] R. Dutt, A. Khare and U. Sukhatme, Am. Jour. Phys. 56 (1988) 163. 
[55] G. Levai, Jour. Phys. A22 (1989) 689. 

[56] G.A. Natanzon, Vestnik Leningrad Univ. 10 (1971) 22, Teoret. Mat. 
Fiz. 38 (1979) 146. 

143 



[57] J.N. Ginocchio, Ann. Phys. 152 (1984) 203; ibid. 159 (1985) 467. 

[58] A. Khare and U. Sukhatme, Jour. Phys. A26 (1993) L901. 

[59] D. Barclay, R. Dutt, A. Gangopadhyaya, A. Khare, A. Pagnamenta and 
U. Sukhatme, Phys. Rev. A48 (1993) 2786. 

[60] A. Shabat, Inv. Prob. 8 (1992) 303. 

[61] V. Spiridonov, Phys. Rev. Lett. 69 (1992) 298. 

[62] V. Spiridonov, Univ. de Montreal preprint UdeM-LPN-TH-94-92. 

[63] A. Degasperis and A. Shabat, Univ. de Roma preprint No. 919 (1992) 

[64] G. Darboux, C.R. Academy Sc. (Paris) 94 (1882) 1456. 

[65] P. Abraham and H. Moses, Phys. Rev. A22 (1980) 1333. 

[66] D. Pursey, Phys. Rev. D33 (1986) 1098, 1103, 2267; M. Luban and D. 
Pursey, ibid. D33 (1986) 431. 

[67] K. Chadan and P.C. Sabatier, Inverse Problems in Quantum Scattering 
Theory, Springer Verlag (1977). 

[68] M.M. Nieto, Phys. Lett. B145 (1984) 208. 

[69] C.V. Sukumar, Jour. Phys. A18 (1985) 2917, 2937. 

[70] A. Khare and U. Sukhatme, Jour. Phys. A22 (1989) 2847. 

[71] A. Khare and U. Sukhatme, Phys. Rev. A40 (1989) 6185. 

[72] W.-Y. Keung, U. Sukhatme, Q. Wang and T. Imbo, Jour. Phys. A22 
(1989) 987. 

[73] W. Kwong and J. Rosner, Prog. Theor. Phys. Supp. 86 (1986) 366. 

[74] C.V. Sukumar, Jour. Phys. A19 (1986) 2297; ibid. A20 (1987) 2461. 

[75] Q. Wang, U. Sukhatme, W.-Y. Keung and T. Imbo, Mod. Phys. Lett. 
A5 (1990) 525. 



144 



[76] W. Kwong, H. Riggs, J. Rosner and H.B. Thacker, Phys. Rev. D39 
(1989) 1242. 

[77] T. Imbo and U. Sukhatme, Phys. Rev. Lett. 54 (1985) 2184. 

[78] F. Cooper and P. Roy Phys. Lett. A143 (1990) 202. 

[79] A. Comtet, A. Bandrauk and D.K. Campbell, Phys. Lett. B150 (1985) 
159. 

[80] A. Khare, Phys. Lett. B161 (1985) 131. 

[81] For a good review on this field and references therein see, R. Dutt, A. 
Khare and U. Sukhatme, Am. Jour. Phys. 59 (1991) 723. 

[82] R. Dutt, R. De, R. Adhikari and A. Comtet, Phys. Lett. A152 (1991) 
381. 

[83] R. Adhikari, R. Dutt, A. Khare and U. Sukhatme, Phys. Rev. A38 

(1986) 1679. 

[84] K. Raghunathan, M. Seetharaman and S. Vasan, Phys. Lett. B188 

(1987) 351. 

[85] D. Barclay and C.J. Maxwell, Phys. Lett. A157 (1991) 351. 

[86] R. Dutt, A. Khare and Y.P Varshni, Phys. Lett. A123 (1987) 375. 

[87] P. Roy, R. Roychoudhury and Y.P. Varshni, Jour. Phys. A21 (1988) 
1589. 

[88] D. DeLaney and M.M. Nieto, Phys. Lett. B247 (1990) 301. 

[89] Y.P. Varshni, Jour. Phys. A25 (1992) 5761. 

[90] A. Khare and Y.P. Varshni, Phys. Lett. A142 (1989) I. 

[91] S. Fricke, A.B. Balantekin, P. Hatchell and T. Uzer, Phys. Rev. A37 

(1988) 1686. 

[92] Y. Murayama, Phys. Lett. A136 (1989) 455. 



145 



[93] U. Sukhatme and A. Pagnamenta, Phys. Lett. A151 (1990) 7. 

[94] A. Inomata and G. Junker, Lectures on Path Integration, editors: H. 
Cerdeira et al, World Scientific, Singapore (1992); Symposium on Adv. 
Topics in Quantum Mechanics, editors: J.Q. Liang et al, Science Press, 
Beijing (1992). 

[95] R. Dutt, A. Gangopadhyaya, A. Khare, A. Pagnamenta and U. 
Sukhatme, Phys. Lett. A174 (1993) 363. 

[96] R. Dutt, A. Gangopadhyaya, A. Khare, A. Pagnamenta and U. 
Sukhatme, Phys. Rev. A48 (1993) 1845. 

[97] A. Inomata, G. Junker and A. Suparmi, Jour. Phys. A26 (1993) 2261. 

[98] A. Das and W. Huang, Phys. Rev. D41 (1990) 3241. 

[99] A. Jayannavar and A. Khare, Phys. Rev. D47 (1993) 4796. 

[100] B. Thaller, The Dirac Equation, Springer Verlag (1992). 

[101] F. Cooper, A. Khare, R. Musto, and A. Wipf, Ann. Phys. 187 (1988) 
1. 

[102] Y. Nogami and F. Toyama, Phys. Rev. A47 (1993) 1708. 

[103] C.V. Sukumar, Jour. Phys. A18 (1985) L697. 

[104] E. D'Hokcr and L. Vinet, Phys. Lett. B137 (1984) 72. 

[105] H. Yamagishi, Phys. Rev. D29 (1984) 2975. 

[106] V. Rubakov and V. Spiridonov, Mod. Phys. Lett. A3 (1993) 1337. 

[107] A. Khare, Jour. Phys. A25 (1992) L749, Jour. Math. Phys. 34 (1993) 
1277. 

[108] A. Khare, A.K. Mishra and G. Rajasekaran, Mod. Phys. Lett. A8 
(1993) 107. 

[109] A. Khare, A.K. Mishra and G. Rajasekaran, Int. Jour. Mod. Phys. A8 
(1993) 1245. 

146 



[110] A. Jevicki and J. Rodrigues, Phys. Lett. B146 (1989) 55. 

[Ill] M. Shifman, A. Smilga and A. Vainshtein, Nucl. Phys. B249 (1988) 
79. 

[112] J. Casahorran and S. Nam, Int. Jour. Mod. Phys. A6 (1991) 2729; J. 
Casahorran, Phys. Lett. B156 (1991) 1925: P. Roy, R.Roychoudhury 
and Y.P Varshni, Jour. Phys. A21 (1988) 3673. 

[113] P. Roy and R. Roychoudhury, Phys. Rev. D32 (1985) 1597. 

[114] P. Panigrahi and U. Sukhatme, Phys. Lett. A178 (1993) 251. 

[115] J. Pappademos, U. Sukhatme and A. Pagnamenta, Phys. Rev. A48 
(1993) 3525. 

[116] J. von Neumann and E. Wigner, Zeit. Phys. 30 (1929) 465. 

[117] F. Stillinger and D. Herrick, Phys. Rev. All (1975) 446. 

[118] L.E. Ballentine, Quantum Mechanics, Prentice Hall, PP. 205 (1990). 

[119] L.E. Gendenshtein and I.V. Krive, Sov. Phys. Usp. 28 (1985) 695. 

[120] R. Haymaker and A.R.P. Rau, Am. Jour. Phys. 54 (1986) 928. 

[121] D. Lancaster, Nuovo Cimento A79 (1984) 28. 

[122] G. Stedman, Euro. Jour. Phys. 6 (1985) 225. 

[123] A. Lahiri, P. Roy and B. Bagchi, Int. Jour. Mod. Phys. A5 (1990) 1383. 

[124] O. de Lange and R. Raab, Operator Methods in Quantum Mechanics, 
Oxford Univ. Press (1991). 

[125] A. Khare and U. Sukhatme, Physics News (India) 22 (1991) 35. 

[126] L. Landau and E. Lifshitz, Nonrelativistic Quantum Mechanics 3rd 
Volume, Pergamon Press, N.Y. (1977). 

[127] G. Junker, Jour. Phys. A23 (1990) L881. 

[128] R. De, R. Dutt and U. Sukhatme, Jour. Phys. A25 (1992) L843. 

147 



[129] A. Gangopadhyaya, P. Panigrahi and U. Sukhatme, Univ. of Montreal 
preprint UdeM-LPN-TH-93-172 (1993). 



[130 
[131 

[132 
[133 

[134 
[135 
[136 
[137 
[138 
[139 
[140 

[141 
[142 
[143 
[144 
[145 

[146 
[147 



V.A. Kostelecky, M.M. Nieto, D.Truax, Phys. Rev. D32 (1985) 2627. 

A. Khare and R.K. Bhaduri, Inst, of Physics, Bhubaneswar preprint 
IP-BBSR/93-65; |hep-th/ 9310101 . 



A. Khare and R.K. Bhaduri, Jour. Phys. A27 (1994) in press. 

M.M. Nieto, Am. Jour. Phys. 47 (1979) 1067; M.M. Nieto and L.M. 
Simmons,ibid. 47 (1979) 634. 

R. Montemayor and L.D. Salem, Phys. Rev. A40 (1989) 2170 

V. Spiridonov, Mod. Phys. Lett. A7 (1992) 1241. 

L.C. Biedenharn, Jour. Phys. A22 (1989) L873. 

A.M. Macfarlane, Jour. Phys. A22 (1989) 4581. 

G.L. Lamb, Elements of Soliton Theory, Wiley, N.Y. (1980). 

A. Das, Integrable Models, World Scientific, Singapore (1989). 

See for example, P.M. Morse and H. Feshbach, Methods of Theoretical 
Physics, Vol.1, Chapter 5, McGraw Hill, N.Y. (1953). 

F. Calogero, Jour. Math. Phys. 10 (1969) 2191. 

J. Wolfes, Jour. Math. Phys. 15 (1974) 1420. 

B. Sutherland, Jour. Math. Phys. 12 (1971) 246. 
F. Calogero, Jour. Math. Phys. 12 (1971) 419. 

For a list of solvable potentials and other references in the field see 
M. Olshanetsky and A. Perelomov, Phys. Rep. 71 (1981) 314; ibid. 94 
(1983) 6. 

F. Haldane, Phys. Rev. Lett. 60 (1988) 635. 
B. Shastry, Phys. Rev. Lett. 60 (1988) 639. 

148 



[148; 

[149 
[150 
[151 

[152 

[153 

[154 
[155 
[156 
[157 
[158 
[159 

[160 

[161 
[162 

[163 
[164 
[165 
[166 



H. Frahm, Jour. Phys. A26 (1993) L473. 

A. Polychronakos, Phys. Rev. Lett. 69 (1992) 703. 

D. Z. Freedman and P. F. Mende, Nucl. Phys. B344 (1990) 317. 

N. Froman and P. O. Froman, JWKB Approximation, Contributions 
to the Theory, North Holland, Amsterdam (1965). 

J. B. Krieger and C. Rosenzweig, Phys. Rev. 164 (1967) 171; Jour. 
Chem. Phys. 47 (1967) 2942. 

M. Seetharaman and S. Vasan, Jour. Phys. A17 (1984) 2485, 2493; 
ibid. A18 (1985) 1041. 

D. Barclay, A. Khare and U. Sukhatme, Phys. Lett. A183 (1993) 263. 

S. Miller and R. Good, Jr., Phys. Rev. 91 (1953) 174. 

S. Chaturvedi and R. Raghunathan, Jour. Phys. A19 (1986) 1775. 

D. Baye, Phys. Rev. Lett. 58 (1987) 2738. 
R. Amado, Phys. Rev. A37 (1988) 2277. 

C. Gardner, J. Greene, M. Kruskal and R. Miura, Comm. Pure App. 
Math. 27 (1974) 97. 

W. Eckhaus and A. Van Harten, The Inverse Scattering Transforma- 
tion and the Theory of Solitons, North-Holland, Amsterdam (1981). 

S. Chern and C. Peng, Manuscr. Math. 28 (1979) 207. 

F. Calogero and A. Degasperis, Spectral Transforms and Solitons, 
North-Holland, Amsterdam, (1982). 

N. Sourlas, Physica 15D (1985) 115. 

E. Gildner and A. Patrascioiu, Phys. Rev. D16 (1977) 1802. 
S. Ma and G. Mazenko, Phys. Rev. Bll (1975) 4077. 

S. Raby, Unpublished Notes and Lectures, Los Alamos (1981). 

149 



[167] A. Salam and J. Strathdee, Phys. Rev. Dll (1975) 1521. 

[168] E. Gozzi, M. Reuter and W. Thacker, Phys. Lett. A 183 (1993) 29/ 

[169] C. Bender et al, Phys. Rev. D37 (1988) 1472. 

[170] A. Caldeira and A. Leggett, Phys. Rev. Lett. 46 (1981) 211; Ann. Phys. 
(N.Y.) 149 (1983) 374. 

[171] S. Chakravarty, Phys. Rev. Lett. 49 (1982) 681. 

[172] A. Bray and M. Moore, Phys. Rev. Lett. 49 (1982) 1545. 

[173] S. Chakravarty and A. Leggett, Phys. Rev. Lett. 52 (1984) 5. 

[174] G. Schon and A. Zaikin, Phys. Rep. 198 (1990) 237. 

[175] S. Glasstone, K. Laidler and H. Eyring, Theory of Rate Processes, 
McGraw-Hill, N.Y. (1941). 

[176] W. Miller, Jour. Chem. Phys. 61 (1974) 1823; ibid. 62 (1975) 1899. 

[177] N. van Kampen, Jour. Stat. Phys. 17 (1977) 71. 

[178] I. Affleck, Phys. Rev. Lett. 46 (1981) 388. 

[179] D. Boyanovsky, R. Willey and R. Holman, Nucl. Phys. B376 (1992) 
599. 

[180] K. Schonhammer, Zeit. Phys. B78 (1990) 63. 

[181] A. Gangopadhyaya, P. Panigrahi and U. Sukhatme, Phys. Rev. A47 
(1993) 2720. 

[182] S. Coleman, Aspects of Symmetry, Cambridge Univ. Press, N.Y. (1985). 

[183] E. Merzbacher, Quantum Mechanics, Wiley, N.Y. (1970). 

[184] F. Hildebrand, Advanced Calculus for Applications, Prentice-Hall, En- 
glewood Cliffs (1976). 

[185] T. Imbo and U. Sukhatme, Am. Jour. Phys. 52 (1984) 140. 



150 



[186 

[187 

[188 

[189 
[190 
[191 

[192 
[193 

[194 
[195 
[196 
[197 

[198 

[199 
[200 
[201 
[202 
[203 



L. Mlodinow and N. Papanicolaou, Ann. Phys. 128 (1980) 314; ibid. 
131 (1981) 1. 

E. Witten in Recent Developments in Gauge Theories, ed. G.' tHooft 
et al., Plenum, N.Y. (1980). 

C. Bender, L. Mlodinow and N. Papanicolaou, Phys. Rev. A25 (1982) 
1305. 

L. Mlodinow and M. Shatz, Jour. Math. Phys. 25 (1984) 943. 
J. Ader, Phys. Lett. A97 (1983) 178. 

R. Gangopadhyay, G. Ghosh and B. Dutta-Roy, Phys. Rev. A30 (1984) 
594. 

L. Yaffe, Rev. Mod. Phys. 54 (1982) 407. 

M. Sinha-Roy, R. Gangopadhyay and B. Dutta-Roy, Jour. Phys. A17 
(1984) L687. 

J. Miramontes and C. Pajares, Nuovo Cimento 84 (1984) 10. 

S. Hikami and E. Brezin, Jour. Phys. A12 (1979) 759. 

U. Sukhatme and T. Imbo, Phys. Rev. D28 (1983) 418. 

T. Imbo, A. Pagnamenta and U. Sukhatme, Phys. Rev. D29 (1984) 
1669. 

T. Imbo, A. Pagnamenta and U. Sukhatme, Phys. Lett. A105 (1984) 
183. 

T. Imbo and U. Sukhatme, Phys. Rev. D31 (1985) 2655. 

R. Dutt, U. Mukherji and Y.P Varshni, Jour. Phys. B18 (1985) 3311. 

C. S. Lam and Y.P. Varshni, Phys. Rev. A4 (1971) 1875. 

G. Stanciu, Jour. Math. Phys. 8 (1967) 2043. 

A. Khare and C.N. Kumar, Mod. Phys. Lett. A8 (1993) 523. 



151 



[204; 
[205 
[206 
[207 
[208 
[209 
[210 
[211 

[212 
[213 

[214 
[215 

[216 

[217 
[218 
[219 
[220 
[221 
[222 
[223 



Y. Aharonov and A. Casher, Phys. Rev. A19 (1978) 2461. 

R. Jackiw and C. Rebbi, Phys. Rev. D13 (1976) 3358. 

R. Jackiw and J.R. Schriffer, Nucl. Phys. B190 (1981) 253. 

D. K. Campbell and A. Bishop, Nucl. Phys. B200 (1982) 293. 

E. Witten and D. Olive, Phys. Lett. B78 (1978) 97. 

A. d'Adda and P. di Vecchia, Phys. Lett. B73 (1978) 162. 

F. M. Toyama, Y. Nogami and Z. Zhao, Phys. Rev. A47 (1993) 897. 

J. Bjorken and S. Drell, Relativistic Quantum Mechanics, McGraw Hill 
(1964) Chapter 4. 

R. Feynman and M. Gell-Mann, Phys. Rev. 109 (1958) 193. 

L.M. Brown, Phys. Rev. Ill (1958) 193, in Lectures in Theoretical 
Physics, Boulder, Interscience Publishers, Inc., NY. 4 (1962) 324. 

R. Jackiw, Phys. Rev. D29 (1984) 2375. 

R. Musto, L. O'Raifeartaigh and A. Wipf, Phys. Lett. B175 (1986) 
433. 

R. Hughes, V.A. Kostelecky and M.M. Nieto, Phys. Rev. D34 (1986) 
1100, Phys. Lett. B171 (1986) 226. 

K. Fujikawa, Phys. Rev. D21 (1987) 2848. 

S.G. Rajeev, Ann. Phys. (NY.) 173 (1987) 249. 

J. Schwinger, Phys. Rev. 82 (1951) 664. 

T. Imbo and U. Sukhatme, Phys. Rev. D33 (1986) 3147. 

W. Frank, D. Land and R. Spector, Rev. Mod. Phys. 43 (1971) 36. 

H. Narnhofer, Acta Physica Austriaca 40 (1974) 306. 

L. Lathouwers, Jour. Math. Phys. 16 (1975) 1393. 

152 



[224; 
[225 
[226 
[227 
[228; 

[229; 

[230 
[231 

[232 
[233 
[234 
[235 
[236 

[237 
[238 

[239 
[240 



N. Meyer- Vernet, Am. Jour. Phys. 50 (1982) 354. 
H. Friedrich and D. Wintgen, Phys. Rev. A32 (1985) 3231. 
M. Robnik, Jour. Phys. A12 (1979) 1175, ibid. 19 (1986) 3845. 
F. Capasso et al, Nature 358 (1992) 5656. 

Handbook of Mathematical Functions, edited by M. Abramowitz, Irene 
A. Stegun, National Bureau of Standards Applied Math. Series, Nr. 55 
(1964). 

Tables of Coulomb Wave Functions, National Bureau of Standards Ap- 
plied Mathematics Series, Vol. 17 (1952). 

J. Spence and J. Vary, Phys. Lett. B254 (1991) 1. 

T. Cowan et al., Phys. Rev. Lett. 56 (1986) 444; P. Salabura et al, 
Phys. Lett. B245 (1990) 153; W. Koenig et al., Phys. Lett. B218 (1989) 
12. 

J. Leinaas and J. Myrheim, Nuovo Cimento B37 (1977) 1. 

H. Green, Phys. Rev. 90 (1953) 270. 

D. Volkov, Zh. Eksp. Teor. Fiz 38 (1960) 519. 

O. Greenberg and A. Messiah, Phys. Rev. B138 (1965) 1155. 

Y. Ohnuki and S. Kamaefuchi, Quantum Field Theory and Parastatis- 
tics, Tokyo Univ. Press (1982). 

F. Ardalan and F. Mansouri, Phys. Rev. D9 (1974) 334. 

J. Beckers and N. Debergh, Mod. Phys. Lett. A4 (1989) 1209, 2289; 
Jour. Math. Phys. A31 (1990) 1523 ; Jour. Phys. A23 (1990) L751, 
L1073 ; Nucl. Phys. B340 (1990) 777. 

S. Durand, R. Floreanini, M. Mayrand and L. Vinet, Phys. Lett. B233 
(1989) 158 ; S. Durand and L. Vinet, Jour. Phys. A23 (1990) 3661. 

V. Spiridonov, Jour. Phys. A24 (1991) L529. 



153 



[241] A. Andrianov and M. Ioffe, Phys. Lett. B255 (1991) 543 ; A. Andri- 
anov, M. Ioffe, V. Spiridonov and L. Vinet, ibid. B272 (1991) 297. 

[242] V. Merkel, Mod. Phys. Lett. A5 (1990) 2555; ibid. A6 (1991) 199. 

[243] S. Durand, M. Mayrand, V. Spiridonov and L. Vinet, Mod. Phys. Lett. 
A6 (1991) 3163. 

[244] A.K. Mishra and G. Rajasekaran, Pramana (Jour. Phys.) 36 (1991) 
537; ibid. 37 (1991) 455(E). 

[245] O. Greenberg, Phys. Rev. Lett. 64 (1990) 705. 

[246] F. Berezin, The Method of Second Quantization, Nauka, Moscow 
(1965). 

[247] F. Iachello, Phys. Rev. Lett. 44 (1980) 772; Physica D15 (1985) 85 
and references therein; R. Casten, Physica D15 (1985) 99. 

[248] J. Cardy, Physica D15 (1985) 123 and references therein. 

[249] Y. Shapir, Physica D15 (1985) 129 and references therein. 

[250] Y. Nambu, Physica D15 (1985) 147 and references therein. 

[251] K. Efetov, Sov. Phy. JETP 55 (1985) 514; Adv. Phys. 32 (1983) 53. 

[252] E. Brezin, D. Gross and C. Itzykson, Nucl. Phys. B235 (1984) 24. 

[253] M. Atiyah and I.M. Singer, Bull. Am. Math. Soc. 69 (1963) 422. 

[254] Y. Alhassid, F. Gursey and F. Iachello, Phys. Rev. Lett. 50 (1983) 873; 
Ann. of Phys. 148 (1983) 336; J. Wu, Ph.D. thesis, Yale Univ. (1985), 
unpublished. 

[255] M. Shifman, Int. Jour. Mod. Phys. A4 (1989) 2897. 

[256] A. Turbiner and A. Usheridze, Phys. Lett. A126 (1987) 181, A. Tur- 
biner, Sov. Phys. JETP 67 (1988) 230 ; Comm. Math. Phys. 118 (1988) 
467. 

[257] D.P Jatkar, C.N. Kumar and A. Khare, Phys. Lett. A142 (1989) 200. 



154 



[258] R.B. Abbott and W.J. Zakrzewski, Z. Phys. C20 (1983) 227. 

[259] A. Khare and J. Maharana, Zeit. Phys. C3 (1984) 191. 

[260] B.S. Shastry and B. Sutherland, Phys. Rev. Lett. 70 (1993) 4029, B.S. 
Shastry, ibid. 69 (1992) 164. 

[261] T. Fukui and N. Aizawa, Phys. Lett. A180 (1993) 308; T. Fukui, 
Yukawa Institute, Kyoto preprint (1993). 

[262] C.N. Kumar, Jour. Phys. A20 (1989) 5347. 

[263] T. Sil, A. Mukherjee, R. Dutt and Y.P Varshni, Phys. Lett. A184 
(1993) in press. 

[264] V.A. Kostelecky, V.I. Man'ko, M.M. Nieto, and D. Truax, Phys. Rev. 
A48 (1993) 951. 

[265] R. D. Amado, R. Bijker, F. Cannata and J. Dedonder, Phys. Rev. Lett. 
67 (1991) 2777. 



155 



FIGURE CAPTIONS 



Figure 2.1: The energy levels of two supersymmetric partner potentials V\ 
and Vi- The figure corresponds to unbroken SUSY. The energy levels are 
degenerate except that V\ has an extra state at zero energy Eq = 0. The 
action of the operators A and A* in connecting eigenfunctions is shown. 

Figure 2.2: The infinite square well potential V — of width L = n and 
its supersymmetric partner potential 2cosec 2 x in units H = 2m = 1. The 
ground state of the infinite square well has energy 1. Note the degenerate 
higher energy levels at energies 2 2 ,3 2 ,4 2 , ... 

Figure 4.1: Self-similar superpotentials W(x) for various values of the de- 
formation parameter q. The curve labeled H.O. (harmonic oscillator) cor- 
responds to the limiting case of q = 1. Note that only the range x > is 
plotted since the superpotentials are antisymmetric W(x) = — W(— x). 

Figure 4.2: Self-similar potentials V\{x) (symmetric about x = 0) corre- 
sponding to the superpotentials shown in Fig. 4.1 . The curves are taken 
from ref. f59|. 



Figure 4.3: A double well potential Vi(x) (solid line) and its single well 
supersymmetric partner V%{x) (dotted line). Note that these two potentials 
are shape invariant [^j with a scaling change of parameters. The energy 



levels of V\{x) are clearly marked. 

Figure 5.1: Diagram showing how all the shape invariant potentials of Table 
4.1 are inter- related by point canonical coordinate transformations. Poten- 
tials on the outer hexagon have eigenfunctions which are hypergeometric 
functions whereas those on the inner triangle have eigenfunctions which are 
confluent hypergeometric functions. These two types are related by suitable 



limiting procedures. The diagram is from ref. ||129| 



Figure 7.1: Selected members of the family of potentials with energy spectra 
identical to the one dimensional harmonic oscillator with u = 2. The choice 
of units is h = 2m = 1. The curves are labeled by the value Ai, and cover 
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the range < Ai < oo. The curve Ai = oo is the one dimensional harmonic 
oscillator. The curve marked Ai = is known as the Pursey potential []66 
and has one bound state less than the oscillator. 

Figure 7.2: Ground state wave functions for all the potentials shown in Fig. 
7.1, except the Pursey potential. 

Figure 7.3: A schematic diagram showing how SUSY transformations are 
used for deleting the two lowest states of a potential V\{x) and then re- 
inserting them, thus producing a two-parameter (Ai, A2) family of potentials 
isospectral to V\{x). 

Figure 7.4: The pure three-soliton solution of the KdV equation as a function 
of position (x) and time (t). This solution results from constructing the 
isospectral potential family starting from a reflectionless, symmetric potential 
with bound states at energies E\ = —25/16, E 2 = —1, E 3 = —16/25. Further 
details are given in the text and ref. [[75]]. 

Figure 9.1: A "deep" symmetric double well potential Vi(x) with minima at 
x = ±xo and its supersymmetric partner potential V^x). 

Figure 9.2: Supersymmetric partner potentials V\{x) and Vzfa) correspond- 
ing to two choices of the parameter xq for the potentials given in eq. ( |392| ). 
For a detailed discussion, see the text and ref. Q . 



Figure 9.3: The energy splitting t = E\ — E as a function of the separation 
2x Q of the superposed Gaussians in the ground state wave function ipo(x). 



This figure, taken from ref. f34|| , shows the remarkable accuracy of the SUSY- 



based energy splitting computations for a double well potential. 

Figure 12.1: The harmonic oscillator potential V\(i) and its partner potential 
^2(1) as given by eq. ([48 8|) . These potentials arise from a singular super- 
potential W\ = x — 1/x corresponding to the choice <fi = ipi(x), e — E\ in 
Table 12.1 . The energy levels for both potentials are shown. Notice the 
negative energy state of Vim at energy —2, and the partial degeneracy of the 
eigenvalue spectra coming from the even parity of the potentials ||114 . 
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Figure 12.2: The Morse and its supersymmetric partner potential [eq. ( |494| )1 
coming from the singular superpotential W± given in eq. ( 493 ). The singu- 
larity at x = —In 3.5 breaks V 2 (i) into two disjoint pieces, and the eigenvalue 



spectra have no degeneracy [114 



Figure 12.3: Potentials V(r) (solid lines) possessing one bound state in the 
continuum (BIC) obtained by starting from a free particle and constructing 
a one-parameter (A) family. The BIC wave function Uo(r) is at energy E — 1, 
and is shown by the dashed line. Fig. (a) corresponds to A = 0.5 and fig. 
(b) corresponds to a much larger value A = 5 . 

Figure 12.4: (a) A potential (solid line) with one BIC at energy E = 0.25 
obtained by starting from an attractive Coulomb potential (dotte d line), (b) 
A plot of the wave functions corresponding to the potentials in part (a). 

Figure 12.5: (a) An example of a potential with two BICs at energies E = 1 
and Ei = 4. (b) The wave functions at energies E = 1 (dashed line) and 
E 1 = 4 (dotted line) |TJ . 



TABLE CAPTIONS 



Table 4.1: All known shape invariant potentials in which the parameters a 2 
and ai are related by a translation (a 2 = a>i + a). The energy eigenvalues and 
eigenfunctions are given in units h = 2m = 1. The constants A, B, a, u>, I are 
all taken > 0. Unless otherwise stated, the range of potentials is — oo < x < 
oo, < r < oo. For spherically symmetric potentials, the full wave function 
is ip n im{j' 1 6 1 (()) = ipni{ r )Yim{9 ,</>)• Note that the wave functions for the first 
four potentials (Hermite and Laguerre polynomials) are special cases of the 
confluent hypergeometric function while the rest (Jacobi polynomials) are 
special cases of the hypergeometric function. Fig. 5.1, taken from ref. | ]1 29j ] , 
shows the inter-relations between all the SIPs in the table via point canonical 
coordinate transformations. 
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Table 6.1: Comparison of the lowest order WKB and SWKB predictions 
for the bound state spectrum of the Ginocchio potential for different values 
of the parameters A, v and several values of the quantum number n. The 
exact answer is also given. Units corresponding to % = 2m = 1 are used 
throughout. 

Table 9.1: Comparison of the three lowest energy eigenvalues obtained by a 
variational method with the exact results. 



Table 12.1: Choice of the constant e and the corresponding solution <fi{x) in 
eq. (|480|) determines the choice of the superpotential which produces any 
given non-singular potential V(x) [eq. ( 481 )]. The table shows how certain 
choices of e give rise to singular superpotentials, negative energy eigenstates 
and a breakdown of the degeneracy theorem. We have taken <p(x = — oo) = 
for convenience [|114||. 



Table 12.2: One and two parameter families of potentials with bound states 
in the continuum. These families are generated by applying the theorem 
described in the text to scattering states Uq and U\ at energies E and Ex in 
succession 11115 
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